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SUMMARY 


The  extreme  sensitivity  of  superconducting  magnetic  gradio- 
meters  affords  means  of  measuring  gradients  of  fluctuating  magnetic 
fields  generated  by  ocean  waves.  The  fields  fluctuate  in  concert 
with  motions  of  ocean  waves,  so  that  reading  fluctuations  of 
magnetic  fields  above  an  ocean  tells  the  motion  of  seawater  below 
its  surface.  As  a first  step  in  developing  requisite  experimental 
and  analytical  techniques,  we  plan  to  use  a superconducting  mag- 
netic gradiometer  to  measure  fluctuating  gradients  of  magnetic 
fields  generated  at  a fixed  point  above  the  surface  by  waves 
passing  the  oceanographic  tower  operated  by  the  Naval  Undersea 
Center . 

The  response  of  a superconducting  magnetic  gradiometer  to 
ambient  magnetic  fields  having  gradients  that  are  sensibly  con- 
stant over  the  length  of  its  axis  (typically  25  cm  or  so)  is 
a sum  of  two  terms.  One  term  is  proportional  to  the  magnetic 
field  at  the  midpoint  of  the  axis,  and  the  other  term,  to  a 
gradient  of  the  magnetic  field  at  the  midpoint.  The  term  pro- 
portional to  the  magnetic  field  results  from  slight  differences 
in  area  and  orientation  of  the  loops  forming  the  pickup  cir- 
cuit of  the  gradiometer  and  vanishes  for  a perfectly  balanced 
gradiometer.  The  term  proportional  to  a gradient  of  the  mag- 
netic field  depends  on  five  independent  elements  of  a symmetric 
matrix  having  a vanishing  trace  that  represents  gradients  of 
the  ambient  field  at  the  midpoint. 

The  five  elements  define  location,  orientation,  and  moment 
of  a dipole  positioned  on  a sphere  of  unit  radius  about  the 
midpoint  that  gives  gradients  at  the  midpoint  equal  to  gradients 
of  the  ambient  field  there.  The  construct  of  an  equivalent 
dipolo  concisely  describes  the  gradient  response.  For  example, 
fluctuations  in  location,  orientation,  and  moment  of  an  equi- 
valent dipole  describe  gradientti  from  fluctuating  magnetization 
currents  in  magnetic  objects  as  well  as  from  fluctuations  in 
position  and  orientation  of  a gradiometer. 

★ 

The  tower  is  located  about  one  mile  off  Mission  Beach  near 
San  Diego,  California,  in  18  ra  of  water. 
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Moreover,  it  provides  a rationale  for  two  techniques  that 
afford  means  of  suppressing  gradient  fluctuations  from  nearby 
magnetic  objects;  namely,  (1)  using  a gradiometer  axis  having 
coplanar  pickup  loops  (transverse  configuration)  and  aligning 
it  for  a null  response  and  (2)  using  a magnetic  dipole  to 
cancel  steady  gradients.  Transverse  gradients  of  the  magnetic 
field  of  a dipole  vanish  in  certain  directions,  so  that  aligning 
a gradiometer  measuring  transverse  gradients  for  a null  response 
suppresses  fluctuating  gradients  from  magnetic  objects.  Placing 
an  opposing  dipole  at  the  location  of  an  equivalent  dipole 
nullifies  gradients  at  the  field  point  and  so  suppresses  appar- 
ent gradient  fluctuations  owing  to  changes  in  gradiometer  orien- 
tation . 

Preliminary  tests  of  suppression  techniques  conducted  v;iti' 
our  gradiometer,  which  has  two  coplanar  piclcup  loops  spaced 
25  cm  between  centers,  demonstrate  that  their  effectiveness  in 
limited  only  by  precision  of  requisite  alignments.  Aligning 
the  gradiometer  within  a few  degrees  ('v  0.1  radian)  for  a null 
response  to  gradient  fluctuations  from  a spherical  iron  shell 
having  a radius  of  0.74  m and  placed  4.5  m north  of  the  gradio- 
meter reduces  RMS  fluctuations  in  the  frequency  band  0.001  to 
0.01  Hz  by  about  a factor  of  ten.  Crude  positioning  of  a coil 
approximating  an  opposing  magnetic  dipole  at  the  location  of 
an  equivalent  dipole  reduces  steady  gradients  from  the  sphere 
by  a factor  of  100. 

Design  of  the  gradiometer,  insofar  as  practical,  eliminates 
sources  of  instrument  noise  driven  by  fluctuations  of  ambient 
temperature  and  pressure  and  of  the  earth's  magnetic  field  by 
using  materials  within  the  helium  bath  that  are  free  of  para- 
magnetic impurities.  A rectangular  bloc)c  of  high-purity  crys- 
talline silicon,  for  example,  supports  the  two  niobium  wire 

picJcup  loops  on  a lateral  face.  Spectra  that  characterize 
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noise  in  the  frequency  range  5 X 10  to  20  Hz  of  the  gradio- 
meter operating  in  a magnetically  quiet  environment  are  effec- 
tively white  at  frequencies  above  about  0.1  Hz  with  a spectral 
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density*  of  0.03  (pT/m)^/Hz  and  are  of  the  form  S(f)  = S(f^) 
(f^/f)^  at  frequencies  below  0.1  Hz.  A value  of  1.3  for  y 
with  S(f  ) = 0.002  (pT/m)^/Hz  at  f^  = 1 Hz  is  representative. 
Intrinsic  noise  of  the  superconducting  quant  ujm  interference 
device  (SQUID)  used  in  the  instrument  to  sense  magnetic  flux 
dominates  spectra  observed  in  a quiet  environment  at  freqencie 
below  0.1  Hz. 

As  shown  in  the  following  figure,  spectral  density  of 
inherent  noise  of  the  instrument  is  well  below  the  level  of 
spectral  densities  of  gradients  expected  from  oceanic  internal 
waves  passing  the  oceanographic  tower  operated  by  the  Naval 
Undersea  Center.  To  measure  gradients  from  waves  passing  the 
tower,  we  plan  to  jut  the  gradiometer  over  water  on  a rigid 
nonmagnetic  cantilever  extending  horizontally  25  meters  from 
the  centerline  of  the  tower  off  its  south  face.  Measurements 
of  steady  gradients  from  magnetization  currents  in  the  steel 
structure  of  the  tower  show  that  gradients  at  positions  in  a 
vertical  plane  of  symmetry  of  the  structure  are  represented 
by  gradients  of  equivalent  dipoles  located  in  the  plane  of 
symmetry.  Aligning  the  gradiometer  with  a precision  of  0.01 
radians,  or  0.6  degrees,  for  a null  response  from  the  equiva- 
lent dipole  Iceeps  noise  from  fluctuating  magnetization  current 
in  the  tower  below  the  level  of  instrument  noise.  Structural 
design  of  the  cantilever  limits  translational  and  rotational 
RI4S  fluctuations  of  the  instrument  to  0.1  mm  and  a few  seconds 
of  arc  ( 10~^  radians)  in  the  bandwidth  0.002  to  0.005  Hz  and, 
together  with  a coil  positioned  on  the  cantilever  to  reduce 
steady  gradients  from  the  tower,  Iceeps  noise  from  irregular 
motions  of  the  instrument  in  the  steady  gradients  below  the 
level  of  instrument  noise. 


*We  use  MKS  units  throughout  and  conform  to  the  international 
convention  suggesting  use  of  the  unit  Tesla,  which  is  a 

Weber/m^  and  so  1 pT/m  = 10~^^T/m  = 10~^°G.cm  - 10"^gamma/m. 
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f,coh 

iD 


f,  Hz 

Three  typical  spectra  of  fluctuating  gradients  of  magnetic 
fields  expected  7 m above  the  surface  from  internal  waves  passing 
the  oceanographic  tower  operated  by  the  Naval  Undersea  Center  as 
compared  to  the  inherent  noise  spectrum  of  the  instrument.  Spectra 
are  estimates  based  on  spectral  measurements  of  isotherm  displace- 
ments made  during  August  1972.  Serial  numbers  marking  spectra  tell 
the  month,  day,  and  local  time  at  the  beginning  of  corresponding 
time  series;  for  example,  08211540  means  August  21  at  15:40  hours 
(PST)  . 
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Section  I 
INTRODUCTION 

The  unprecedented  sensitivity  of  superconducting  magnetic 
gradiometers  affords  means  of  measuring  fluctuating  gradients 
of  magnetic  fields  generated  by  oceai  waves.  Seawater  moving 
across  the  earth's  magnetic  field  drives  electric  currents 
that  produce  weak  magnetic  fields  above  the  oceans.  The  fields 
fluctuate  in  concert  with  motions  of  ocean  waves,  so  that 
reading  fluctuations  of  magnetic  fields  above  an  ocean  tells 
the  motion  of  seawater  below  its  surface.  Podney  (1975)  among 
others  describes  the  fluctuating  magnetic  field  that  results 
from  a wave  progressing  horizontally  in  a stratified  ocean. 

Our  aim  is  to  tell  motions  of  seawater  from  measurements 
of  fluctuating  magnetic  fields  over  oceans.  As  a first  step 
in  developing  requisite  experimental  and  analytical  techniques, 
we  plan  to  use  a superconducting  magnetic  gradiometer  to  mea- 
sure fluctuating  gradients  of  magnetic  fields  generated  at  a 
fixed  point  above  the  surface  by  waves  passing  the  oceanographic 
tower  operated  by  the  Naval  Undersea  Center.  The  tower  is  loca- 
ted about  one  mile  off  Mission  Beach  near  San  Diego,  California 
in  18  m of  water.  Successful  completion  of  the  first  step  pro- 
vides the  experience  and  information  necessary  for  planning, 
carrying  out,  and  interpreting  measurements  over  the  open  ocean 
using  an  airborne  gradiometer.  Measurements  at  a fixed  point 
above  the  surface,  in  effect,  provide  the  code  for  reading 
motions  of  seawater  in  Measurements  made  over  an  open  ocean 
from  aircraft. 

The  tower  provides  a stable  platform  for  jutting  the  instru 
ment  over  water  on  a cantilever  support.  Surface  and  internal 
waves  passing  the  tower  have  simple  spectra  that  are  well 
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characterized  by  past  measurements  using  thermistor  cnains  and 
current  meters  (Cox,  1962;  La  Fond,  1962;  Winant  and  Olson, 
1976;  and  Ziemer , 1976)  and  provide  a readily  accessible  source 
for  first  measurements.  Nonetheless,  iragnetizat ion  currents 
in  the  steel  structure  of  the  tower  malce  the  task  of  measuring 
fluctuating  gradients  from  the  waves  an  exacting  one.  To  pre- 
pare for  carrying  out  measurements  at  the  tower,  ve  completed 
the  following  tasks: 

o Specification,  procurement,  and  acceptance  tests  of 
a suitable  superconducting  magnetic  gradiometer* 

o Land  based  trials  of  the  gradiometer  conducted  to 
estimate  spectra  of  inherent  instrument  noise  in  a 
quiet  environment  and  to  determine  effectiveness  of 
techniques  designed  to  suppress  noise  from  nearby 
magnetic  objects 

o Specification  of  fluctuating  radients  of  magnetic 
fields  expected  from  waves  passing  the  tower. 

o Use  of  a fluxgate  gradiometer  to  measure  gradients 
of  the  steady  magnetic  field  produced  by  magnetiza- 
tion currents  in  the  steel  structure  of  the  tower 
(Gillespie  and  Podney,  1976) 

o Design  of  a rigid  nonmagnetic  cantilever  re- 
quired to  project  the  instrument  over  water** 

Ziemer  and  Gillespie  (1976)  describe  early  planning  for 
measurements  at  the  tower  and  present  a chronology  of  the 
design  specification,  procurement,  and  acceptance  tests  of 
the  gradiometer.  Here,  we  describe  use  of  the  instrument 

* Superconducting  Technology,  Inc.  constructed  the  instrument. 

**  Mechanics  Research,  Inc.  developed  a structural  design  to 
meet  stability  requirements  set  by  Physical  Dynamics,  Inc. 
(Haire  and  Van  Lerberg,  1976) 
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to  measure  gradients  of  ambient  magnetic  fields  and  develop  a 
procedure  for  operating  at  the  tower  that  suppresses  noise 
from  magnetization  currents  in  the  steel  structure.  The  pro- 
cedure unites  results  of  the  last  three  tasks. 

We  first  describe  design  features  of  the  instrument  in 
Section  II.  Insofar  as  practical,  the  aesigri  eliminates  sources 
of  instrument  noise  driven  by  fluctuations  of  ambient  tempera- 
ture and  pressure  and  cf  the  earth's  magnetic  field  by  using 
materials  within  the  helium  bath  that  are  free  of  paramagnetic 
impurities.  Two  loops  of  niobium  wire  mounted  on  a lateral 
face  of  a rectangular  block  (2  X 2 X 12  in.)  of  high-purity 
crystalline  silicon,  for  example,  form  the  pickup  circuit  of 
the  gradiometer. 

Section  III  gives  a formulation  describing  instrument 
response  to  ambient  magnetic  gradients  that  are  sensibly  con- 
stant over  the  distance  separating  centers  of  the  pickup  loops. 
Because  the  loops  are  coplanar,  the  gradiometer  responds  to 
a transverse  gradient  of  the  magnetic  field  at  the  point  mid- 
way between  centers  of  the  loops.  Slight  differences  in  area 
and  orientation  of  the  loops  give  an  imbalance  response  pro- 
portional to  the  magnetic  field  at  the  midpoint  as  well.  The 
gradient  field  of  a magnetic  dipole  located  on  a sphere  of 
unit  radius  about  the  midpoint  describes  the  part  of  the  re- 
sponse owing  to  ambient  gradients. 

Section  IV  describes  procedures  for  making  mechanical  and 
electronic  adjustments  to  nullify  the  part  of  the  response  owing 
to  imbalance.  Moving  three  small  niobium  disks  mounted  on 
slides  near  one  loop  adjusts  its  effective  area  and  orientation 
and  balances  response  of  the  loops  to  a uniform  magnetic  flux. 
Electronically  substracting  a fraction  of  the  magnetic  field 
components  measured  by  a triaxial  fluxgate  magnetometer  near 
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the  midpoint  from  the  qradiometer  response  provides  r iier 
balancing.  Procedures  for  balancing  a gradiometer  in  both 
uniform  and  nonuniform  ambient  magnetic  fields  are  presented. 

In  Section  V,  we  present  the  first  measurements  of  spectra 
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that  characterize  noise  in  the  frequency  range  5 X 10  to  20 
Hz  of  a superconducting  magnetic  gradiometer  operating  in  a 
magnetically  quiet  environment.  Spectra  are  effectively  white 
at  frequencies  above  about  0.1  Hz  with  a spectral  density*  of 
0.03  (pr/m)^/Hz  and  are  of  the  form  S(f)  = S(f^)(f^/f)^  at 
frequencies  below  0.1  Hz.  A value  of  1.3  for  y with  S(f^) 

0.002  (pT/m)^/Hz  at  f^  = 1 Hz  is  representative  of  values 
observed  during  quiet  periods.  Intrinsic  noise  of  the  super- 
conducting quantum  interference  device  (SQUID)  used  in  :he 
instrument  to  sense  magnetic  flux  then  dominates  observed 
spectra  at  frequencies  below  0.1  Hz.  Use  of  a dc-type  SQUID 
(Clarke  et.  al. , 1975)  would  improve  instrument  performance 

nearly  tenfold. 

In  Section  VI,  we  give  a rationale  for  two  techniques 
that  provide  means  of  suppressing  noise  from  nearby  magnetic 
objects;  namely,  (1)  aliening  the  gradiometer  to  obtain  a null 
response,  and  (2)  using  a magnetic  dipole  to  cancel  steady  gra- 
dients. We  describe  noise  from  nearby  magnetic  objects  in  terms 
of  gradient  fluctuations  corresponding  to  small  changes  in  loca- 
tion, orientation,  and  moment  of  a magnetic  dipole  located  on 
a sphere  of  unit  radius  about  the  midpoint  of  the  gradiometer 
axis.  Because  transverse  gradients  of  the  magnetic  field  of 
a dipole  vanish  in  certain  directions,  response  of  a gradio- 
meter measuring  transverse  gradients  vanishes  at  certain  orien- 
tations. Aligning  it  for  null  response  suppresses  noise  from 
magnetic  objects.  Placing  an  oppositely  directed  dipole  at  the 
location  of  an  equivalent  dipole  nullifies  steady  gradients  at 
the  midpoint  and  so  suppresses  apparent  gradient  fluctuations 
owing  to  changes  in  gradiometer  orientation. 

* We  use  MKS  units  throughout  and  conform  to  the  international 
convention  suggesting  use  of  the  unit  Tesla,  which  is  a 
Weber/m^  and  so  IpT/m  = 10"^^T/m  = 10"^°G/cm  = 10  gamma/m. 
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Section  VII  presents  results  of  preliminary  tests  of  noise 
suppression  techniques.  Our  field  tests  show  that  the  construct 
of  an  equivalent  dipole  gives  a useful  description  of  the  re- 
sponse to  gradients  that  are  sensibly  constant  over  the  length 
of  a gradiometer  axis.  Preliminary  tests  using  crude  means 
of  aligning  the  gradiometer  reduce  noise  power  by  a factor  of 
lOh  and  suggest  that  reduction  of  noise  is  limited  by  align- 
merit  precision  alone. 

Finally,  Section  VIII  describes  a procedure  for  operating 
at  the  tower  that  both  uses  the  techniques  to  suppress  noise 
from  magnetization  currents  in  the  steel  structure  and  gives 
a maximum  response  to  gradients  from  internal  waves.  We  plan 
to  jut  the  gradiometer  over  water  on  a rigid  nonmagnetic  canti- 
lever extending  horizontally  25  meters  from  the  centerline  of 
the  tower  off  its  south  face.  Figure  1 shows  that  spectral 
densities  of  gradients  expected  from  internal  waves  are  well 
above  the  level  of  instrument  noise.  Aligning  the  gradiometer 
for  a null  response  with  a precision  of  0.01  radians,  or  0.6 
degrees,  keeps  noise  from  fluctuating  magnetization  currents 
in  the  tower  below  the  level  of  instrument  noise.  Limiting 

translational  and  rota'-ional  RMS  fluctuations  of  the  instrument 
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to  0.1  mm  and  a few  seconds  of  arc  ( 10  radians)  in  the  band- 
width 0.002  to  0.005  Hz  keeps  noise  from  irregular  motions  of 
the  instrument  in  the  steady  gradient  field  of  the  tower  below 
the  level  of.  instrument  noise. 
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Section  II 

INSTRUMENT  DESCRIPTION 

As  shown  in  Figure  2a,  the  instrument  comprises  a gradio- 
meter  probe  immersed  in  a bath  of  liquid  helium  contained 
in  the  interior  vessel  of  a Dewar  made  of  a fiber  glass 
laminate.  A wrapping  of  about  50  alternating  layers  of 
fiber  glass  cloth  and  aluminized  Mylar  covering  the  interior 
vessel  reduces  radiative  heat  flow  across  the  vacuum  separa- 
ting the  inner  vessel  and  outer  jacket  of  the  Dewar.  The 
outer  jacket  is  26"  in  diameter  and  4C"  in  length,  and  the 
Dewar  weighs  about  250  lbs.  Liquid  helium  in  the  reservoir, 
Vvhich  has  a capacity  of  160  liters,  boils  off  to  the  atmo- 
sphere at  a rate  somewhat  less  than  3 liters  per  day. 

The  probe  supports  a rectangular  block  (2"  X 2"  X 12") 
of  high— purity,  crystalline  silicon.  Two  loops  of  niobium 
wire  mounted  on  a lateral  face  of  the  block  form  the  pickup 
circuit  of  the  gradiometer.  Three  small  niobium  disks 
mounted  on  slides  near  the  top  loop  provide  means  of  adjusting 
the  effective  area  and  orientation  of  the  loop  in  order  to 
balance  response  of  the  loops  to  a uniform  magnetic  flux. 

A niobium  capsule  fixed  midway  between  loop  centers  contains 
a torodial,  point-contact  type  SQUID  that  senses  differential 
magnetic  flux  threading  the  pickup  loops. 

The  silicon  block  is  housed  at  the  foot  of  the  probe  in 
a tube  3"  in  diameter  and  15"  long,  which  fits  into  a 
matching  socket  fixed  to  the  bottom  of  the  helium  reservoir. 
Both  the  housing  and  its  socket  are  made  of  a laminate 
reinforced  with  Kevlar*  cloth.  Four  fiber  glass  rods  with 
reinforcing  spacers  attach  the  foot  of  the  probe  to  its 
top  plate,  which  bolts  to  the  top  of  the  Dewar.  The  top 

* Kevlar  is  a trade  name  for  an  organic  fiber  manufactured 
by  the  Du  Pont  Co. 
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plate  of  the  probe  holds  an  rf  electronics  box  used  to  drive 
the  SQUID,  three  removable  brass  micrometers  used  to  adjust 
slides  during  balancing,  and  a three-axis  fluxgate  magneto- 
meter used  both  to  orient  the  gradiometer  axis  and  to  pro- 
vide vernier  balancing,  as  described  subsequently. 

A.  DESIGN  FEATURES 

Because  susceptibility  of  paramagnetic  substances  is 
inversely  proportional  to  temperature  at  temperatures  of  a 
few  degrees  Kelvin,  we  use,  where  practical,  materials  with- 
in the  helium  bath  that  are  free  of  paramagnetic  impurities 
in  order  to  eliminate  noise  from  magnetization  currents 
driven  by  fluctuations  of  ambient  temperature  and  pressure 
and  of  the  earth's  magnetic  field.  In  addition  to  its  high 
purity,  crystalline  silicon  has  a high  thermal  conductivity 
that  smooths  thermal  gradients  in  the  neighborhood  of  the 
gradiometer  pickup  loops. 

Although  the  laminate  reinforced  with  Kevlar  cloth  is 

effectively  free  of  paramagnetic  impurities,  the  fiber 

glass  laminate  is  contaminated.  Susceptibility  of  the  fiber 

glass  laminate,  resulting  from  paramagnetic  impurities,  is 
* 4 o 

about  2 X 10  at  4.2  K.  Nonetheless,  alternatives  to  a 
fiber  glass  interior  vessel,  such  as  quartz,  are  as  yet  im- 
practical for  the  large  capacity  Dewars  required  for  operating 
at  remote  field  sites. 

Fluctuations  of  the  earth's  magnetic  field  also  drive 
eddy  currents  in  the  aluminized  Mylar  layers  .rapped  around 
the  interior  vessel  to  reduce  radiative  heat  flow.  Because 
dielectric  materials  are  poor  reflectors,  conducting  layers 
are  required  to  reflect  radiant  heat.  Noise  resulting  from 
fluctuating  eddy  currents,  however,  does  not  markedly  in- 
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crease  with  decreasing  frequency  because  amplitude  of  eddy 
current  fluctuations  is  proportional  to  excitation  frequency. 

B.  SUPERCONDUCTING  CIRCUIT 


Figure  2b  is  a diagram  of  the  superconducting  circuit. 

Two  coplanar  pickup  loops,  each  t,45  cm  in  diameter  and 
spaced  25  cm  between  centers,  are  oppositely  wound  in  series, 
so  that  fluctuations  of  the  super  current,  1^^,  in  the  pickup 
circuit  are  driven  by  fluctuations  in  the  net  ambient  flux 
threading  the  pickup  loops.  Fluctuations  in  the  net  flux 
are  coupled  to  the  SQUID  sensor  through  an  rfi  transformer 
containing  a normal  metal  shield  between  its  superconducting 
windings  that  attenuates  interference  at  radio  frequencies 
(above  19  KHz) . The  flux  coupled  to  the  sensor  is  the  pro- 
duct of  the  super  current,  I2/  in  the  coupling  circuit  and 
the  mutual  inductance  between  the  sensor  and  the  field 
coil.  The  SQUID  sensor  itself  is  effectively  an  open  circuit 
for  low  frequency  fluctuations  (Zimmerman,  1971)  . 


Because  the  net  voltage  around  each  superconducting  cir- 
cuit vanishes,  the  flux  coupled  to  the  sensor,  is  pro- 

portional to  the  net  flux  threading  the  pickup  loops,  4>  ; 

ir 

namely. 


. r 1^/2 

- = l<V^p^  ^"2] 

^p  (l+Jl^)  (1+Jl2)-Kt^^l^2 


(la) 


where  J-2  ” ^p^  ^s*  ^1^  and  L2  are 

self  inductances  of  the  circuit  components  labeled  in  Figure 
2b; and  the  coefficients  K and  K are  a measure  of  coupling 

L.  £3 

between  windings  of  the  rfi  transformer  and  between  the  field 
coil  and  the  SQUID  sensor,  respectively.  Equation  la  shows 


that  the  flux  coupled  to  the  sensor  is  greatest  for 

2 -1/2  ^ 
= (1-K  ) , at  fixed  values  of  L and  L , and  that  its 

t p s 

greatest  value  is  given  by 


/ "s' 

1 K .1/2 

\"p) 

1 s P 

' 2 

m 

1+ (1-K^^) 

(lb) 


In  our  instrument,  L = 300  nH,  L = 0.05  nH,  K = 0.7, 

p s s 

and  windings  in  the  rfi  transformer  are  adjusted  to  give  a 
maximum  flux  ratio  of 

= 2.26  X lO'^  (2) 

m 

with  = 0 . 8 . 

From  the  diameter  and  spacing  of  the  piclcup  loops,  we  find 
that  an  ambient  gradient*  of  5.32  pT/m  produces  a net  ambient 
flux  of  one  quantum,  at  the  pickup  loops.  Equation  2 then 

tells  us  that  a gradient  of  1 pT/m  at  most  applies  a flux  of 
4.25  X lo'"*  4>Q/(pT/m)at  the  SQUID  sensor.  A change  in  flux  of 
one  quantum  applied  to  the  SQUID  gives  an  output  of  12.7  volts. 

At  ;iest,  then,  we  expect  a calibration  constant  of  5.4  mV/(pT/m). 
Our  measured  value  is  4.5  mV/(pT/m)  ^ 15%,  which  indicates  that 
the  instrument  is  operating  with  nearly  the  maximum  flux  ratio. 


* We  use  MKS  units  throughout  and  conform  to  the  international 

convention  suggesting  use  of  the  unit  Tesla,  which  is  a 

2 -12  -10  —3 

Weber/m  and  so  IpT/m  = 10  T/m  = 10  G/cm  = 10  gamma/m 

and  Oq  = 2.07  X 10  ^pTm^ . 
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Section  III 

INSTRUMENT  RESPONSE  TO  AMBIENT  MAGNETIC  GRADIENTS 

We  describe  the  response  of  an  axis  of  a superconducting 
magnetic  gradiometer  to  ambient  magnetic  gradients  in  terms 
of  the  differential  magnetic  flux  threading  the  two  pickup 
loops  whose  centers  are  separated  by  a distance  2s.  Areas 
of  the  loops  are  A^  and  and  respective  normals  to  the 

plane  of  each  loop  are  directed  along  unit  vectors  n^^  and  n2 , 
as  shown  in  Figure  3.  Pickup  loops  are  connected  in  opposi- 
tion, so  that  the  instrument  responds  to  the  differential 
magnetic  flux  threading  the  loops;  namely,  O = - <J)2- 


The  relations 


= J d??j_nj^.b(R+s+Pj^»t) 


and 


= J d? 2^2 ' 


express  magnetic  flux  threading  each  loop  in  termj  of  integrals 
of  the  density  of  magnetic  flux,  b (R+s+'^  2^ , t)  and  b (R-s+'t52  » t)  , 
at  points  within  each  loop.  Here,  R is  a position  vector  loca- 
ting the  point  midway  between  loop  centers;  vectors  s and  -s 
locate  loop  centers  with  respect  to  the  midpoint;  and  -p2 
are  radius  vectors  locating  points  within  a looo  contour;  and 
integrals  extend  over  the  area  of  each  loop.  By  expanding 
flux  densities  in  Taylor's  series  about  the  midpoint,  we  express 
the  flux  threading  each  loop  as 


^ ^ 


oo  m 


oo 

1 

(s. 

m = 1 

ml 

(s-V) 

m-k 

(m-k)  I 

k: 

?7  \ ill 

nj^-b(R,t) 


oo 

2 ^2 

ml 

CO  m 

yy 

(-l)"''^(s.V)' 

m=1  k-1 

(m-k) I kl 

n2 • b (R, t) 


kA 


dp2(p2-V)'^n2-b(R,t) 


where  b(R,t)  is  the  magnetic  flux  density  at  the  point  midway 
between  loor  centers.  For  symmetrical  pickup  loops,  integral 
terms  of  odd  order  in  Equations  3c  and  3d  vanish  because  inte- 
grands are  periodic. 

If  area  and  orientation  of  the  loops  are  identical  so  that 
~ ^2  ~ ^ *^1  ~ ~ relation 

$ 1 ^ ^ 

= - y- n.b(R,t) 

2sA  s“  (2m-l); 


+ y_jZ7~~ / dpnf.V)  ^^n;b(R,t) 

(m-k)+l]  i (2k)  : { 

M 

gives  the  response,  and  so  to  first  order 
$/2sA  = (s.V)fi.b(R,t) . 
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Hiaher  order  terms  are  negligible  provided  gradients  of  the  am- 
bient field  are  sensibly  constant  over  the  distance  between  pickup 

2 

loops.  Specifically,  we  require  that  (s/R)  «1,  where  R is  a 
characteristic  scale  of  the  ambient  field.  In  the  field  of  a 
dipole,  for  example,  R is  distance  from  the  dipole. 

Because  of  slight  differences  in  area  and  orientation  of 
the  pickup  loops, however , the  response  includes  an  imbalance 
or  error  term.  Using  loop  number  two  as  a reference  and  de- 
fining A = A2  and  ft  = ^2'  then  express  the  response  of  a 
gradiometer  axis,  Z('s,ft;R,t)  = $/2sA,  to  first  order  as 

Z(9,n;R,t)  ^ r(.s,n;R,t)  + ^-b(lf,t)  (4a) 

where 

r(ft,n;R,t)  = (s-V)n.b(R,t)  (4b) 

and 


To  a first  order  approximation,  then,  the  gradient  term  T(s,n;R,t) 
of  the  response  depends  on  gradients  at  the  point  midway  bet’.veen 
loop  centers,  and  the  imbalance  or  error  term,  "?-b,  is  proportional 
to  the  magnetic  field  at  the  midpoint. 
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by 


We  represent  gradients  of  magnetic  fields  in  free  space 
a matrix  G having  elements  g^j(R,t),  where 


X , 


, A 

MX, 


b) 


(i,j  = 1,2,3) 


(5a) 


in  an  orthogonal  basis  |x^[  . An  element  g^^  represents  the 
gradient  in  a direction  x^  of  the  component  of  magnetic  field 
in  a direction  x..  Consequently,  the  relation 


r(s,fi)  = 


g . . s . n . 
13  1 3 


(5b) 


1/3 


expresses  the  gradient  in  a direction  s of  the  component  of 
magnetic  field  in  a direction  n in  terms  of  matrix  elements 
in  a basis  | } / where  s^  and  n^  are  components  of  the  unit 

vectors  s and  n in  the  same  basis. 

By  choosing  a basis  fixed  at  the  point  midway  between 
loop  centers  with  its  x^  axis  aligned  along  the  gradiometer  axis 
s and  its  x^  axis  aligned  along  n X s,  then,  we  obtain  the  expansion 

r(s,njR,t)  = n2  g23(R/t)  + n^  g23(R/t)  (6) 

for  the  gradient  term  of  the  response.  For  coplanar  pickup 
loops  (transverse  configuration),  n^  vanishes  and  n2  = 1;  for 
coaxial  pickup  loops  (longitudinal  configuration),  n2  vanishes 
and  n^  = 1.  A perfectly  balanced  gradiometer  axis  having  a 
transverse  configuration,  then,  measures  skew  elements  of  a 
gradient  matrix,  and  one  having  a longitudinal  configuration 
measures  diagonal  elements. 

T)ie  axis  of  our  gradiometer  has  a transverse  configuration, 
so  that  its  response  to  ambient  gradients  that  are  sensibly 
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constant  over  the  distance  between  centers  of  the  pickup  loops 
(25  cm)  is  the  sum  of  a transverse  gradient  at  the  point  mid- 
way between  pickup  loops  and  an  imbalance  term  that  is  pro- 
portional to  the  magnetic  field  at  the  midpoint,  as  expressed 
by  the  relation 

E(s,n;R,t)  = g23(H,t)  + 6*b(R,t)  (7) 

The  imbalance  part  of  the  response  resulting  from  small 
differences  in  orientation  and  area  of  the  pickup  loops,  in 
effect,  is  the  response  of  a magnetometer  that  measures  the 
component  of  magnetic  field  in  a direction  6 at  the  midpoint. 

The  gradient  part  of  the  response  is  equivalent  to  a gradient 
of  the  field  of  a magnetic  dipole  located  on  a sphere  of 
unit  radius  about  the  midpoint. 

A.  EQUIVALENT  DIPOLE  OF  A GRADIENT  FIELD 

Because  a luagnetic  field  in  free  space  is  both  nondiver- 
gent  and  irrotational , its  gradients  are  elements  of  a 
symmetric  matrix  that  has  a vanishing  trace.  A matrix  that 
is  both  real  and  symmetric  is  diagonal  in  an  orthogonal 
basis  comprised  of  principal  axes  with  diagonal  elements 
or  eigenvalues  X2»  and  X^*  Since  its  trace  vanishes, 

^1  ^2  ^3  ~ 

We  express  gradie  -s  at  a point  of  a magnetic  field, 
then,  in  terms  of  their  eigenvalues,  X^,  and  eigenvectors 
which  define  directions  of  principal  axes,  by  using 
the  matrix  product 


G = RAR, 


(8) 


where  the  matrix  G represents  gradients  in  a fiducial  basis 
I I . The  matrix  R represents  a rotation  from  the  fiducial 
basis  to  a basis  comprised  of  eigenvectors  e^,  in  which  gra- 
dients are  represented  by  the  diagonal  matrix  / having  ele- 
ments and  A^.  Components  of  eigenvectors  in  the 

fiducial  basis  form  the  columns  of  the  matrix  R;  namely, 
r^j^  = A real  matrix  representing  a rotation  is 

orthogonal,  so  that  RR  = RR  = I,  where  tilde  marks  a trans- 
posed matrix  and  I is  the  unit  matrix. 

By  adding  a magnetic  dipole  to  a gradient  field,  we  find 
that  the  relation 

r = R(A  - (9) 

gives  gradients  at  a field  point.  Here,  is  a diagonal  matrix 
whose  elements  are  eigenvalues  of  the  gradients  of  the  dipole 
field  at  the  pOo.nt,  and  the  matrix  R,  specifies  a rotation  from 

cl 

principal  axes  of  the  ambient  gradient  field  to  principal  axes 
of  the  dipole  gradient  field  ac  the  point.  Equation  9 tells 
us  that  gradients  of  an  ambient  field  are  equal  to  gradients 
of  a dipole  field  at  a point  whenever  R^=I  and  A^=A.  At  each 
point  of  an  ambient  field,  then,  we  use  the  conditions  R =•!  and 
A^=A  to  define  location,  orientation,  and  moment  of  a dipole 
positioned  on  a sphere  of  unit  radius  about  the  point  so  that 
gradients  of  th  ’ dipole  field  at  the  point  are  equal  to  ambient 
gradients . * 

Gradients  at  each  point  of  an  ambient  field,  then,  are 
equivalent  to  the  gradient  field  of  a magnetic  dipole  located 
on  a sphere  of  unit  radius  surrounding  the  point.  Because 
the  sense  of  an  eigenvector  is  indeterminate,  the  dipole  has 
three  images,  so  that  gradients  at  a point  of  an  ambient  field 
are  equivalent  to  the  gradient  field  of  a dipole  located  at 

* Equality  of  eigenvalues  at  a point  is  always  possible  because  - 
eigenvalues  are  ordered.  A,  A^  A_,  and  their  sum  vanishes. 

Phe  Appendix  describes  means  of  determining  location,  orienta- 
tion, and  moment  of  an  equivalent  dipole  from  gradients  of 
the  ambient  field. 
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any  one  of  four  positions  on  the  sphere.  Radius  of  the  sphere 

4 

is  arbitrary  provided  the  ratio  m/r  is  constant,  where  m is 
the  moment  of  the  dipole  and  r,  its  distance  from  the  field 
point.  Location,  orientation,  and  moment  of  equivalent  dipoles 
change  continuously  from  point  to  point. 


B.  GRADIENT  FIELD  OF  A MAGNETIC  DIPOLE 


As  illustrated  in  Figure  4,  we  specify  location  of  a mag- 
netic dipole,  m’,  with  respect  to  the  gradiometer  basis  | | 

in  terms  of  polar  angles  cj)  and  6 defining  direction  of  the 
position  vector  r and  specify  its  orientation  with  respect 
to  the  position  vector  in  terms  of  the  polar  angle  x ^nd 


azimuthial  angle  The  basis 


A 


is  defined  so  that  the 


A A . 

y^  axis  points  along  the  position  vector;  the  ^2  sxis,  along 
r X m;  and  the  y^^  axis,  along  (r  X nt)  X r. 


As  shown  in  the  Appendix,  the  matrix 


Gy(g,x)  = g 


cosx 

0 

sinx 


0 

cosx 

0 


sinx 

0 

-2cosx 


represents  gradients  at  position  r of  a dipole  fiel^  in  the 


basis  |y^ 


, where 


-7 

Q = 2_  and  u = 4tt  x 10  H/m. 

g 4 o 

47Tr 


Equation  10a  tells  us  that  longitudinal  gradients  result  from 
the  component  of  the  dipole  along  r and  that  transverse  gra- 
dients result  from  the  component  perpendicular  to  r.  None- 
theless, each  component  produces  both  longitudinal  and  trans- 
verse gradients  in  the  gradiometer  basis. 


(10a) 


(10b) 
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Three  consecutive  rotations  represented  by  the  matrix 


e,  n)  = + tt/2)  X(Q)  Z (Q) 

bring  the  gradiometer  basis  into  coincidence  with  the  basis 
, where 


ZW  = 


cos  <l> 
sin  ip 
0 


-sin  ij) 
cos  Ip 
0 


and  represents  a rotation  through  an  angle  ijj  about  the 
axis  of  the  gradiometer  basis  and 


X(iJ;)  = 


1 

0 

0 


0 0 
cos  tjj  -sin  Ip 
sin  cos  i|» 


and  represents  a rotation  about  the  x^  axis 


Consequently,  the  matrix 


G (<}),  e,  Q;  g,x)  = fi)G  {g,x)R{4>r 

X y 

represents  gradients  of  the  dipole  field  in  the  gradiometer 
basis.  We  then  find  that 

= g cos  X “ 3/2  {(J),  6)j 

+g  sinx  |cos  A2(<{),  0)  - sin  Q A2((l),  9)j  , 

where  I denotes  the  unit  matrix, 


(11a) 


(lib) 


(11c) 


(12a) 


(12b) 
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(12c) 


0) 


2 2 
2cos  4)Sin  0 

2 

sin24)sin  0 
cos4isin20 


2 

sin2(J)Sin  0 

2 2 
2sin  (j)sin  0 

sin(j)sin20 


cos4)sin20 

sin4)sin20 

2cos^0 


A^((p,  6)  = 


-sin24>sin0 

cos2(f)sin0 

-sin4)cos0 


cos2(|)sin0 

sin2(J)sin0 

COS(J)COS0 


-sin4)cos0 

cos4>cos0 

0 


(12d) 


and 


0) 


cos  4>sin20  (1/2)  sin2(})sin20 

2 

( 1/2 ) sin24>sin20  sin  (J>sin20 
cos4>cos20  sin(j)cos20 


cos(l)cos20 

sin(})cos20 

-sin20 


(12e) 


The  first  term  of  Equation  12b  results  from  the  component  of 
the  dipole  along  r,  and  the  second  term,  from  the  component 
perpendicular  to  1r.  We  see  that  each  component  produces  both 
longitudinal  and  transverse  gradients  in  the  gradiometer  basis. 


C.  GRADIENT  RESPONSE 

From  Equation  12b,  then,  we  find  that  the  relation 


_3 

^23  ~ 2 ^ cosx  sin<|)  sin20 

+ g sinx  (cosJ2cos<|)  cos0-sinJ^  sim})  cos20) 

expresses  the  gradient  term  of  the  instrument  response  in  terms 
of  the  five  parameters  <p,  0,  Q,  X/  and  g that  specify  location, 
orientation,  and  moment  of  a dipole.  Because  gradients  at  a 
point  are  equivalent  to  the  gradient  field  of  a dipole  located 
on  a sphere  of  unit  radius  about  the  point.  Equation  13  com- 
pletely specifies  the  gradient  response. 


(13) 
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The  construct  of  an  equivalent  dipole  furnishes  a convenirnt 
transformation  of  the  five  independent  elements  of  a matrix 
representing  gradients  of  a magnetic  field  in  free  space  to 
five  parameters  {()),  6,  Q,  x»  9)  that  provide  ready  means 

of  describing  the  gradient  response.  For  example,  Equation 
13  tells  us  that  a dipole  located  along  the  x^  axis  of  the 
gradiometer  basis  (())  = 0 cuid  6 = 71/2)  gives  a null  response 
whatever  its  orientation.  Moreover,  oppositely  located  and 
oppositely  directed  dipoles  ((P-*  (p  + v , 9— and  Q-+ 

.give  the  same  gradient  response. 

The  gradient  field  of  a dipole  represents  gradients  of  a 
magnetic  field  at  a point  in  free  space  even  though  the  actual 
source  of  the  field  is  not  a dipole.  Because  of  the  principal 
of  superposition,  separate  equivalent  dipoles  represent  gra- 
dients from  distinct  sources.  The  gradient  field  of  a dipole 
describing  gradients  resulting  from  a sum  of  distinct  sources 
is  simply  the  sum  of  fields  of  equivalent  dipoles  representing 
gradients  from  each  source,  but  the  dipoles  themselves  are  not 
simply  additive  unless  they  are  located  colinearly.  Nonethe- 
less, fields  of  a sum  of  equivalent  dipoles  conveniently  de- 
scribe contributions  of  distinct  sources. 

D.  INSTRUMENT  CALIBRA'  '"ON 


To  calibrate  the  gradiometer,  we  use  a current  oscillating 
in  a coil  placed  equidistant  from  the  center  of  each  pickup 
loop  with  its  dipole  moment  parallel  to  the  line  joining  centers 
of  the  pickup  loops.  We  face  the  plane  containing  pickup 
loops  toward  the  coil  to  obtain  a maximum  response  and  keep 
the  distance  from  the  midpoint  of  the  gradiometer  axis  to  the 
center  of  the  coil  much  larger  than  both  the  spacing  between 
centers  of  the  pickup  loops  (25  cm)  and  the  radius  of  the 
coil  ( 21  cm) . 
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(* 


The  coil,  in  effect,  is  a dipole  located  along  the  axis 
of  the  gradiometer  basis  and  directed  parallel  to  the  axis, 
so  that  ())  = 3ti/2,  0 = tt/2,  = 37r/2 , and  x = tt/2  . The  expression 

3 y m ( t ) 
o 


obtained  from  Equation  13  then  gives  the  gradient  produced  at 
the  gradiometer,  where  m(t)  is  the  dipole  moment  of  the  coil 
and  1 is  its  distance  from  the  gradiometer.  Equation  13  also 
shows  that  misalignment  of  the  coil  and  gradiometer  adds  terms 
of  second  order  in  angular  deviations  and  so  produces  negligi- 
ble error  for  reasonably  precise  alignment.  Errors  result 
largely  from  uncertainties  in  the  distance  of  the  coil  from 
the  gradiometer. 

By  measuring  current  in  the  coil  and  output  voltage  of 
the  gradiometer  at  several  values  of  r,  we  find  that  a gra- 
dient of  1 nT/m  gives  an  output  voltage  of  4.5  V + 15%.  Cali- 
bration is  independent  of  frequency  and  of  the  quantum  state 
of  the  SQUID  sensor. 

We  also  use  the  dipole  used  for  calibration  to  align  the 
triaxial  fluxgate  magnetometer  with  respect  to  the  gradiometer 
basis.  Turning  the  gradiometer  about  its  axis  until  a null 
i response  is  obtained  aligns  the  plane  containing  pickup  loops 

t with  the  plane  defined  by  the  dipole  moment  of  the  calibration 

coil  and  the  gradiometer  axis.  Oscillating  magnetic  field 
components  along  the  x^  and  ^2  axes  of  the  gradiometer  basis 
then  vanish,  so  we  align  corresponding  axes  of  the  fluxgate 
I magnetometer  to  obtain  null  responses  as  well. 
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Section  IV 


GRZ^DIOMETER  BALANCING 

Both  mechanical  and  electronic  means  of  balancing  the 
gradiometer  are  available.  We  balance  mechanically  by  using 
micrometers  on  the  top  plate  of  the  probe  to  move  three  small 
niobium  disks  mounted  on  slides  near  the  top  pickup  loop. 
Because  a superconducting  disk  distorts  the  magnetic  field 
in  its  vicinity,  movement  of  the  disks  adjusts  the  effective 
orientation  and  area  of  the  loop.  Axes  of  the  disks  are 
mutually  orthogonal,  so  that  movement  cf  each  disk,  in  effect, 
adjusts  the  corresponding  component  of  the  imbalance  vector 


Because  precision  of  disk  movements  is  limited,  however, 
we  use  the  triaxial  fluxgate  magnetometer  mounted  on  the  top 
plate  of  the  probe  to  provide  vernier  balancing  electroni- 
cally. We  balance  electronically  by  subtracting  an  adjustable 
fraction  of  the  magnetic  field  components  measured  by  the 
fluxgate  magnetometer  from  the  gradiometer  response. 

To  make  mechanical  and  electronic  balance  adjustments,  we 
rotate  the  gradiometer  axis  in  the  earth's  steady  magnetic 
field.  Because  gradients  of  the  earth's  steady  magnetic 
field  are  small HIO  pT/m) , rotation  in  the  earth's  field  fur- 
nishes a time  varying  field  that  is  nearly  uniform,  so  that  the 
gradiometer  response  effectively  results  from  imbalance  alone. 

A.  BALANCING  IN  A UNIFORM  MAGNETIC  FIELD 

Rotating  a gradiometer  about  an  axis  in  a direction 
in  a uniform  magnetic  field,  ”B,  produces  a response  owing  to 
imbalance  of  its  pickup  loops  that  is  expressed  in  terms  of 
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the  angle  of  rotation  Q by  the  relation 


6-B  = (B 


■fi) 


A 


A a1  -» 

Qx  (Bxfi)j*6^ 


cosQ  + 


A 

(BxQ) 


6 sinO, 
o 


where  is  the  initial  position  of  the  imbalance  vector  6. 
The  imbalance  vector  is  fixed  to  the  gradiometer,  so  that 
6^  depends  on  initial  orientation  of  the  gradiometer.  As  is 
evident  from  Equation  15^ we  define  a fiducial  basis  using 

A 

the  vectors  B and 


As  before,  we  choose  a right-handed  basis  fixed  to  the 
point  midway  between  centers  of  the  pickup  loops  with  its 
5^2  axis  along  their  line  of  centers  and  its  ^2  axis  normal 
to  the  plane  containing  pickup  loops.  We  denote  components 
of  the  imbalance  vector  in  the  gradiometer  basis  by  6^^,  62, 
and  6 ^ « 

In  our  balancing  procedure,  the  axis  of  rotation  is 
vertical,  and  the  earth's  magnetic  field  provides  an  effec- 
tively uniform  field.  We  use  two  initial  orientations  of 
the  gradiometer  basis  in  order  to  make  sequential  and  inde- 
pendent adjustments  of  components  of  the  imbalance  vector. 

First,  we  initially  orient  the  gradiometer  with  the 
axis  directed  vertically  upwards;  the  ^2  ^xis,  northward; 
and  the  x^  axis,  eastward.  Rotating  the  gradiometer  about 
its  $2  axis  through  an  angle  then  gives  the  response 

6*B  = -62Bsin3  + (62sin  Q+  62Cosf2 ) BcosB  , 

where  6 is  the  local  dip  angle  of  the  earth's  magnetic  field 
and  B is  its  magnitude.  Equation  16  tells  us  that  a cosine 
response  corresponds  to  the  62  component  of  the  imbalance 
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(15) 


(16) 


vector  and  that  a sine  response  corresponds  to  the  6^^  com- 
ponent. The  6 2 component  determines  a reference  level  of 
the  sinusoidal  response.  To  nullify  components  6^^  and  62 
of  the  imbalance  vector,  we  first  adjust  the  niobium  disk 
having  its  axis  parallel  to  until  rotation  produces  a 
cosine  response  and  then  adjust  the  disk  having  its  axis 
parallel  to  “$2  'J^til  rotation  produces  a flat  response. 

Next,  we  initially  align  the  gradiometer  with  the  “x^ 
axis  parallel  to  the  earth's  magnetic  field  and  the  ^2 
axis  directed  westward.  Rotating  the  gradiometer  about  a 
vertical  axis  through  an  angle  Q then  gives  the  response 

-6.1b  = (5,cos8  + 6^sin6)  Bsing  + |6»sinf2 

1 J I.Z 

-(6j^sin3  - 52Cos6)cosJ^  BcosB 


Equation  17  tells  us  that  a cosine  response  corresponds  to 
components  6^  and  <53  of  the  imbalance  vector  and  that  a sine 
response  corresponds  to  the  62  component.  To  complete  balanc- 
ing the  gradiometer,  then,  we  only  adjust  the  disk  having  its 
axis  parallel  to  until  rotation  produces  a flat  response, 
because  the  first  procedure  nullifies  components  6^^  and  62- 

B.  LIMITS  ON  GRADIOMETER  BALANCE 


Both  limited  precision  of  disk  movements  and  mutual  in- 
ductances between  disks  limit  the  balance  attainable  mechani- 
cally, and  so  prevent  attainment  of  precisely  flat  responses 

during  rotations.  Mechanical  adjustments  reduce  components 

-5  -1* 

of  the  imbalance  vector  to  values  of  the  order  of  10  m 
To  reduce  imbalance  further,  we  electronically  subtract 
a fraction  of  magnetic  field  components  measured  by  the  flux- 


* Balance  of  a gradiometer  is  also  commonly  expressed  in  terms 
of  a common  mode  rejection  ratio  (CMRR)  that  is  the  magnitude 
of  the  imbalance  vector  multiplied  by  the  distance  between 
centers  of  its  pickup  loops.  For  our  gradiometer,  an  im- 
balance of  10"5  m“^  gives  a CMRR  of  2.5  ppm. 
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gate  magnetometer  from  the  gradiometer  response  during  rota- 
tions, after  reaching  the  limit  of  mechanical  adjustments. 
Combined  mechanical  and  electronic  adjustments  reduce  com- 
ponents of  the  imbalance  vector  to  values  of  about  4 X 10  ^ 


Further  refinement  of  gradiometer  balance  presently  is 
limited  by  hysteresis  observed  during  rotation  when  the  im- 
balance is  somewhat  less  than  10  ^ m The  record  in  Figure 
5 illustrates  the  effect.  The  trace  records  gradiometer  re- 
sponse during  a clockwise  rotation  followed  by  a counter  clock- 
wise rotation  with  the  gradiometer  axis  vertical.  Rotations 
are  made  stepwise  in  sixteen  increments  of  22.5°,  as  shown  in 
Figure  6#  starting  with  the  "^2  axis  pointing  northward.  Spikes 
in  Figure  5 mark  movements  between  incremental  positions  and 
result  from  eddy  currents  excited  in  aluminized  Mylar  layers 
used  to  insulate  the  helium  reservoir.  Plateaus  mark  values 
recorded  while  dwelling  at  incremental  positions.  As  is 
evident,  values  recorded  at  incremental  positions  are  not 
unique  and  depend  on  the  sense  of  rotation;  namely;  the  value 
recorded  at  a position  during  a clockwise  rotation  differs 
from  the  value  recorded  at  the  same  position  during  a counter 
clockwise  rotation. 

Figure  7 shows  the  hysteretic  signature  of  the  effect. 

The  trace  records  response  during  a stepwise  oscillation  of 
the  '^2  axis  about  position  5 marked  in  Figure  6.  Oscillation 
begins  with  a clockwise  movement  of  two  increments  from  position 
5 to  position  7,  then  reverses  with  a counter  clockwise  movement 
to  position  5,  continues  with  a counter  clockwise  movement  of 
two  increments  to  position  3,  returns  to  position  5 with  a clock- 
wise movement,  and  repeats  the  cycle  thereafter.  As  is  evident 
in  Figure  7,  the  value  recorded  at  oosition  5 depends  on  the 
direction  of  approach;  namely,  the  value  recorded  during  a 
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clockwise  crossing  differs  from  the  value  recorded  during  a 
counter  clockwise  crossing.  The  difference  depends  on  ampli- 
tude of  the  oscillation. 

We  observe  a similar  hysteretic  signature  when  oscillations 
are  made  with  a crescent  wrench  taped  to  the  top  of  the  Dewar 
and  so  suspect  that  the  hysteretic  effect  results  from  ferro- 
magnetic contaminants  in  the  vicinity  of  the  pickup  loops. 

Although  hysteresis  precludes  further  reduction  of  gradiometer 
imbalance,  the  residual  imbalance  is  insignificant  for  our 
immediate  purposes. 

Eddy  currents  excited  in  the  aluminized  Mylar  insulation 
during  rotations  do  not  limit  refinement  of  gradiometer  balance 
but  do  require  that  rotations  be  made  stepwise.  Continuous 
rotation,  in  effect,  averages  values  at  spikes  and  plateaus 
and  so  is  misleading  for  purposes  of  refining  balance.  The 
response  resulting  from  eddy  currents, is  consistent  with 
the  description  given  by  the  relation 

I = C = CBcosB  aT— ^ sinfJ  , (18) 

® dt  \dt/ 

where’s  is  the  earth's  magnetic  field; "A,  the  effective  area 
of  eddy  current  loops;  C,  a coupling  constant;  and  is  the 
angle  of  rotation,  reckoned  positive  for  a clockwise  rotation. 

The  vector  representing  effective  area  of  eddy  current  loops 
is  normal  to  the  plane  containing  pickup  loops.  Allowing 
for  a difference  in  rates  of  movement  between  increment  posi- 
tions, Figure  5 shows  that  spikes  resulting  from  eddy  currents 
^j-0  greatest  at  tt/2  and  3tt/2,  least  at  0 and  it,  and  of  opposite 
seise  for  clockwise  and  counter  clockwise  rotations,  as  described 
by  Equation  18. 
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Gradients  of  the  earth's  magnetic  field  are  a more 
fundamental  limit  on  refinement  of  gradiometer  balance.  Gra- 
dients of  the  earth's  magnetic  field  are  of  the  order  of  10  pT/m 
and  so  give  a response  comparable  to  the  response  owing  to  im- 

_7 

balance  when  magnitude  of  the  imbalance  vector  is  about  5 X 10 
m ^ or  less.  In  the  absence  of  the  limit  imposed  by  hysteresis, 
reduction  of  gradiometer  imbalance  below  about  10  m is 
limited  by  gradients  of  the  earth's  magnetic  field. 

C.  BALANCING  IN  A NONUNIFORM  MAGNETIC  FIELD 

Because  a matrix  representing  gradients  of  a magnetic  field 
in  free  space  is  symmetric,  it  is  diagonal  in  a basis  comprised 
of  principal  axes.  Transverse  gradients  vanish  along  principal 
axes,  so  that  a perfectly  balanced  gradiometer  having  a trans- 
verse configuration  gives  a null  response  when  its  axis  is 
aligned  along  a principal  axis  of  an  ambient  gradient  field. 

The  Appendix  describes  both  a means  of  determining  directions 
of  principal  axes  of  a gradient  field  from  values  of  its  matrix 
elements  measured  in  a fiducial  basis  and  response  of  a per- 
fectly balanced  gradiometer  during  rotation  about  an  arbitrary 
axis . 

To  refine  balance  of  our  gradiometer  in  the  presence  of  a 
nonuniform  magnetic  field,  then,  we  align  its  axis  along  a 
principal  axis  of  the  ambient  gradient  field  and  rotate  it 
about  its  axis.  Because  transverse  gradients  vanish  along 
a principal  axis,  the  response  during  rotation  results  from 
imbalance  alone.  We  then  adjust  components  of  the  imbalance 
vector  to  obtain  a flat  response  during  rotation.  Precise 
determination  of  directions  of  principal  axes,  however,  re- 
quires a finely  balanced  gradiometer,  so  we  devisa  an  iterative 
procedure. 
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We  suppose  that  directions  of  principal  axes  of  the  ambient 
gradient  field  are  known  approximately.  For  example,  two  prin- 
cipal axes  of  the  earth's  dipolar  magnetic  field  lie  in  a ver- 
tical plane  containing  the  magnetic  field  vector,  and  the  third 

* 

principal  axis  is  normal  to  the  plane.  If  ambient  gradients 
exceed  about  10  nT/m,  then  we  use  the  gradiometer  itself  to 
estimate  approximate  directions  of  principal  axes,  because 
magnitude  of  its  imbalance  vector  without  refinement  is  at 
most  about  lO”^  m”^.  Rotating  the  gradiometer  about  its  axis 
when  it  is  aligned  with  the  local  magnetic  field  vector  gives 
an  estimate  of  the  magnitude  of  ambient  gradients,  since  the 
response  owing  to  imbalance  then  vanishes. 

The  relation 

with 

tanilo  = V13/Y23 

gives  the  response  owing  to  ambient  gradients  when  the  gradio- 
j^0ter  is  rotated  through  an  angle  about  its  axis . Coeffi- 
cients Yjl3  ^23  elements  of  the  matrix  representing 

the  ambient  gradient  field  in  the  basis  defined  by  initial 
orientation  of  the  gradiometer.  If  the  gradiometer  axis  is 
aligned  along  a principal  axis  of  the  ambient  gradient  field, 
then  Yi3  and  Y23  venish,  so  that  the  response  owing  to  ambient 
gradients  during  rotation  is  flat.  If  polar  angles  defining 
direction  of  the  axis  of  rotation  in  a fiducial  basis  deviate 
by  amounts  A(|)  and  A0,  respectively,  from  polar  angles  defining 

*The  Appendix  describes  the  gradient  field  of  a magnetic  dipole 
in  terms  of  its  eigenvalues  and  directions  of  its  principal 
axes . 


(19a) 


(19b) 
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direction  of  the  principal  axis  corresponding  to  eigenvalue 
say,  then 

Y^3  ^1/2  1^3X3  - (2X^  + X3)cos2TpJ  sin0p(Act)) 

-1/2  (2X^  + X3)  sin24'p  (A0) 


and 


Y2s'^  -1/2(2Xi  + X3)  sin24'pSin0p(A(j)) 

+ 1/2  [^3X3  + (2Xj^  + X3)cos24'p]  (A0) 

where  $ , 0 and  are  Euler  angles  (defined  in  Figure  Al) 
specifying  a rotation  from  the  fiducial  basis  to  the  basis 
comprised  of  principal  axes  and  X^,  X2,  and  X3  are  corresponding 
eigenvalues . 

We  choose  a fiducial  basis jxj  with  its  X,  axis  directed 
vertically  downwards,  its  X^^  axis,  northward,  and  its  X2  axis, 
eastward.  To  align  the  gradiometer  axis  with  a principal  axis, 
we  start  with  the  gradiometer  basis  coincident  with  the  fiducial 
basis  and  first  rotate  the  gradiometer  through  an  angle  0 
about  the  axis  and  then  through  an  angle  (j)  about  the  X3 
axis,  where  the  polar  angles  (p  and  9 approximately  define 
direction  of  the  principal  axis  in  the  fiducial  basis.  We  then 
rotate  the  gradiometer  about  its  axis  and  mark  the  amplitude 
of  its  response.  We  next  make  an  incremental  rotation  5(j) 
about  the  “X3  axis  and  again  rotate  the  gradiometer  about  its 
axis  and  mark  the  amplitude  of  its  response.  We  continue 
making  incremental  rotations  about  th?.  ^'3  axis  followed  by 
rotations  about  the  gradiometer  axis  until  we  determine  the 
polar  angle  + <54*  which  amplitude  of  the  response  is  least. 


(19c) 


(19d) 
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With  the  gradiometer  axis  fixed  at  the  polar  angle  (|)  + 6(f), 
we  carry  out  a similar  procedure  to  determine  the  polar  angle 
0 + 66  at  which  amplitude  of  the  response  is  least  and  so 
locate  the  polar  angles  cf)  + 6c()  and  0 + 60  defining  direction 
of  the  principal  axis  more  precisely. 

The  response  resulting  from  imbalance  is  sensibly  constant 
for  incremental  changes  6c{)  and  ci0  in  the  polar  angles,  but 
the  response  resulting  from  ambient  gradients  changes  markedly 
for  small  deviations  of  the  axis  of  rotation  from  the  princi- 
pal axis.  Nonetheless,  imbalance  eventually  masks  changes 
in  amplitude  of  the  response  to  rotation  as  the  axis  of  rota- 
tion approaches  the  principal  axis.  When  incremental  changes 
in  polar  angles  defining  direction  of  the  axis  of  rotation 
produce  imperceptible  changes  in  amplitude  of  the  response, 
we  then  adjust  components  6^^  and  62  of  the  imbalance  vector 
to  reduce  the  response  owing  to  imbalance,  following  a pro- 
cedure similar  to  that  used  in  a uniform  field. 

After  refining  'he  balance,  we  refine  alignment  of  the 
gradiometer  axis  with  the  principal  axis  by  again  making  in- 
cremental changes  in  polar  angles  until  we  determine  the 
direction  of  the  axis  of  rotation  for  which  amplitude  of  the 
response  is  least.  We  then  again  refine  the  balance  and  re- 
peat the  procedure  until  further  refinements  are  imperceptible. 

The  iterative  procedure  affords  means  of  reducing  components 
6j^  and  62  of  the  imbalance  vector.  To  reduce  the  6^  component, 
we  rotate  the  gradiometer  about  its  axis  through  an  angle 
X after  aligning  its  axis  with  the  principal  axis.  The  relation 

g23(x)  = 1/2  [3X3  + (2X^  + X3)cos2'l'p]  sin2x 


(20) 
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gives  the  response  resulting  from  ambient  gradients.  The 
response  resulting  from  imbalance,  however,  is  sinusoidal 
with  respect  to  x rather  than  Consequently,  we  adjust 

the  6 2 component  of  the  imbalance  vector  during  rotation 
until  the  response  is  proportional  to  sin  2x- 

Because  our  present  means  of  rotating  the  gradiometer  only 
provides  for  rotation  about  a horizontal  and  a vertical  axis, 
the  suggested  procedure  for  balancing  in  a nonuniform  ambient 
magnetic  field  is  as  yet  untried. 
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Section  V 


INSTRUMENT  PERFORMANCE 


Past  reports  (Clarke,  1974;  Wynn  et  al. , 1975;  Zimmerman 
and  Frederick,  1971)  have  specified  performance  of  super- 
conducting magnetic  gradiometers  in  terms  of  a noise  level 
without  regard  to  spectral  shape.  Recent  measurements  (Clarke 
et  al . , 1975)  of  noise  spectra  of  a superconducting  quantum 
interference  device  (SQUID) , however,  show  that  their  noise 
power  density  increases  approximately  in  inverse  proportion 
to  frequency  as  frequency  falls  below  about  1 Hz.  Moreover, 
fluctuations  of  ambient  temperature  and  pressure  and  of  the 
earth's  magnetic  field  drive  fluctuating  magnetization  currents 
and  eddy  currents  in  materials  within  the  Dewar  that  generate 
noise  whose  power  density  characteristically  increases  with 
decreasing  frequency  as  well.  Characterizing  instrument 
performance  in  terms  of  a noise  level,  then,  is  inadequate 
in  the  frequency  range  of  interest  (f^l  Hz) . 

In  what  follows,  we  first  present  spectral  data  in  the 

-4 

frequency  range  5 X 10  Hz  to  20  Hz  that  characterize  per- 
formance of  the  instrument,  which  is  designed  insofar  as  prac- 
tical to  eliminate  sources  of  noise  driven  by  fluctuations  of 
ambient  temperature  and  pressure  and  of  the  earth's  magnetic 
field,  and  demonstrate  that  its  performance  in  a magnetically 
quiet  environment  is  limited  by  noise  of  its  SQUID  sensor  at 
frequencies  below  about  0.1  Hz.  We  then  compare  performance 
demonstrated  using  a torodial  point-contact  type  SQUID  sensor 
with  performance  expected  using  a chin-film  dc  type  SQUID,  in 
order  to  mark  limits  on  enhancing  instrument  performance. 
Finally,  we  delineate  improvements  in  instrument  design. 
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A.  GRADIOMETER  NOISE  SPECTRA  IN  A MAGNETICALLY  QUIET  ENVIRONMENT 


To  characterize  performance  of  the  gradiometer  in  the  fre- 

_4  * 

quency  range  5 X 10  Hz  to  20  Hz,  we  made  a number  of  records 
of  its  output,  ranging  from  12  to  24  hours  in  duration,  for 
a variety  of  ambient  conditions  during  an  observational  period 
of  several  weeks  at  the  La  Posta  Astrogeophysical  Observatory, 
which  is  located  in  the  mountains  about  70  miles  east  of  San 
Diego,  California,  at  an  altitude  of  1188  m MSL.  Coordinates 
of  the  observatory  are  116°  25'  6"  west  longitude  and  32°  40' 

39"  north  geodetic  latitude  (Bleiweiss  and  Wefer,  1975) . We 
present  spectral  data  based  on  a five  hour  long  segment  of  a 
record  that  is  representative  of  the  lowest  instrument  noise 
consistently  observed. 

Figure  8 shows  the  spectra]  power  density  in  the  frequency 
range  5 X lo”"^  Hz  to  20  Hz  characterizing  performance  of  the 
instrument  during  a quiet  period  observed  at  night.  We  esti- 
mate the  spectrum  below  0.1  Hz  using  a standard  fast  Fourier 
Transform  algorithm  and  average  ten  successive  estimates  made 
from  1/2  hour  segments  of  the  bime  series.  We  obtain  spectral 
data  at  frequencies  above  0.1  hV^^  using  a real-time  spectrum 
analyzer  having  an  averaging  cap^ility. 

For  frequencies  below  0.1  Hz,  w^  assume  a spectrum  of  the 
form  S(f)  = S(f^)(f  /f)^  and  then  usle  a least  squares  criterion 
to  fit  averaged  spectral  data  to  a straight  line  in  logarithmic 
coordinates  to  obtain  values  of  y and  S(f^).  Values  corresponding 
to  the  spectrum  shown  in  Figure  8 are  y = 1.3  and  S(f^)  = 0.002 
(pT/m)^/Hz  at  f^  = 1 Hz,  for  f<0.1  Hz.  Spectral  power  densities 

* We  use  a low-pass  filter  having  its  half-power  point  at  0.3 
Hz  to  preclude  aliasing  and  digitally  record  on  magnetic 
tape  once  a second  to  obtain  data  at  frequencies  below  0.1  Hz. 

We  obtain  spectral  data  at  frequencies  above  0.1  Hz  using  a 
real-time  spectrum  analyzer  having  an  averaging  capability. 


30 


in  the  frequency  range  0.1  Hz  to  20  Hz  correspond  to  a white 

2 

noise  level  of  0.03  (pT/m)  /Hz.  Additional  estimates  made 
2 3 

near  10  Hz  and  10  Hz  correspond  to  the  same  white  noise 
level . 

We  present  spectra  characterizing  instrument  performance 
during  quiet  times  in  order  to  demonstrate  that  the  inherent 
noise  spectrum  of  the  gradiometer  is  in  practice  limited  by 

the  intrinsic  noise  spectrum  of  the  SQUID  sensor  alone  at 

★ 

frequencies  below  0.1  Hz.  Data  points  in  Figure  9 depict 
the  spectrum  of  intrinsic  noise  of  a torodial,  point-contact 
type  SQUID  having  the  same  design  as  that  used  in  the  gradio- 
meter. The  solid  line  delineates  the  spectrum  of  gradiometer 
noise  shown  in  Figure  8.  Figure  9 shows  that  the  spectrum 
of  gradiometer  noise  closely  approaches  that  of  the  SQUID  at 
frequencies  below  0.1  Hz.  Parameters  characterizing  the  form 

of  the  SQUID  spectrum  have  the  values  y = 1.1  and  S(f  ) = 2.4 

-3  2 ° 

X 10  (pT/m)  /Hz  at  f^  = 1 Hz  as  compared  to  the  values  y 

=1.3  and  S(f  ) = 2.1  X lO"^ (pT/m) ^/Hz  at  f„  = 1 Hz  of  the 
o o 

gradiometer  s[  ’ctrum. 

At  frequencies  above  0.1  Hz,  however,  the  gradiometer  noise 
level  is  appreciably  higher  than  the  SQUID  noise  level.  The 
increase  in  gradiometer  white  noise  over  that  of  the  SQUID 
largely  results  from  Johnson  noise  coming  from  the  normal 
metal  shield  in  the  transformer.  We  expect  the  shield  to 
increase  the  noise  level  by  a factor  of  two  or  so,  which  is 
comparable  to  the  observed  increase  of  nearly  a factor  of 
four.  Without  an  rfi  transformer,  we  would  expect  to  attain 
a white  noise  level  limited  by  intrinsic  white  noise  of  the 

* Data  presented  in  Figures  9 and  10  are  provided  by  courtesy 
of  John  Clarke  and  his  research  group  at  the  University  of 
California,  Berkeley,  California. 


31 


SQUID  sensor.  Presuming  a maximum  flux  ratio  of  (<p  /(^  ) 

-3  s P n 

= 4.52  X 10  without  an  rfi  transformer  (obtained  from 

Equation  lb  with  K = 1) , then,  we  would  expect  a white  noise 

^ - 3 2 

level  for  the  gradiometer  of  about  10  (pT/m)  /Hz,  which  is 

-2 

comparable  to  the  noise  level  of  3 X 10  (pT/mV'VITr  reported 
Wynn  et  al . , (1975)  that  largely  results  from  intrinsic 
noise  of  their  sensor. 

Although  intrinsic  noise  of  the  SQUID  sensor  dominates 
the  noise  spectrum  of  the  gradiometer  below  0.1  Hz  during 
quiet  observational  periods,  we  observe  occasional  increases 
in  slope  and  level  of  spectra  below  0.1  Hz  that  apparently 
are  associated  with  semidiurnal  effects.  We  expect  to  charec- 
terize  sources  of  noise  driven  by  fluctuating  ambient  condi- 
tions by  spectrally  analyzing  fluctuations  of  ambient  tempera- 
ture and  pressure  and  of  the  earth's  magnetic  field  recorded 
at  the  field  site. 

B.  INSTRUMENT  IMPROVEMENT 

Because  intrinsic  noise  of  its  SQUID  sensor  limits  instru- 
ment performance,  replacing  the  torodial  point-contact  type 
SQUID  in  the  instrument  with  a thin-film  dc  type  SQUID  (Clarke 
al . , 1975)  is  an  immediate  means  of  enhancing  instrument 
performance.  Measurements  of  noise  spectra  of  the  two  types 

of  SQUID  sensors  show  that  spectral  density  of  the  thin-film 

2 

dc  type,  in  units  of  (f^/Hz,  is  appreciably  smaller  than  spec- 
tral density  of  the  torodial  point-contact  type  at  frequencies 
below  about  0.1  Hz  (Clarke , personal  communication,  1975). 

Moreover,  self-inductance  of  a thin-film  dc  type  SQUID  is  about 
1 nH  compared  to  0.05  nH  for  a torodial  point-contact  type. 

Equation  lb  then  tells  us  that  coupling  of  flux  to  a thin-film 
dc  type  sensor  is  greater  than  coupling  to  a torodial  point-contact 
type  sensor,  in  the  same  superconducting  circuit.  Namely,  equi- 

valent  noise  energy  in  the  pickup  loops.  /L  , is  proportional 

P P 
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. 2 
to  noise  energy  of  the  SQUID  sensor,  4’g/Lg/  for  fixed  values 

of  the  coupling  coefficients  K and  K . The  ratio  (p^/L  pro- 

S L.  s s 

vides  a f igure-of-iner it  for  comparing  performance  of  SQUID 
sensors  (Classen , 1975) . 

Data  points  in  Figure  10  depict  the  spectrum  of  intrinsic 

2 

noise  of  a thin-film  dc  type  SQUID  in  units  of  marked 

on  the  right  hand  coordinate  scale.  The  light  dashed  curve 

delineates  the  spectrum  of  intrinsic  noise  of  a torodial 

2 

point-contact  type  sensor,  in  units  of  4>^/Hz,  as  depicted  by 

data  points  in  Figure  9.  The  solid  line  marks  the  spectrum 

2 

of  inherent  noise  of  the  gradiometer  in  units  of  (pT/m)  /Hz, 

marked  on  the  left  hand  coordinate  scale,  as  depicted  by  data 

points  in  Figure  8.  The  heavy  dashed  curve  marks  the  spectrum 

of  inherent  noise  expected  of  a gradiometer  operating  with  a 

thin-film  dc  type  SQUID.  We  obtain  the  expected  gradiometer 

2 

noise  spectrum  in  units  of  (pT/m)  /Hz  from  the  spectrum  of 

intrinsic  noise  of  a thin-film  dc  type  sensor  in  units  of 

2 -3 

$^/Hz  by  using  the  conversion  factor  1.56  X 10  <J>^/(pT/m), 

which  accounts  for  the  change  in  self-inductance  of  the  SQUID 

in  the  superconducting  circuit. 

-3 

^t  frequencies  belcw  10  Hz,  then,  intrinsic  noise  of 
the  t.. in-film  dc  type  SQUID  limits  expected  performance  of 
the  gradiometer,  and  at  frequencies  above  10  ^ Hz,  Johnson 
noise  from  the  normal  metal  shield  in  the  rfi  transformer 
limits  performance.  Figure  10  shows  that  the  expected  en- 
hancement of  instrument  performance  is  substantial  at  fre- 
quencies ' alow  about  0.1  Hz.  For  example,  in  the  frequency 
band  from  0.001  to  0.01  Hz,  RMS  instrument  noise  is  reduced 
from  0.17  pT/m  to  0.017  pT/m,  a tenfold  enhancement  of  per- 
formance . 
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one 


We  suggest  two  improvements  in  instrument  design  — 
trivial  and  the  other  substantial.  First,  scoring  aluminized 
Mylar  layers  wrapped  around  the  interior  vessel  of  the  Dewar 
reduces  strength  of  eddy  currents  excited  in  the  superinsula- 
tion. Second,  a triaxial  superconducting  magnetometer  should 
be  incorporated  at  the  midpoint  of  the  gradiometer  axis  both 
to  provide  nore  accurate  means  of  balancing  and  aligning  the 
gradiometer,  ami  to  afford  measurement  of  the  magnetic  field 
as  well  as  its  gradients  at  a point. 
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Section  VI 


TECHNIQUES  FOR  SUPPRESSING  NOISE  FROM  NEARBY  MAGNETIC  OBJECTS 

Fluctuating  gradients  of  magnetic  fields  resulting  from 
ionospheric  currents  are  too  small  to  be  directly  observed 
at  the  ground.  Nonetheless,  ionospheric  currents  excite 
magnetization  and  eddy  currents  in  magnetic  objects  near  a 
gradiometer  and  so  indirectly  produce  local,  fluctuating 
gradients.  Moreover,  objects  magnetized  by  the  earth's 
magnetic  field  produce  sharp,  steady  gradients  in  their 
vicinity.  Slight,  irregular  movements  of  a gradiometer  near 
magnetic  objects,  then,  produce  gradient  fluctuations  as 
well . * 

Fluctuating  gradients  owing  to  nearby  magnetic  objects 
can  mask  ambient  gradient  fluctuations  of  interest.  Con- 
sequently, we  consider  techniques  for  suppressing  gradient 
fluctuations  resulting  from  small,  irregular  movements  of 
a gradiometer  relative  to  a steady  gradient  field  and  from 
fluctuating  magnetization  currents. 

As  we  have  shown,  gradients  of  a magnetic  field  at  a 
point  in  free  space  are  equivalent  to  the  gradient  field  of 
a magnetic  dipole  located  on  a sphere  of  unit  radius  about 
the  point.  Changes  in  location,  orientation,  and  moment  of 
an  equivalent  dipole  then  describe  gradient  fluctuations  owing 
to  magnetic  objects.  Two  techniques  for  suppressing  fluctua- 
tions are  available;  (1)  aligning  the  gradiometer  to  obtain 
a null  response  and  (2)  using  a magnetic  dipole  to  cancel 
gradients  of  an  equivalent  dipole. 

* Because  of  residual  gradiometer  imbalance,  irregular  move- 
ments in  a uniform  magnetic  field  also  produce  gradient 
fluctuations,  but  the  fluctuations  are  comparitively  negli- 
gible for  ambient  gradients  exceeding  10  pT/m. 
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A.  SUPPRESSION  OF  NOISE  FROM  IRREGULAR  MOVEMENTS 


Small  displacements  of  a gradiometer  in  a steady  gradient 
field  produce  small  changes  in  location,  orientation,  and 
moment  of  the  equivalent  dipole  describing  gradients  of  the 
ambient  field  in  the  neighborhood  of  the  gradiometer.  Changes 
in  location  result  both  from  virtual  rotations  of  the  dipole 
accompanying  small  displacements  at  fixed  gradiometer  orien- 
tation and  from  actual  gradiometer  rotations,  but  changes  in 
orientation  and  moment  result  from  displacement  alone. 

From  Equation  12a  that  gives  gradients  in  a gradiometer 
basis  in  terms  of  the  gradient  field  of  a dipole,  we  find  that 
the  relation 


G + 6G  = R R^  (G  + 6G 
xo  X o 6 yo  y 6 o 


(23a) 


expresses  gradients  in  the  gradiometer  basis  following  a small 
displacement  from  an  initial  point  marked  by  subscript  zero, 
where 


6G. 


= 


vo 


(g  . 

y'^o 


+ ^/2) 


(23b) 


The  matrix  R^  represents  an  infinitesimal  rotatio  that  we 
express  in  terms  of  angular  deviations  in  location  and  orien- 
tation of  the  dipole  by  the  sum 


R^  = I + 60cosn^e^  - (<5<|)Sin0^cosJ^^  + 60sinfJj^)e^ 


+ + 64)  (cos0Q-sin0^sinfi^)  j 


(23c) 
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From  Equation 
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we  then  find 
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gives  first  order  changes  of  gradients  in  the  gradiometer  basis 
following  a small  displacement  in  a steady  gradient  field,  which 
is  represented  by  the  matrix  G^^  at  the  initial  point.  Equation 
24  shows  that  changes  in  gradients  are  proportional  to  the 
initial  values  represented  by  the  matrix  G^^. 


To  suppress  gradient  fluctuations  resulting  from  irregular 
motions,  then,  we  first  judiciously  choose  initial  orientation 
of  the  gradiometer  basis  with  respect  to  the  equivalent  dipole. 
We  align  the  gradiometer  so  that  the  equivalent  dipole  lies 
along  the  axis  with  its  axis  in  the  plane  containing  pick- 
up loops,  which  is  normal  to  the  axis  of  the  gradiometer 
basis.  Namely,  we  take  = 0,  0^  = tt/2,  and  = tt/2,  so 
that 


G 


XO 


-2cosXq 

0 


0 sinx^ 

cosx  0 

o 

0 cosXq 


(23d) 


(24) 


(25) 
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Slight  deviations  in  initial  alignment  give  fluctuations  of 
second  order  during  irregular  motions  and  so  are  negligible. 


Gradient  changes  resulting  from  small  displacements  or 
rotations  from  the  initial  orientation  are  then  given  by  the 
expression 


6G  = + 6yG  + 6fi(e  G - G e„) 

''  g xo  xo  X xo  xo  X 


X 


+ 60  (e  G - G e ) + 6(He  G - G e ) 
y xo  xo  y z xo  xo  z 


(26a) 


where 


G = g 
xo  ^o 


2sinx^ 

0 

cosx^ 


0 

-sinx^ 
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cosXj 
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•sinx^  J 


(26b) 
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(26c) 
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(26d) 


and 


G G ““  G G 
Z XO  xo  z 


0 

-3cosx, 


-3cosx^ 

0 

sinx„ 


sinx. 


(26e) 


As  a result,  we  see  that  changes  in  the  transverse  gradient, 
^22'  the  gradiometer  result  from  effective  rota- 

tions about  the  gradiometer  axis  alone;  namely,* 


* The  same  result  is  obtained  by  dif  ferentiating  Equation  13 
that  gives  the  gradient  response  in  terms  of  parameters 
specifying  location,  orientation  i and  moment  of  an  equivalent 

dipole . 
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(27) 


6g23  = sinx^  . 

Magnitude  of  a gradient  fluctuation  is  then  the  product  of 
an  effective  angular  deviation  and  a transverse  gradient 
resulting  from  the  component  of  the  equivalent  dipole  along 
the  gradiometer  axis  (the  ^3  axis) . The  effective  angular 
deviation  is  a sum  of  an  angular  deviation  6(})^owing  to 
virtual  rotation  of  the  equivalent  dipole  accompanying  dis- 
placement and  an  actual  angular  deviation  owing  to  rota- 

tion of  the  gradiometer  basis. 

To  suppress  gradient  fluctuations  resulting  from  effective 
rotations  about  the  gradiometer  axis,  we  position  a dipole 
along  the  axis  of  the  gradiometer  basis  with  its  axis 
parallel  to  the  x^  axis  and  adjust  its  moment  so  that  it  can- 
cels the  component  of  the  equivalent  dipole  along  the  ^3  axis. 
Specifically,  we  position  the  dipole  so  that  = 0,  6^  = it/2, 

= :’tt/2,  and  x*^  = tt/2  and  adjust  its  moment  so  that  g^  = g 
o o 

sinx^.  Changes  in  the  transverse  gradient  measured  by  the 
gradiometer  that  result  from  small  displacements  in  the  steady 
gradient  field  of  the  dipole  are  then  expressed  by  the  relation 

where  6(})^  = 61})^  + 61})^  . 

Consequently,  small  displacements  in  the  steady  ambient 
field  with  the  cancelling  dipole  in  position  give  gradient 
changes  expressed  by  the  sum 

6923  + <Sg^3  = (6<p^-  6(})^) g^sinx^  . 


(28) 


(29) 
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Although  angular  deviations  owing  to  actual  rotation  of  the 
gradiometer  are  the  same  for  both  the  equivalent  dipole  and 
the  cancelling  dipole  and  so  cancel,  angular  deviations  owing 
to  virtual  rotations  accompanying  displacement  differ  to  the 
extent  that  curvature  of  the  ambient  field  deviates  locally 
from  curvature  of  the  cancelling  dipole  field. 

Because  we  can  position  a dipole  used  to  cancel  gradients 
at  a point  at  different  distances  from  the  point  by  adjusting 
its  moment*,  we  can  adjust  the  difference  in  angular  deviations 
owing  to  virtual  rotations  accompanying  small  displacements. 

Wb  place  the  dipole  close  to  the  point  of  cancellation,  ir  the 

equivalent  dipole  representing  the  ambient  gradient  field 
rotates  markedly  during  a small  displacement,  and  far  from 
the  point,  if  it  rotates  slightly.  Namely,  we  choose  the 
distance  from  the  point  of  cancellation  to  make  curvatures 
of  the  dipole  and  ambient  magnetic  fields  comparable. 

Although  initial  misalignment  of  the  cancelling  dipole 
gives  gradient  fluctuations  of  second  order  in  angular  devia- 
tions and  so  is  negligible  for  reasonably  precise  alignment, 
small  displacements  and/or  rotations  of  the  cancelling  dipole 
itself  relative  to  the  ambient  field  produce  additional  gra- 
dient fluctuations  of  first  order  and  so  degrade  its  effective- 
ness. Fluctuations  of  the  moment  of  the  cancelling  dipole, 
however,  produce  gradient  fluctuations  of  second  order.  To 
be  effective  than,  positioning  of  a dipole  to  cancel  ambient 
transverse  gradients  must  be  steady  with  respect  to  the  ambient 
field  and  reasonably  precise,  but  its  moment  may  fluctuate. 
Ideally,  a cancelling  dipole  should  be  rigidly  fixed  to  a 
magnetic  object  producing  a local  gradient  field. 


* The  ratio  m/r 


4 


is  constant. 
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B.  SUPPRESSION  OF  NOISE  FROM  FLUCTUATING  MAGNETIZATION 


In  the  absence  of  magnetic  objects,  the  ambient  magnetic 
field.  B(t),  is  effectively  uniform  and  comprises  the  earth's 
steady  magnetic  field,  ^ , and  a small  fluctuating  magnetic 
field,  ’b(t),  resulting  from  electric  currents  in  the  iono- 
sphere. A magnetic  object  distorts  the  ambient  field  and  so 
produces  steady  gradients  of  the  magnetic  field  in  its  vici- 
nity as  well  as  gradients  that  fluctuate  as  the  uniform  am- 
bient field  fluctuates.  An  object  having  a large  magnetic 
permeability,  for  example,  warps  the  surrounding  magnetic 
field  until  resulting  field  lines  are  nearly  normal  to  its 
surface . 


Gradients  of  the  magnetic  field  in  the  vicinity  of  a 
stationary  magnetic  object  fluctuate  as  the  undi&turh'd  am- 
bient field,  which  magnetizes  the  object,  chan o vs  in  strength 
and  direction  in  response  to  ionosphere  curre.its.  Magni- 
tudes of  gradients  are  proportional  to  the  strength  of  the 
ambient  field  divided  by  a characteristic  length  o*‘  -in  object 
and  decrease  with  distance  from  an  object  as  a of  its 

characteristic  length  divided  by  distance. 


For  example,  a spherical  iron  shell  magnetized  by  a uni- 
form magnetic  field,  B(t) , produces  an  induced  magnetic  field 
outside  the  shell  that  is  described  by  the  field  of  a magnetic 
dipole  located  at  the  center  of  the  sphere.  The  dioole  points 
along  the  uniform  magnetizinq  field,  and  its  moment  m(t)  is 
given  bv  the  expres.sion 


iTi(t) 


R^6(t) 

o 


9 


(30a) 
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where 


6 = 3 + X 


(Rl/R^) 


l+!3/2x) 


Here, 
and  X its 
bility  of 


is  the  outer  radius  of  the  shell;  R^,  its 
maanetic  susceptibility.  Because  maqnet 
iron  is  large  (x~1000)  , Tfi(t)  - 4TrR^B(t). 


inner  radius 
ic  suscepti- 


Gradients  of  the  magnetic  field  surrc unding  an  iron  sphere, 
then,  are  everywhere  described  by  the  gradient  field  of  a mag- 
netic dipole  that  is  fixed  at  the  center  of  the  sphere  and 
changes  its  direction  and  moment  in  response  to  changes  in 
direction  and  strength  of  the  uniform  ambient  magnetic  field. 
The  equivalent  dipole  representing  gradients  at  each  point  of 
the  field  is  the  actual  dipole  and  so  lies  along  a radial 
vector  emanating  from  the  center  of  the  sphere. 


Consequently,  we  find  from  Equation  13  that  the  transverse 
gradient  measured  by  the  gradiometer  at  a position  defined  by 
a radial  vector  r is  expressed  as 

g23(r,t)  = g ( t)  I siJix  (t)  j^cosQ  ( t)  cos(}>cos0  - sinfi  ( t)  sin()>cos2e 
- -|cosx  ( t)  sin({)sin20  I , 


where 


X(t)  and  f2(t)  specify  direction  of  the  dipole  axis  with  respect 
to  the  radial  vector  f , and  and  0 s|5ecify  direction  of  the 
radial  vector  in  the  gradiometer  basis.  Equation  31a  tells  us 


(30b) 


(31a) 


(31b) 


42 


that  aligning  the  axis  of  the  gradiometer  basis  along  a 
radial  vector  (<})  = 0 and  0 = tt/2)  gives  a null  response  what- 
ever the  dipole  orientation  an  ‘‘o  completely  suppresses 
noise  from  fluctuating  magne  . ;ion  currents  in  the  sphere. 
Slight  misalignment  of  the  gradiometer , however,  gives  the 
response 

Ag23  "fft)  = g (t)  sinx  ( t)  ^A<j)sinQ  ( t) -AGccsQ  ( t)|  , 

where  A<1)  and  A6  specify  misalignment  of  the  axis  and  radial 
vector  r. 

Nonetheless,  the  fluctuating  part  of  the  earth's  magnetic 
field  is  much  smaller  than  the  steady  part,  so  that 


g(t)  = + (5g(t),  with  (Sg(t)/g^  = b//(t)/B^ 


X(t)  = Xq  + «5x(t),  with  6x(t)  = 


X B ) (B^  X b,  (t)  ) 

o o 


B sinx 
o ^o 


B 


and 


fi(t)  = + 6Q(t),  with  6n(t)sinxQ  = 


(?  X S ) "6,  (t) 

O 


BoSinXo 


B 


o 


Here,  x and  fi  give  the  direction  of  B with  respect  to  r, 
^o  o o 

' - (--)'  ' 

R \ r / 


% = 


o 


b,,(t)  = 


B^.b(t) 


B‘ 


®o  ' 


(31c) 


(32a) 

(32b) 

(32c) 

(32d) 

(32e) 
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which  is  the  component  of  the  fluctuating  part  of  the  field 

parallel  to  B , and 
^ o 


bj_(t) 


o 


X to  X B ) 
o 


t 


which  is  the  component  perpendicular  to  Consequently, 

allowing  for  slight  misalignment.  Equation  31c  shows  that  gra- 
dient fluctuations  resulting  from  fluctuating  magnetization 
currents  in  an  iron  sphere  are  of  second  order  and  so  negli- 
gible when  the  axis  of  the  gradiometer  is  aligned  along  a 
radial  vector. 


For  an  intricately  shaped  object,  we  surmise  that  an  equi- 
valent dipole  describing  gradients  at  a point  near  the  object 
changes  not  only  its  orientation  and  moment  but  also  its  loca- 
tion as  strength  and  direction  of  the  earth's  magnetic  field 
fluctuates.  Gradient  fluctuations  resulting  from  fluctuating 
magnetization  currents  are  then,  in  effect,  equivalent  to 
fluctuations  resulting  from  irregular  motions.  Fluctuating 
magnetization  currents  in  an  intricately  shaped  magnetic 
object  produce  virtual  rotation  of  an  equivalent  dipole. 

Aligning  the  gradiometer  so  that  the  equivalent  dipole 
determined  from  the  steady  gradient  field  of  a magnetic  object 
lies  along  the  axis  of  the  gradiometer  basis  with  its  axis 
in  the  plane  containing  pickup  loops  nullifies  gradient  fluctua- 
tions resulting  from  fluctuations  of  the  dipole  orientation 
and  moment  produced  by  fluctuations  of  the  earth's  magnetic 
field.  Virtual  rotations  of  the  equivalent  dipole  produced  by 
f I'^tituat ions  of  the  earth's  magnetic  field  then  give  gradient 
changes  expressed  by  the  relation 


(32f) 
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(33) 


'5go-j  = 64)  g sinv  1 
^23  vrn^o  '^o  ‘ 

where  64;^  is  the  apparent  angle  of  rotation  about  the  gra- 
diometer  axis. 
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Section  VII 

PRELIMINARY  TESTS  OF  NOISE  SUPPRESSION  TECHNIQUES 

Our  field  tests*  of  noise  suppression  techniques  are 
necessarily  preliminary  because  the  waoden  gimbal  used  to 
orient  the  gradiometer  provides  only  vertical  and  horizontal 
axes  of  rotation  and  so  both  restricts  gradiometer  orienta- 
tion (three  axes  of  rotation  provide  unrestricted  orientation) 
and  limits  precision  of  alignment  to  about  0.1  radian  or  a 
few  degrees.  We  report  results  of  tests  of  two  techniques: 

(I'l  use  of  a current  in  a coil  approximating  a magnetic  di- 
pole to  cancel  the  steady  gradient  field  at  a point  near  a 
magnetic  object  and  (2)  aligning  the  gradiometer  to  suppress 
gradient  fluctuations  from  fluctuating  magnetization  currents 
in  an  iron  sphere.  In  each  case,  we  first  summarily  describe 
the  technique,  give  a detailed  procedure,  and  then  present 
results  of  preliminary  tests. 

A.  CANCELLING  STEADY  GRADIENTS 


To  cancel  steady  gradients  at  a point  near  a magnetic 
object,  we  place  a coil  at  the  position  of  an  equivalent 
dipole  corresponding  to  gradients  at  the  point  and  adjust 
current  in  the  coil  to  nullify  the  gradients.  We  determine 
positions  of  equivalent  dipoles  by  measuring  gradients  in  a 
fiducial  basis,  calculating  corresponding  eigenvalues  and 
eigenvectors,  and  transforming  to  a description  in  terms  of 
locati.on,  orientation,  and  moin<=>nt  of  associated  equivalent 
dipoles.  After  placing  the  coil  at  a position  determined 
from  gradients  measured  in  the  fiducial  basis,  we  drive  an 
cscillating  current  in  the  coil  and  adjust  its  location  and 

* We  do  tests  at  the  La  Posta  Astrogeophysical  Observatory 
located  in  the  mountains  about  70  miles  east  of  San  Diego, 
California.  Coordinates  of  the  observatory  are  116°  25'  6" 
west  longitude  and  32°  40'  39"  north  geodetic  latitude 
(Bleiweiss  and  Wefer,  lO'^S;  . 
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orientation  to  obtain  a null  response  at  the  r. radiometer  with 
the  gradiometer  basis  aligned  along  principal  axes  of  the 
ambient  gradient  field.  We  >hen  tilt  the  gradiometer  away 
from  principal  axes  and  adjust  a steady  current  in  the  coil 
until  rotating  the  gradiometer  about  its  axis  gives  a null 
response . 

1.  Procedure 


To  determine  the  five  independent  steady  gradients  Yj^j 
defined  in  a fiducial  basis  |z^[  , we  rotate  the  gradiometer 
about  each  one  of  three  axes  forming  a linearly  independent 
triad.  Rotating  the  gradiometer  through  an  angle  about  its 
axis  gives  the  sinusoidal  response  expressed  by  the  relation 


g22  = U cos  Q - V sin  (34a) 

Amplitudes  U and  V depend  on  initial  orientation  of  the  gradioreter 
basis.  Measuring  amplitude  and  phase  of  the  siiiusoidal  re- 
sponse developed  during  a complete  rotation  determines  the 
coefficients  U and  V corresponding  to  an  initial  orientation. 

For  the  three  axes  of  rotation,  we  choose  the  axis  of 
the  fiducial  basis  and  two  axes  defined  by  the  two  sets  of 
polar  angles  (0,tt/2)  and  (0,0)  with  respect  to  the  fiducial 
basis.  When  the  gradiometer  basis  is  initially  coincident 
with  the  fiducial  basis,  coefficients  U and  V are  equal  to  the 
gradients  Y23  ^13'  ^®sP®<^tively . Rotating  the  gradiometer 

about  its  axis  then,  in  effect,  measures  the  gradients  Y23 
Yj^2'  When  the  axis  of  the  gradiometer  is  aligned  with  an  axis 
of  rotation  defined  by  polar  angles  0 and  tt/2,  the  coefficients 
are  given  by  the  relations 

U = Y23<^°s0  - y(’Y22"Y33)  sin20  (34b) 
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and 


Rotating  the  gradiometer  then  gives  a measure  of  y the 

difference  (Y22“"^33^  ' Y23  ^^13  determined.  Finally, 

when  the  gradiometer  axis  is  aligned  along  an  axis  of  rotation 
defined  by  polar  angles  0 and  0,  the  coefficients  are  given 
by  the  relations 


U = Y23*^°s®  t Y2^2^^^* 


and 


V = 


2 (Y11-Y33) sin20 


Rotating  the  gradiometer  then  gives  a measure  of  the  difference 
(Yi3~Y33)»  once  Y23'  "*^13'  ^12  determined.  Because  the 

sum  of  longitudinal  gradients  vanishes,  the  three  rotations  give 
a measurement  of  the  independent  gradients  Y^^^  "^33'  "*^12'  "^13' 
and  Y23  defined  in  the  fiducial  basis. 


From  the  five  gradients,  we  determine  three  eigenvalues 

^o^^^spo^ding  eigenvectors  'e^,  “e^)  , and 

three  Euler  angles  (4>,  G,  T)  that  specify  a rotation  R (<I>, 

P 

0,  4')  from  a fiducial  basis  to  the  basis  comprised  of  prin- 
cipal axes  defined  by  the  eigenvectors*.  Namely,  we  establish 
the  relation 


G = R Air 
P P 

thcit  gives  the  matrix  G describing  gradients  in  the  fiducial 
basis  in  terms  of  Rp  and  the  diagonal  matrix  A whose  elements 
are  the  eigenvalues  X^,  and  X^.  Gradients  in  the  fiducial 

basis,  however,  are  also  equal  to  gradients  of  the  field  of  a 
magnetic  dipole. 

* The  Appendix  delineates  calculation  of  eigenvalues,  eigen- 
vectors, and  Euler  angles  from  gradients  of  an  ambient  field. 
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(34c) 


(34d) 


(34e) 


(35) 


The  relation 

G = R(({),  0,  f^)T(a)  AY  (a)'R(ct),  0, 


gives  gradients  of  the  field  of  a magnetic  dipole  located 
with  respect  to  a fiducial  b.*sis  by  polar  angles  (j>  and  0 
defining  direction  of  its  position  vector  F and  oriented  with 
respect  to  the  position  vector  by  the  azimuthal  angle  Q and 
polar  angle  x,  as  defined  in  Figure  4.  The  matrix  R(<}),  0, 
represents  three  consecutive  rotations,  as  defined  by  Equation 
11a,  that  bring  the  fiducial  basis  into  coincidence  with  the 


basis 


A 


defined  in  Fiaure  4,  and 


Y(a)  = 


cosa 

0 

sina 


0 

1 

0 


-sina 

0 

cosa 


A 


and  represents  a rotation  through  an  angle  -a  about  the  axis 
that  brings  the  basis  into  coincidence  with  principal  axes 

of  the  gradient  field  of  the  dipole  when  tan  2a  = (2/3)  tany . 
Eigenvalues  of  the  gradient  field  of  a dipole  equal  eigenvalues 
of  an  ambient  gradient  field  when 


tanx 


J 


^3] 


0<x<fT 


and 


4irr^ 


^2  ^1  ^3, 


1/2 


where  m is  the  moment  of  the  dipole  and  r,  its  distance  from 
the  field  point. 


(35a) 


(35b) 


(35c) 


(35d) 


48 


Consequently,  Equations  35  and  35a  tell  us  that  gradients 
of  a dipole  field  equal  gradients  of  an  ambient  field  when 
location  ((j>,  6),  orientation  (x,  / and  moment  of  the  dipole 

are  determined  by  the  relation 

R (<5>,  0,  'i')  = P((|),  0,  ro'Y(a) 

IT 

together  with  Equations  35c  and  35d.  Polar  angles  d>  and  0 

m -.-n 

specifying  direction  of  the  dipole  axis  in  the  fiducial  basis 
are  determined  by  the  relation 

R(4),  0,  fJ)Y(x)  = Z(4>^  + tt/2)X(0j^), 

where  matrices  Z(ip)  and  X(ip)  are  given  by  Equations  11b  and 
11c. 


From  Equations  35a  and  36a,  we  then  find  that  gradients  of 
the  field  of  an  arbitrarily^  positioned  dipole  are  expressed  in 
the  basis  comprised  of  principal  axes  of  an  ambient  gradient 
field  by  the  relation 


r = Y 'R  RYA -YRR  Y 
e e d e e 


f 


where  elements  of  the  diagonal  matrix  are  eigenvalues  of  the 
gradient  field  of  the  dipole  and  siibscript  e marks  matrices 
corresponding  to  an  equivalent  dipole.  If  a dipole  is  near 
the  position  of  an  equivalent  dipole,  then  R = RgR^^  / Y =Y^Y^ , 
and  A^  = (g/g^)(A  + 6A)  , so  that 


r = (g/g  )Y  R.Y^Y  (A  + 6A)Y  Y^R^Y^ 
e e66e  eo6e 


(36a) 


(36d) 


(37a) 


(37b) 
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where  the  ratio  g/g^  gives  magnitude  of  gradients  o^  the  di- 
pole field  relative  to  ambient  gradients  and  6A  is  proportional 

to  5x- 

The  matrix  product  infinitesimal 

rotation  that  we  express  in  terms  of  angular  deviations  in 
location  and  orientation  of  the  dipole  by  che  sum 

VsVe  ‘ ^ ’'♦y  "y  “♦z  ' 

where 

= 6<{)(cose^sinag+sine^sinQgCOsa^) 

+ 66cosQ  cosa  + Si^sina  / 

6 e ^ 

64^^  = 64>sin6^sinQg  - 66sinQg-6a  , 

54^2  ~ 64)  (cosOgCOsa^  - sinB^sinQ^sina^) 

-60cosQ  sina_  + 6Qcosa  , 
e e e 

elementary  matrices  e^,  and  are  given  by  Equation  23d, 

and  we  find  from  Equation  36d  that 

6Q  = (64)j^-64))c;os0^+6  0j^sin(4>g-4>n,e)sin0g 

and 

6x  = 6 0j^  jcos0gcosn^sin(4)j^^-4>g)-sinfi^cos(4)^g-4)g)j 
-60sinf2g  + (64)j^-64>)  sinB^cosSi^ 


(37c) 

(37d) 

(37e) 

(37f) 

(37g) 

(37h) 
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Consequently,  we  find  from  Equation  37b  that 


r = (g/g^)  (A  + 6 A + 6T) 


(38a) 


and,  to  first  order  in  angular  deviations,  that  the  relation 


(g/g^) 6T  = (g/g^) 


(A^-A2)6ti^^  ( A^-A  6ijjy 

0 ( A2~A2) 


(A3-A^)6iJ;y  (X^-X^)6ip^ 


(38b) 


gives  transverse  gradients  in  the  basis  comprised  of  principal 
axes  of  an  ambient  gradient  field.  Specifically,  the  exores- 
sion 


^23^  (g/ge^ 

gives  response  of  the  gradiometer  when  the  gradiometer 
basis  coincides  with  the  basis  comprised  of  principal  axes 
of  an  ambient  gradient  field,  and  the  expression 

^23  ~ [ (A2-A3)  6i(;^cost|;  - ( A^-A^^)  6iJ;^simJ;j  (38d) 

gives  response  when  the  gradiometer  is  rotated  through  an 
angle  ij;  about  its  axis.  Equation  38b  tells  us  that  true  align- 
ment of  a dipole  at  the  position  of  an  equivalent  dipole  gives 
a null  response  whatever  the  dipole  moment.  Moreover,  large 
dipole  moments  provide  means  of  discerning  minute  misalignment. 


I 
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To  position  and  align  the  moment  of  a coil  approximating 
a magnetic  dipole  precisely  at  the  location  and  orientation 
of  an  equivalent  dipole  of  an  ambient  gradient  field,  then, 
we  first  align  the  gradiometer  basis  coincident  with  the 
basis  comprised  of  principal  axes  of  the  ambient  gradient 
field,  as  described  subsequently,  and  place  the  coil  at  the 
location  and  orientation  of  an  equivalent  dipole  determined 
from  measurements  of  gradients  in  a fiducial  basis.  We  then 
drive  the  coil  with  current  oscillating  at  a frequency  of 
say  5 Hz  and  adjust  location  and  orientation  of  the  coil  un- 
til the  peak  at  5 Hz  in  the  spectrum  of  gradiometer  response 
(as  displayed  by  a real-time  spectrum  analyzer) ' falls  below 
the  noise  level  of  the  instrument.  The  plane  defined 
by  the  dipole  axis,  m,  and  position  vector,  r,  is  then  nor- 
mal to  the  eigenvector  ^2-  Next,  we  rotate  the  gradiometer 
90  degrees  about  its  axis,  so  that  pickup  loops  face  the 
dipole,  and  again  adjust  location  of  the  coil  until  the  spec- 
tral peak  falls  below  the  level  of  instrument  noise.  Princi- 
pal axes  of  the  gradient  field  of  the  dipole  are  then  coinci- 
dent with  principal  axes  of  the  ambient  gradient  field,  so 
that  the  matrix  6T  vanishes,  but  the  matrix  6A  does  not  neces- 
sarily vanish.  Finally,  we  return  the  gradiometer  to  the 
fiducial  basis  and  adjust  a steady  current  in  the  coil  until 
rotating  the  gradiometer  about  its  axis  gives  a minimum  re- 
sponse so  that  g/9g  = “1-  The  response  during  rotation  is 
then  proportional  to 

To  align  the  gradiometer  basis  precisely  coincident  with 
the  basis  comprised  of  principal  axes  of  an  ambient  gradient 
field,  we  initially  align  the  basis  with  principal  axes  de- 
termined from  measurements  of  gradients  in  a fiducial  basis, 
so  that  the  matrix 

r = A - 6t 


(39a) 
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represents  gradients  in  the  gradiometer  basis,  where  elements 
of  the  diagonal  matrix  A are  eigenvalues  of  the  ambient  gra- 
dient field  and 


6 T 


0 


Uj-Xi)6ty 


( 6t 

( A2-A2)  <5t 


y 

X 


0 


Angular  deviations  are  given  by  the  relations 


6t 

= 64>sin0sin4'  + 

60COSH' 

X 

6t 

= 64>sin0cos4'  - 

60sinH' 

y 

and 


6t  = 6<I>cos0  + 6^ 
z 

where  641,  60,  and  6H'  are  deviations  in  the  Euler  angles 
4),  0,  and  4^  that  define  a rotation  from  the  fiducial  basis 
to  the  basis  comprised  of  principal  axes  of  the  ambient  gra- 
dient field. 

Rotating  the  gradiometer  about  its  axis  through  an  angle 
\p  then  gives  the  response  expressed  by  the  relation 

g23(4/)  = (A3-A2)  6t^cosi|x  - ( A^^-A^)  6TySint|j 

Equation  40a  tells  us  that  true  alignment  of  the  gradiometer 
axis  with  the  eigenvector  ^3  gives  a null  response  during  rota- 
tion of  the  gradiometer  about  its  axis,  so  we  adjust  angular 
deviations  64>  and  60  until  rotation  gives  a null  response. 


(39b) 

(39c) 

(39d) 

(39e) 


(40a) 
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Similarly,  rotating  aboat  the  ^2  of  the  gradiometer  basis 

gives  the  response  expressed  by  the  relation 

~ (ST^cosijj  + ( A^-A  6T^sinip 

Equation  40b  tells  us  that  true  alignment  of  the  X2  axis  of 
the  gradiometer  basis  with  the  eigenvector  ^2  gives  a null 
response  during  rotation  and  so  affords  means  of  adjusting 
the  angular  deviation  64'. 

2 . Results  of  Preliminary  Tests 

The  wooden  gimbal  used  to  orient  the  gradiometer  provides 
only  vertical  and  horizontal  axes  of  rotation  and  so  pre- 
cludes aligning  the  gradiometer  basis  coincident  with  the 
basis  comprised  of  principal  axes  of  an  ambient  gradient 
field.  Although  we  can  not  use  the  present  gimbal  to  test 
thoroughly  the  procedure  outlined  above  for  cancelling  am- 
bient gradients,  our  preliminary  tests  show  that  a steady 
current  in  a coil  placed  approximately  in  the  position  of 
an  equivalent  dipole  reduces  ambient  gradients  by  about  a 
factor  of  100. 

To  facilitate  positioning  the  coil,  we  first  chose  to  cancel 

ambient  gradients  resulting  from  an  iron  sphere,  which  was  placed 

4.5m  north  of  the  gradiometer  and  has  a radius  of  0.74  m.  Magni 

4 

tude  of  the  earth's  magnetic  field  is  about  4.5  X 10  nT  at 
the  test  site,  so  that  the  dipole  induced  at  the  center  of 
the  sphere  by  the  earth's  magnetic  field  produces  gradients 
having  magnitudes  proportional  to 

4 

= 133  nT/m 


(40b) 
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The  relation 


^23  ~ ^ ^sinxcosSlcosOcoscj)- (ainxsinncos20  + |cosxsin29)  sin<})j  (41a) 

gives  the  gradiometer  response  to  the  induced  dipole,  where 
polar  angles  (|)  and  0 c -ecify  direction,  in  the  gradiometer 
basis,  of  the  position  vector  Ic  pointing  from  the  dipole  to 
the  gradiometer,  as  shown  in  Figure  4,  and  angles  x n 

define  direction  of  the  dipole  axis,  which  is  parallel  to  the 
earth's  magnetic  field,  with  respect  to  the  position  vector. 

Because  of  rough  terrain  at  the  test  site,  the  center  of  the 
sphere  is  somewhat  above  the  center  of  the  gradiometer,  but 
we  place  the  sphere  so  that  the  position  vector  and  gradio- 
meter axis  lie  approximately  in  a vertical  plane  containing 
the  e^'.rth's  magnetic  field  vector.  Then  = 3tt/2  , and  so 

g23  ^ ‘ } sin2  ( 9-a)  sin<^  , (41b) 

where  the  expression 

Aj^-A^  = g (4  + 5cos^X)  (41c) 

gives  the  difference  in  eigenvalues  A^^  and  A^  of  gradients  of  the 
induced  dipole  field  and  the  angle  a, determined  by  the  relation 
tan  2a  = (2/3)tanx,  gives  inclination  of  the  eigenvector  to  the 
position  vector,  as  shown  in  Figure  11, and  x = it- (<t)^+0^)  . 

The  local  dip  angle  of  the  earth's  magnetic  field,  4>^,  is 
about  55°  at  the  site,  and  the  dip  angle  of  the  position  vec- 
tor, 0j,  is  about  5°  or  so. 

Rotating  the  gradiometer  about  its  axis  (0<4)<2tt)  and 
measuring  amplitude  of  the  response  gives  a gradient  magnitude 
of  56  nT/m,  and  so  we  conclude  from  Equation  41b  that  0-a  ^ 11°, 
which  tells  us  that  the  gradiometer  axis  is  tilted  somewhat 


55 


from  vertical.  Because  alignment  procedures  at  the 
site  are  crude  at  present,  we  exoect  angular  deviations 
of  the  order  of  0.1  radians  or  a few  degrees. 

To  cancel  gradients  of  the  induced  dipole  field,  we  place 
a coil  about  5.6  m north  of  the  gradiometer  and  align  its 
axis  parallel  to  the  earth's  magnetic  field  vector.  The  coil 
comprises  1000  turns  of  wire  wound  on  an  aluminum  frame  about 
0.97  m in  diameter.  After  making  slight  adjustments  in  lateral 
position  and  orientation  of  the  coil,  we  find  that  a current 
of  0.6  A in  th  ) coil  reduces  amplitude  of  the  response  to  rota- 
tion of  the  gradiometer  about  its  axis  from  56  nT/m  to  0.49  nT/m, 

or  by  a factor  of  100.  For  a current  of  0.6  A,  the  moment  of 

2 

the  coil  is  444  A m , and  so 
3y  m 

g = = 135  nT/m  , 

^ 4TTr 

with  r = 5.6  m,  which  is  nearly  equal  to  the  strength  of  gra- 
dients produced  by  the  sphere,  g = 133  nT/m. 

Because  a magnetic  dipole  parallel  to  the  earth  s magnetic 
field  at  the  center  of  an  iron  sphere  describes  the  gradient 
field  surrounding  the  sphere,  placing  a coil  at  the  position 
of  an  equivalent  dipole  of  the  gradient  field  is  straightfor- 
ward. For  irregularly  shaped  magnetic  objects,  however,  de- 
termination of  location  and  orientation  of  equivalent  dipoles 
requires  measurement  of  gradients  in  a fiducial  basis.  As  a 
first  test  of  cancellation  techniques  for  irregularly  shaped 
objects,  we  placed  a collection  of  oil  drums  and  steel  clanks 
approximately  5 m north  of  the  gradiometer. 
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We  define  a fiducial  basis  by  rotating  the  gradiometer 
about  its  axis,  when  vertical,  until  a null  response  is 
obtained,  so  that  q22  vanishes  in  the  fiducial  basis.  To 
measure  gradients  in  the  fiducial  basis,  we  rotP'-e  the  gra- 
diometer first  through  an  angle  tj)  about  its  axis  (the 
axis  of  the  gradiometer  basis)  and  then  through  an  angle  6 
about  the  gradiometer  basis,  so  that  the 

expression 

^23  " [^12^°^^'^  ■ ■ Yi3sinc()cose 

gives  its  response  in  terms  of  gradients,  in  the  fiducial 

basis  and  angles  of  rotation.  Rotations  about  the  gradiometer 
axis,  0 = 0,  and  about  horizontal  axes  corresponding  to  t()  = 0 

and  (P  = tt/4  then  determine  gradients  y^3  and  y^2 
2y  + y23*  Next,  we  rotate  the  gradiometer  through  an  angle 
If;  about  the  x^  axis  of  the  gradiometer  basis,  so  that  the 
expression 

■323  ' + 2Y33)sin2* 

gives  its  response  in  terms  of  gradients  in  the  fiducial  basis 
and  the  angle  of  rotation  and  so  determines  the  sum  y^^  + 2y^3* 
Since  y2^  vanishes,  the  procedure  measures  the  four  nonvanishing 
gradients  in  the  fiducial  basis. 

From  gradients  measured  in  the  fiducial  basis,  we  find,  as 
described  in  the  Appendix,  location,  orientation,  and  gradient 
strength  g of  the  two  pairs  of  equivalent  dipoles  corresponding 
to  the  collection  of  magnetic  objects.  The  gradient  strength 
of  the  objects  is  about  67  nT/m,  and  the  plane  containing  axes 
of  equivalent  dipoles  is  no  longer  vertical  as  for  the  iron 
sphere  but  is  tilted  about  30°  from  vertical  about  the  x^^  axis 

of  the  fiducial  basis. 


(42a) 


(42b) 
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Because  the  gimbal  used  to  orient  the  coil  provides  only 
a horizontal  axis  of  rotation,  it  precludes  accurately  aligning 
the  axis  of  the  coil  in  a tilted  plane.  Crude  placement  of 
the  coil  at  the  position  of  an  equivalent  dipole,  however, 
reduced  gradients  by  about  a factor  of  30. 


B.  ALIGNING  FOR  NULL  RESPONSE 

To  suppress  gradient  fluctuations  resulting  from  fluctuating 
magnetization  currents  in  nearby  magnetic  objects,  we  orient 
the  gradiometer  so  that  an  equivalent  dipole  corresponding  to 
the  steady  gradient  field  of  the  objects  lies  along  the 
axis  of  the  gradiometer  basis  with  its  axis  in  the  plane  con- 
taining pickup  loops.  To  align  the  gradiometer,  we  begin  by 
following  tne  procedure  described  for  cancelling  steady  gra- 
dients. Namely,  we  align  the  gradiometer  basis  coincident 
with  the  basis  comprised  of  principal  axes  of  the  steady  gra- 
dient field  and  precisely  position  a coil  at  the  location  and 
orientation  of  an  equivalent  dipole. 

With  a current  oscillating  in  the  coil,  we  then  tilt  the 
gradiometer  about  the  X2  axis  of  the  gradiometer  basis  and 
rotate  it  about  its  axis  until  a maximum  response  is  attained 
so  that  pickup  loops  face  the  coil.  Next,  we  point  the  axis 
of  the  coil  directly  toward  the  gradiometer  and  turn  the  gra- 
diometer about  the  x^  axis  of  the  gradiometer  basis  until  a 
null  response  is  obtained.  The  axis  of  the  coil  then  lies 
along  the  X2  axis  of  the  gradiometer  basis. 

Finally,  we  return  the  axis  of  the  coil  to  its  initial 
position  and  turn  the  gradiometer  about  its  axis  until  a null 
response  is  obtained.  The  coil  then  lies  along  the  x^^  axis 
of  the  gradiometer  basis  with  its  axis  in  the  plane  containing 
pickup  loops.  response  tnen  vanishes  fur  every  orien- 

tation of  the  coil. 
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1.  Procedure 


To  align  the  gradiometer  for  a null  response  from  a current 
oscillating  in  a coil  positioned  at  the  location  and  orienta- 
tion of  an  equivalent  dipole,  we  first  align  the  gradiometer 
basis  coincident  with  the  basis  comprised  of  principal  axes 
and  then  turn  the  gradiometer  through  an  angle  tt/2  - a about 
the  X2  axis  of  the  gradiometer  basis,  so  that  the  coil  then 
lies  along  the  x^^  axis  of  the  gradiometer  basis  with  its  axis 
in  the  plane  containing  pickup  loops.  The  matrix 

R {())  , 0 , i|;  ) = R (4>,  0,  H')y(TT/2  - a) 
a a a a p 

represents  consecutive  rotations  of  the  gradiometer  from  a 
fiducial  basis  first  to  the  basis  comprised  of  principal  axes 
(represented  by  R^)  and  then  about  the  axis  of  the  gradio- 
meter basis  through  the  angle  tt/2  - a.  Angles  <I> , 0,  and  T 
are  Euler  angles  specifying  a rotation  to  the  basis  comprised 
of  principal  axes,  and  the  angles  6^,  and  are  Euler 

angles  specifying  the  rotation  that  aligns  the  gradiometer 
for  a null  response;  namely, 

R (4)  , 0 , ) = Z((}).  + Tr/2)X(0^)Z(i|;^) 

a a a a c a a 

From  Equation  36a  that  specifies  location  and  orientation 
of  an  equivalent  dipole  in  terms  of  Euler  angles  defining  a 
rotation  to  a basis  comprised  of  principal  axes,  we  then  find 
that 


R {(})  , 0 , 4;  ) = R{(}),  0,  n)Y(n/2)  , 

a a a a 

so  that  the  relations 

COS0  = sinOcosJ^ 

a 


(43a) 


(43b) 


(44a) 


(44b) 
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(44c) 


sin0  sin((})  - ({>  ) 

3.  3 

sin0  cos((()  - (})  ) 

d 3 

cosip  = sln0sin((J) 

d 


-cos0cosf2 


and 


sinij;  = cos9sin0  - sin9cos0  cos((J)  - <p  ) 
a a a a 

give  Euler  angles  <{)  , 9 , and  tjj  in  terms  of  angles  , 9 , and 

3 

that  specify  location  and  orientation  of  an  equivalent  dipole. 


When  the  gradiometer  is  aligned  for  a null  response,  gra- 
dients are  given  by  the  relation 


r 


R GR 
a 


a 


wnere  the  relation 


G = R((j),  9,  n)Gy  (g,x)R(4>» 

with  Gy(g,x)  given  by  Equation  10a,  gives  gradients  in  the 
fiducial  basis  resulting  from  current  in  the  coil  positioned 
at  the  location  and  orientation  of  an  equivalent  dipole.  For 
true  alignment,  then. 


= Y(Ti/2)Gy(g,x)Y(TT/2)  = g 


-2cosx  0 

0 cosx 

.-sinx  0 


-sinx  ' 
0 

cosx  _ 


where  x is  the  polar  angle  giving  inclination  of  the  acis  of 
the  coil  to  its  position  vector.  Equation  45c  shows  that  gra- 


(44d) 

(44e) 

(44f) 


(45a) 

(45b) 


(45c) 
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dients  are  then  independent  of  the  azimuthal  angle  and 

that  the  qradiometer  response  vanishes  whatever  the  polar  angle 

X and  so  vanishes  for  every  orientation  of  the  coil. 


For  misalignment  of  the  gradioiteter , ' where 

represents  an  infinitesimal  rotation.  We  then  find  that 


r - a5T 
o - a 


where 


- 

• 

p 

0 sinx 

0 

-2sinx 

0 

3cosx 

6T  =6A 
a X 

sinx  0 

0 

0 

0 

0 

-6Az 

. 0 0 

0 . 

. 3cosx 

0 

2sinx. 

L 

0 3cosx 
3cosx  0 
0 sinx 


5A  = 6())sin6  sini|)  + 59cosi1j  , 
X a a a 


6A  = 5d)sin6  cosiIj  “ 6esint|^  , 
y ^ a a a 


and 


6A  = 5(})Cos9  + 6ip 

Z 5L 

To  make  adjustments  of  6(j>  and  69,  then,  first  we  orient 
the  gradiometer  for  a null  response  and  turn  it  about  iLs  axis 
to  obtain  a maximum  response  from  current  oscillating  in  the 
coil,  so  that  picxup  loops  face  the  coil.  Its  response  is 
then  given  by 

^23  “ ^ [“sinx  + (5(j)sin9^cos<j;^  - 69sin4)^)  3cosxj 

We  then  point  the  axis  of  the  coil  directly  at  the  gradiometer 
so  that  X ^ 0 adjust  6(j)  and  60  to  obtain  a null  response. 


(46a) 

0 

sinx 

(46b 

0 

(46c) 

(46d) 

(46e) 


(47) 
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The  axis  of  the  coil  then  lies  along  the  axis  of  the  gra- 
diometer  basis. 

To  adjust  6i1j,  we  return  the  axis  of  the  coil  to  its  initial 
orientation  and  turn  the  gradiometer  about  its  axis  to  obtain 
a null  response.  The  coil  then  lies  along  the  x^  axis  of  the 
gradiometer  basis  with  its  axis  in  the  plane  containing  pickup 
loops.  Finai.ly,  we  point  the  axis  of  the  coil  in  several  di- 
rections to  ensure  that  the  response  is  null  for  every  orienta- 
tion of  the  coil. 

2 . Results  of  Preliminary  Tests 

Although  we  can  not  thoroughly  test  the  alignment  procedure 
outlined  above  for  suppressing  noise  from  fluctuating  magneti- 
zation currents  in  magnetic  objects  because  the  gimbal  used  to 
orient  the  gradiometer  precludes  making  the  required  alignments 
precisely,  our  preliminary  tests  show  that  crude  alignment  of 
the  gradiometer  reduces  noise  power  from  fluctuating  magneti- 
zation currents  in  an  iron  sphere  by  a factor  of  about  100. 

As  a first  test  of  noise  suppression  afforded  by  aligning 
a gradiometer  for  a null  response,  we  placed  an  iron  sphere, 
having  a radius  of  0.46  m,  at  distances  of  3.7  m and  1.8  m 
north  of  the  gradiometer.  At  each  distance,  we  place  the 
sphere  so  that  the  axis  of  the  gradiometer  and  the  position 
vector  directed  from  the  center  of  the  sphere  to  the  gradio- 
meior  lie  approximately  in  a vertical  plane  containing  the 
earth's  magnetic  field  vector. 

The  magnetic  dipole  induced  at  the  center  of  the  sphere 
by  the  earth's  magnetic  field  then  gives  a gradiometer  response 
described  by  the  relation 

^23  5cos^X)^^^sin2(6  - a)sin4> 


(48a) 
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where  (p  and  6 are  polar  angles  specifying  direction  of  the 
position  vector  in  the  gradiomater  basis;  the  relation  tan  2a  = 
(2/3)tanx  determines  the  angle  a;  and  x ' 

v/here  (j> , is  the  local  dip  angle  of  the  earth  s magnetic  field 

Q 

and  0^,  the  dip  angle  of  the  position  vector.  The  relation 


gives  gradient  strength  of  the  induced  dipole,  where  B = 4.5  X 
lo"^  nT  is  the  strength  of  the  earth's  magnetic  field  at  the 
site,  R is  the  radius  of  the  sphere,  and  r,  its  distance  from 
the  gradiometer.  At  the  distance  r = 3.7m,  g = 70  nT/m,  and  at 
r = i.8m,  g = 1252  nT/m. 

Rotating  the  gradioraeter  about  its  axis  and  measuring  am- 
plitude of  the  response  gives  a gradient  magnitude  of  7.8  nT/m, 
with  r = 3.7m,  and  a magnitude  of  136  nT/m,  with  r = 1.8m. 

The  center  of  the  sphere  lies  somewhat  below  the  gra- 
diometer in  each  case,  9^  - ~5*^,  and  so  for  (})^  = 55  , we  con- 
clude from  Equation  48a  that  in  each  case  9 - a = 3°,  which 
again  tells  us  that  the  gradiometer  axis  is  tilted  somewhat 
from  vertical  and  that  9 = tt/2  - A9  with  A9  = 10  . 

To  examine  the  affect  of  gradiometer  alignment  on  suppressing 
noise  from  fluctuations  in  strength  and  direction  of  the  dipole 
induced  in  the  sphere  that  result  from  fluctuations  in  strength 
and  direction  of  the  earth's  magnetic  field,  we  express  gradio- 
meter response  in  terms  of  components  of  the  fluctuating  di- 
pole in  the  gradiometer  basis.  Namely,  we  write  the  fluctuating 
gradiometer  response  as 

6g23(t)  = 6g^(t)  (|)(cos29  -1)  cos9sin2(l) 

+ 6g2(t)  (j)  |^5cos20  - 1 - 5(cos29  - 1)  cos2<J)J  cos9 
-6g3(t)  (I)  (5cos29  + 3)sin9sin<J)  , 


(48b) 


(49a) 
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where 


6g^(t)  = g 

and  b^(t)  are  components  of  the  fluctuating  part  of  the  earth's 
magnetic  field  in  the  gradiometer  basis.  For  0 = tt/2  - A0, 
then.  Equation  49a  tells  us  that 

<5g23(t)  = (Sg2(t)sin4)  + ^ [6gj^  (t)  5sin2(|)  + (Sg2  ( t ) ( 3-5cos24>  ;J 

Moreover,  Equation  48a  shows  that  rotating  the  gradiometer 
about  its  axis  until  response  to  steady  gradients  vanishes 
gives  the  zero  position  for  the  angle  (|) . For  4>  = 0 , response 
to  gradient  fluctuations  is  suppressed;  namely, 

(5g23{t)  = -<5g2{t)A0  , 

and  for  4>  = tt/2,  response  to  gradient  fluctuation  is  unsup- 
pressed; namely, 

<5g23(t)  = (5g^{t)  + 46g2(t)A0 

To  investigate  the  effect  of  suppressing  noise  by  aligning 
a gradiometer  for  a null  response,  then,  we  compare  spectra 
of  gradiometer  response  for  the  two  alignments  4)  = 0 and  4>  = 
tt/2  with  the  iron  sphere  placed  1.8m  north  of  the  gradiometer. 
Figure  12a  shows  a spectrum  obtained  for  the  alignment  4>  = tt/2, 
and  Figure  12b  shows  a spectrum  obtained  the  following  day  for 
the  alignment  4-  = 0.  The  heavy  solid  line  in  each  figure  de- 
lineates the  spectrum  of  inherent  instrument  noise  observed  in 
a magnetically  quiet  environment,  as  shown  in  Figure  8 . 

Heavy  broken  lines  delineate  the  spectrum  expected  without 


b^(t) 


B 


9 


(49b) 


(49c) 


(50a) 


(50b) 
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effective  suppression,  which  we  determine  from  spectra  of 

fluctuating  components  of  the  earth's  magnetic  field  reported 

by  Davidson , 1964,  by  using  Equation  49b  and  presuming  that 

variations  in  direction  of  the  fluctuating  magnetic  field  are 

isotropic.  Light  broken  lines  mark  spectra  expected  with 

2 . . 

noise  power  reduced  by  a factor  (A0)  , as  specified  by  Equa- 

tion 50a,  with  AO  = 0.1  radians. 

From  Figures  12a  and  12b,  we  first  observe  that  spectral 
density  of  estiif.ated  spectra  increases  with  decreasing  fre- 
quency at  about  the  same  rate  as  spectral  density  of  magnetic 

field  fluctuations;  namely,  spectral  densities  are  proportional 
2 fi 

to  (1/f)  (Davidson,  1964).  Next,  we  see  that  magnitudes 

of  observed  spectra,  in  both  cases,  are  about  at  the  level 

2 

expected  when  gradient  fluctuations  are  reduced  by  a factor  (A3) 
with  A0  = 0.1  radians.  Finally,  the  expected  increase  in  level  of 
observed  spectra  for  (p  = tt/2  is  not  evident  in  Figure  12a. 

Equation  50b  then  suggests  the  conclusion  that  the  verti- 
cal component  of  the  fluctuating  magnetic  field  is  suppressed 
at  the  test  site.  The  conclusion  is  consistent  with  data  from 
magnetic  observatories  (Campbell , 1975)  that  show  the  vertical 
component  of  magnetic  field  fluctuations  at  mid  latitudes  is 
appreciably  smaller  than  the  horizontal  component. 

The  slight  increase  in  level  at  high  frequencies  of  the 
observed  spectrum  shown  in  Figure  12b  results  from  a class 
Pc  5 micropulsation  event  (Campbell , 1967;  Sai to , 1969) 
evident  in  the  time  series  record  shown  in  Figure  13. 

Pulsations  begin  at  about  0400  hrs  (PST)  on  6 February  1976 
and  continue  for^SOOO  seconds  as  a damped  sinusoidal  oscil- 
latxon  having  a period  of  roughly  500  seconds.  Maximum  peak- 
to-peak  amplitude  of  the  pulsation  is  about  350  nT. 
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Spectra  of  gradiometer  response  observed  when  the  iron 
sphere  is  placed  3.7  m north  of  the  gradiometer  are  indistin- 
guishable from  spectra  observed  in  a magnetically  quiet  en- 
vironment. Moving  the  sphere  from  1.8  m to  3.7  m reduces 

g 

noise  power  by  a factor  of  (3.7/1. 8)  or  about  300,  and  so 
we  expect  spectral  levels  comparable  to  inherent  noise  of 
the  instrument. 
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Section  VIII 

OPERATION  AT  THE  OCEANOGRAPHIC  TOWER 

Our  field  te.^ts  show  that  the  construct  of  an  equivalent 
dipole  provides  a useful  means  of  describing  response  of  the 
gradiometer  to  ambient  gradients  that  are  sensibly  constant 
over  the  distance  separating  pickup  loops  (25  cm) . Pre- 
liminary tests  of  noise  suppression  techniques,  developed  by 
describing  gradiometer  response  in  terms  of  an  equivalent 
dipole,  demonstrate  reduction  of  noise  power  by  a factor  of 
100  using  crude  means  of  aligning  the  gradiometer  and  suggest 

thau  reduction  of  noise  is  limited  by  alignment  precision 

_ -> 

alone.  By  improving  alignment  to  a precision  of  10  radians 
or  0.05  degrees,  we  expect  to  reduce  noise  power  by  a factor 
of  10^. 


Here,  we  use  the  construct  of  an  equivalent  dipole  to 
estimate  stalxlity  and  alignment  required  to  suppress  noise 
enough  to  afford  certain  measurement  of  fluctuating  gradients 
of  magnetic  fields  generated  above  the  surface  by  oceanic 
internal  waves  passing  the  oceanographic  tower  operated  by 
the  Naval  Undersea  Center.  The  tower  is  located  about  one 
mile  offshore  near  San  Diego,  California  in  water  18  m deep. 
Although  internal  waves  passing  the  tower  provide  a well 
characterized  and  readily  accessible  source  for  first  measure- 
ments of  gradients  of  magnetic  fields  generated  by  internal 
wav(js,  magnetization  currents  in  the  steel  structure  of  the 
tower  produce  steady  as  well  as  fluctuating  local  gradients 
that  make  the  task  of  measuring  fluctuating  gradients  from 
internal  waves  an  exacting  one. 

In  what  follows,  we  first  describe  gradients  of  the  steady 
magnetic  field  near  the  tower,  give  typical  spectra  of  fluc- 
tuating gradients  expected  from  internal  waves  passing  the 
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tower,  and  then  soecify  stability  and  alignment  required  to 
afford  certain  measurement  of  fluctuating  gradients  from 
internal  waves. 

A.  GRADIENTS  OF  THE  STEADY  MAGNETIC  FIELD  NEAR  THE  TOWER 

To  obtain  requisite  information  for  assessing  effects  of 
gradients  resulting  from  magnetization  currents  in  the  tower 
structure,  we  used  a fluxgate  gradiometer  to  measure  the 
steady  gradient  field  in  the  vicinity  of  the  tcwer  (Gillespie 
and  Rodney,  1976) . We  measured  steady  gradients  at  nine  posi- 
tions due  west  of  the  tower  in  a vertical  plane  containing  its 
centerline.  The  pair  of  coordinates  (Z2,  Z^)  specify  a posi- 
tion, where  the  coordinate  Z2  gives  its  horizontal  distance 
in  meters  due  west  of  the  centerline  of  the  tower  and  the 
coordinate  7.^,  its  height  in  meters  above  the  ocean  bottom. 

The  water  depth  is  about  18  m. 

Table  I lists  eigenvalues  (A^,  A2,  A^)  and  polar  angles 
(4)^,  4>2,  02'*  4>3'  ^3)  giving  directions  of  a set  of  eigen- 
vectors e^,  e^)  of  the  steady  gradient  field  at  the  nine 

positions.  Figure  14  defines  the  fiducial  basis  used  to 
reckon  directions  of  eigenvectors.  The  z^  axis  is  directed 
vertically  downward;  the  z^  axis,  northward;  and  the  §2  axis, 
eastward.  Table  I shows  that  the  eigenvector  62  points  nearly 
northward  and  10  to  20  degrees  above  horizontal.  The  plane 

A 

cor.taining  axes  of  equivalent  dipoles,  which  is  normal  to  62 » 
is  then  tilted  10  to  20  degrees  about  a horizontal  axis  pointing 
eastward  and  so  is  close  to  a vertical  plane  of  symmetry  of 
the  tower  structure.  The  structure  comprises  four  steel  pilings 
set  in  a square  array  aligned  north-south  and  east-west. 
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EIGNEVALUES  X^,  X^,  AND  X^  (IN  UNITS  OF  nT/m)  AND 
POLAR  ANGLES  (IN  DEGREES)  OF  THE  EIGENVECTOR  SET  (e, , ^ 
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Table  II  lists  parameters  specifying  location,  orientation, 
and  gradient  strength  of  the  equivalent  dipole  corresponding 
to  the  set  of  eigenvectors  {e^,  e^)  at  each  position.  It 

shows  that  the  position  vector  locating  equivalent  dipoles 
lies  in  the  vertical  plane  of  symmetry  of  the  tower  structure 
((j)  ^ 90°)  . Figure  15  depicts  location  and  orientation  in  the 
plane  of  symmetry  of  the  equivalent  dipole  for  each  position. 
Dashed  lines  indicate  locations  of  dipole  images.  Lengths 
of  vectors  indicating  orientation  of  dipole  axes  are  propor- 
tional to  respective  dipole  moments.  Gradient  strength,  g, 
of  equivalent  dipoles  decreases  with  increasing  horizontal 
distance,  from  the  tower  centerline  in  proportion  to 

with  n = 2.55. 

In  a vertical  plane  of  symmetry  of  the  tower  structure, 
then,  gradients  from  magnetization  currents  in  the  structure 
are  represented  by  gradients  of  a dipole  located  in  the  plane 
of  symmetry.  Asymmetry  of  the  superstructure  of  the  tower, 
however,  tilts  the  axis  of  the  equivalent  dipole  a few  degrees 
out  of  the  plane  of  symmetry. 

B.  SPECTRA  OF  FLUCTUATING  GRADIENTS  EXPECTED  FROM  INTERNAL 

WAVES  PASSING  THE  TOWER 

Magnetic  fields  generated  above  the  surface  by  progressive 

waves  in  a stratified  ocean  are  circularly  polarized  in  a 

vertical  plane  normal  to  wave  crests,  and  their  magnitudes 

decrease  exponentially  with  height,  h,  above  the  surface  as 
""  kh 

e , where  k denotes  wave  number  (Podney , 1975) . As  a result, 
gradients  parallel  to  wa'^e  crests  vanish,  and  the  matrix 
(Podney,  1976) 

G(h,r,  t;k,(i))  = g(l?,(i))e 


sin  (wt-k*r+n)  0 cos  (wt-k*r+ri) 


cos (wt-k.r+n)  0 -sin (wt-k • r+n ) 


(51) 
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represents  gradients  at  a height  h above  the  surface  in  a basis 
having  its  axis  pointing  vertically  downward;  its  axis, 
in  the  direction  of  wave  propagation  k;  and  its  axis,  parallel 
to  wave  crests.  The  scalar  g(k,(ii)  gives  the  strength  of  gra- 
dients generated  at  the  surface  by  a wave  having  a frequency 

^ -A 

tii  and  wave  vector  k,  and  n(k)  is  a phase  shift. 

From  I'quations  51  and  12b,  we  conclude  that  an  equivalent 
dipole  representing  gradients  generated  at  a point  above  the 
surface  by  an  ocean  wave  circles  a sphere  of  unit  radius  about  the 
point  on  a great  circle  path  in  a vertical  plane  normal  to 
wave  crests.  The  axis  of  the  dipole  is  tangent  to  the  sphere 
in  a vertical  plane  normal  to  wave  crests,  x = = it/2  , 

and  circles  the  sphere  at  a frequency  equal  to  one  half  the 
wave  frequency;  namely,  29 (t)  = wt  - k.r  + n.  Its  moment  is 
proportional  to  wave  amplitude. 


For  surface  waves,  the  relation 


g(k,co)  = ^(k,o)) 


l+f)  ^sin^(})^+(l-f)  ^cos^(})^cos^6^  (52a) 


gives  the  strength  of  gradients  generated  at  the  surface  by  a 
surface  wave  that  has  a wave  height  C(k,co)  and  heads  at  an 
angle  9^^  east  of  magnetic  north.  Here,  B denotes  magnitude  of 
the  earth's  magnetic  field;  cf)^,  its  dip  angle;  a,  the  electrical 
conductivity  of  seawater  ('^4  mhos/m)  ; and 


sinh  kD 


where  D denotes  water  depth.  The  relation 


(52b) 


cot  n = 


■ 


cot(f)^cos9^ 


(52c) 
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gives  the  phase  shift. 


For  internal  waves,  the  relation 


(IJ  ao)\  „ 22 

^)B(sin  (})^+cos  ({)^cos 


(d,k) 


gives  the  strength  of  gradients  generated  at  the  surface  by 
an  nth  mode  internal  wave  that  has  a displacement  ^^(d,k,03) 
at  a depth  d(0-^d<D)  and  heads  at  an  angle  6^^  ^ east  of  magne- 
tic north.  Here, 

D 

I (d,k)  = ,k)e"’^^kd?  , 

^ o 


where  <5>^(C,k)  is  an  eigenfunction  giving  the  profile  of  an 
nth  mode  wave.  Eigenfunctions  vanish  at  the  surface  and 
bottom  of  an  ocean  and  are  orthonormal  with  respect  to  the 
weight  function  N^(z);  namely;  <J>^(0,k)  = $^(D,k)  = 0,  and 

•»„(?, k)*^(C,k)N2(C)dc  = 6^^  , 

o 


where  N(z)  is  the  Brunt-Vaisala  frequency  profile  of  a thermo- 
cline.  The  relation 

tan  = -cot(t)^cos0^^^ 


(53a) 


(53b) 


(53c) 


(53d) 


gives  the  phase  shift. 


We  use  Equation  53a  to  estimate  spectra  of  fluctuating 
gradients  expected  from  it.ternal  waves,  it  shows 
that  strength  of  gradients  generated  at  the  surface  is  pro- 
portional to  wave  displacement  at  depth  and  that  the  coef- 
ficient depends  on  a weighted  integral  of  an  eigenfunction 
over  depth,  as  expressed  by  Equation  53b.  From  spectral 
measurements  of  isotherm  displacements  at  the  top  of  the 
thermocline  together  with  the  measured  thermocline  profile, 
then,  we  compute  spectra  of  fluctuating  gradients  expected 
from  internal  waves  passing  the  tower.*  The  waves  pass 
the  tower  in  packets  of  long-crested  first-mode  waves 
running  close  to  due  east,  so  9^^  ^ tt/2.  Magnitude  of  the 
earth's  magnetic  field  at  the  tower  is  about  47,200  nT 
and  its  dip  angle  is  about  61°  (Gillespie  and  Podney,  1976) 


Figure  1 shows  three  typical  spectra  of  transverse  gra- 
dients  expected  7 m above  the  surface  from  internal  waves  passing 
the  tower  comoared  to  the  Inherent  noice  spectrum  of  the  instru- 
ment. spectra  shown  in  Figure  1 give  the  expected  spectral  re- 
sponse when  pickup  loops  of  the  gradiometer  face  in  the  direction 

of  wave  propagation  or,  more  specifically,  when  the  Xj  axis  o 

I-  ■ i <=  n^T-allel  to  wave  crests.  As  is  evident, 
the  gradiometer  basis  is  parallel  to  wave  ^ ^ ^ ^ 

spectral  densities  of  fluctuating  gradients  expected  from  inter 
nal  waves  are  well  above  the  level  of  instrument  noise. 


* Fluctuating  gradients  from  surface  waves  are  neglibly 
small  at  frequencies  of  internal  waves  at  the  tower. 
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C.  ALIGNMENT  AND  STABILITY  REQUIREMENTS 


We  plan  to  jut  the  gradiometer  over  water  on  a rigid  nun- 
magnetic  cantilever  extending  horizontally  25  meters  from  the 
centerline  of  the  tower  in  a vertical  plane  of  symmetry.  The 
gradient  field  of  a dipole  located  in  the  plane  of  symmetry 
and  having  a gradient  strength,  g,  of  about  200  nT/m  or  2 X 

5 

10  pT/m  gives  gradients  at  the  position  of  the  gradiometer 
from  magnetization  currents  in  the  tower  structure.  We  re- 

4 

present  fluctuations  in  response  of  the  gradiometer,  which 
result  both  from  fluctuations  in  its  position  and  orientation 
and  from  fluctuations  in  magnetization  currents  in  the  tower, 
in  terms  of  fluctuations  in  location,  orientation,  and  gra- 
dient strengtn  of  the  dipole. 

1 . Alignment 

Without  aligning  the  gradiometer  to  suppress  noise  from 
fluctuating  magnetization  currents  in  the  tower,  we  expect 
a spectral  density  of  fluctuating  gradients,  N^,  given  by 
the  relation 

Ni(f)  = g^[b(f)/e]^  = 1.25  X 10"'^(l/f)^'®(pT/m)^/-Hz  , 


where  the  relation 

b^(f)  = 6.31  X 10"®(l/f)2*®  (nT)^/Hz 

gives  the  spectral  density  of  magnetic  field  f luctuatiOiis 
ov/ing  to  changing  ionospheric  currents  as  reported  by  David- 
•lon , 1964.  We  see  from  Figure  1 and  Equation  54a,  then/ 
that  gradient  fluctuations  from  fluctuating  magnetization 
currents  in  the  tower  mask  gradient  fluctuations  expected 
from  internal  waves  unless  the  gradiometer  is  aligned  to 
suppress  gradient  fluctuations  from  magnetization  currents  in 
the  tower. 


(54a) 


(54b) 
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By  aligning  the  gradiometer  for  a null  response  with  a 
precision  of  0.01  radians  or  0.6  degrees,  liowever,  we  re- 
duce the  spectral  density  of  gradient  fluctuations  owing  to 
fluctuating  magnetization  currents  in  the  tower  by  a factor 
of  10~^,  which  then  puts  it  below  the  level  of  instrument 
noise.  We  align  for  a null  response  by  orienting  the  gra- 
diometer so  that  the  dipole  representing  gradients  of  magne- 
tization currents  in  the  tower  is  located  along  the  axis 
of  the  gradiometer  basis  with  its  axis  in  the  plane  containing 
pickup  loops.  The  plane  containing  pickup  loops  is  then  close 
to  a vertical  plane  of  symmetry  of  the  tower  structure. 

In  order  to  align  the  gradiometer  for  a null  response  and 
at  the  same  time  to  face  pickup  loops  in  the  direction  of  in- 
ternal wave  propagation,  we  choose  the  plane  of  symmetry  so 
that  internal  waves  cross  the  plane  at  nearly  normal  inci- 
dence. Namely,  we  extend  the  cantilever  from  either  the 
north  or  the  south  side  of  the  tower.  Aligning  the  gradio— 
meter  for  a null  response  to  suppress  gradient  fluctuations 
from  fluctuating  magnetization  currents  in  the  tower  then 
gives  a near  maximum  response  to  gradient  fluctuations  from 
internal  waves  as  well. 


i 

I 


When  the  gradiometer  is  aligned  for  a null  response,  it 
responds  only  to  displacements  and  changes  in  orientation  that 
produce  an  effective  rotation  about  the  gradiometer  axis  of 
the  equivalent  dipole  representing  gradients  from  magnetiza- 
tion currents  in  the  tower.  Namely,  the  relation 


g23(t)  = 6(j)(t)  g sinx 


(55) 
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then  gives  gradient  fluctuations  resulting  from  irregular  motions, 

where  6(J)(t)  is  the  effective  angle  of  rotation  about  the  gradio- 

meter  axis,  g is  the  gradient  strength  of  the  equivalent  dipole, 

and  the  angle  x qives  inclination  of  its  axis  to  the  position 

vector  locating  the  dipole  in  the  gradiometer  basis.  Our  measure- 

5 o 

ments  give  g = 2 X 10  pT/m  and  x = . From  Equation  55a, 

we  then  expect  a spectral  density  of  fluctuating  gradients  re- 
sulting from  irregular  motion,  N^,  given  by  the  relation 

Ng(f)  = 0.77$(f)  (pT/m)^/Hz  , (55b) 

where  <I>(f)  is  a spectral  density  of  fluctuations  in  the  effective 

angle  of  rotation  about  the  gradiometer  axis,  in  units  of  (seconds 
2 

of  arc)  /Hz. 

To  keep  spectral  density  of  gradient  fluctuations  from  irre- 
gular motions  at  the  level  of  instrument  noise  in  the  bandwidth 
of  gradient  fluctuations  expected  from  internal  waves,  nominally 
2 X 10*^  Hz  to  5 X 10~^  Hz,  the  cantilever  must  limit  angular 

deviations  of  effective  rotations  about  the  gradiometer  axis  to 

-5 

a few  seconds  of  arc  or  about  10  radians  in  the  bandwidth.  We 

require,  then,  that  structural  design  of  the  cantilever  limit 

rotational  and  translational  fluctuations  of  the  instrument 

“5 

relative  to  the  tower  to  10  radians  ( 2 seconds  of  arc)  and 
0.1  mm  in  the  bandwidth  of  gradient  fluctuations  expected  from 
internal  waves.  A displacement  of  0.1  mm  at  a radius  of  25  m 
gives  an  angular  deviation  of  about  1 second  of  arc. 

3 • Cantilever  Design 

Figure  16  shows  an  elevation  and  top  view  of  the  structural 
desian  developed  by  Mechanics  Research,  Inc.  (Haire  and  Van  Lerberg, 
1976)  to  give  the  requisite  stability  under  wave,  wind,  and  thermal 
loads  expected  during  operations  at  the  tower.  The  cantilever  is 
a four-sided  truss  space  frame  extending  about  20  meters  from  the 
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face  of  the  tower.  An  8 m base  section  made  of  tub  alar  aluminum 
attaches  to  the  tower  supports,  and  a 12  m outer  section  made 
of  a fiberglass  laminate  supports  the  instrument  at  the  tip  of 
the  cantilever.  The  instrument  Dewar  fits  in  a fiberglass  gim- 
bal  mount  wnose  three  independent  axes  of  rotation  provide  means 
of  orienting  the  gradiometer  axis  as  necessary.  A catwalk  running 
the  length  of  the  cantilever  provides  access  to  the  instrument. 
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Section  IX 
CONCLUSION 


At  present,  we  describe  the  response  of  a superconducting 
magnetic  gradiometer  in  terms  of  the  magnetic  field  and  its 
gradients  at  the  midpoint  of  the  gradiometer  axis.  The  descrip- 
tion is  a useful  approximation  provided  gradients  of  the  ambient 
magnetic  field  are  sensibly  constant  over  the  length  of  the 
gradiometer  axis.  The  part  of  the  response  proportional  to 
the  magnetic  field  at  the  midpoint  results  from  slight  dif- 
ferences in  area  and  orientation  of  the  loops  forming  the  pick- 
up circuit  of  the  gradiometer  and  vanishes  for  a perfectly 
balanced  gradiometer.  The  part  proportional  to  a gradient  of 
the  magnetic  field  at  the  midpoint  depends  of  five  independent 
elements  of  a matrix  representing  gradients  of  the  ambient  mag- 
netic field  at  the  midpoint. 

To  concisely  describe  the  gradient  response,  we  use  the 
five  elements  to  define  location,  orientation,  and  moment  of 
an  equivalent  magnetic  dipole  that  gives  gradients  at  the  mid- 
point equal  to  gradients  of  the  ambient  field.  Constructing 
an  equivalent  dipole  lets  us  visualize  the  response  to  ambient 
gradients.  For  example,  the  gradient  field  cf  a dipole  fluc- 
tuating about  a mean  location,  orientation,  and  moment  repre- 
sents fluctuating  gradients  from  nearby  magnetic  objects.  A 
dipole  circling  a sphere  of  unit  radius  about  the  midpoint  on 
a great  circle  path  in  a vertical  plane  norm.al  to  wave  crests 
represents  gradients  from  an  ocean  wave.  A swarm  at'  d.jooli  s 
circling  the  sphere  represents  gradients  from  u f.xc  IJ.  of  ran- 
dom waves. 

Results  of  preliminary  tests  of  gradioa'ot or  response  axe 
consistent  with  the  description  based  on  a fin  t order  ap; ■ xi- 
mation.  Nonetheless,  the  results  are  not  precxre  o..ough  t 
delimit  the  range  of  validity  of  the  first  order  descrjotion. 

We  plan  further  tests  using  more  accurate  gimbals  and  a] ignment 
procedures  both  to  determine  a more  precise  description  of 
gradiometer  response  and  to  mark  limits  of  valid; of  the  first 
order  description. 
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FIGURE  CAPTIONS 


Figure  1 Three  typical  spectra  of  fluctuating  gradients  of 
magnetic  fields  expected  7 m above  the  surface  from  internal  waves 
passing  the  oceanographic  tower  operated  by  the  Naval  Under- 
sea Center  as  compared  to  the  inherent  noise  spectrum  of  the 
instrument.  Spectra  are  estimates  based  on  spectral  measure- 
ments of  isotherm  displacements  made  during  ;'ugust  1972. 

Serial  numbers  marking  spectra  tell  the  month,  day,  and  local 
time  at  the  beginning  of  corresponding  time  series;  for  example, 
08211540  means  Aucrust  21  at  15:40  hours  (PST)  . 


Figure  2a  Illustration  of  the  instrument  showing  Dewar  con- 
struction and  principal  components  of  the  gradiometer  probe. 
A wrapping  of  alternating  layer's  of  fiber  glass  cloth  and 
aluminized  Mylar,  which  is  not  shown,  insulates  the  interior 
vessel.  The  Dewar  stands  48"  high,  is  26"  in  diameter,  and 
weighs  about  250  lbs. 


Figure  2b  Diagram  of  the  superconducting  circuit  of  the  gra- 
diometer. Two  coplanar  pickup  loops,  having  a combined  stlf 
inductance  L^,  are  oppositely  wound  in  series  and  connected 
to  the  primary  winding  (having  a self  inductance  of  a 

transformer  containing  a normal  metal  shield  between  its 
superconducting  windings  that  attenuates  radio  frequency 
interference  (rfi) . The  secondary  winding  of  the  transformer 
has  a self  inductance  L2  and  is  connected  to  a field  coil, 
having  a self  inductance  L^,  that  couples  flux  to  a super- 
conducting quantum  interference  device  (SQUID) . The  SQUID 
has  a self  inductance  L and  is  driven  by  a tank  circuit 
connected  to  room  temperature  electronics.  The  rfi  trans- 
former, SQUID,  and  tank  circuit  are  encapsulated  in  a super- 
conducting shield.  A super  current  I^  flows  in  the  pickup 
circuit,  and  a super  current  I2/  in  the  coupling  circuit. 
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Figure  3 Geometry  of  pickup  loops  comprising  an  axis  of  a 
superconducting  magnetic  gradiometer.  Centers  of  the  two 
pickup  loops  are  separated  by  a distance  2s,  and  their  areas 
are  and  A2.  Normals  to  the  plane  of  each  loop  are  directed 
along  unit  vectors  and  n2 , and  radial  vectors  and  ^2 
;.ocate  points  within  a loop  contour.  The  gradiometer  basis 
is  fixed  midway  between  centers  of  the  loops  and  is 
located  with  respect  to  a fiducial  basis  by  the  vector  R. 


Figure  4 Vector  diagram  depicting  location  and  orientation 
of  a magnetic  dipole,  nT,  with  respect  to  the  gradiometer 
basis  • Polar  angles  4>  and  9 specify  direction  of  the 

position  vector  IT  locating  the  dipole.  The  polar  angle  x 
and  azimuthial  angle  Q specify  direction  of  the  dipole  axis 
with  respect  to  the  position  vector.  The  y^  axis  of  the 
basis  points  along  the  position  vector;  the  ^2  axis, 

along  r X m;  and  the  y^  axis,  along  X m)  X r. 

Figure  5 Record  of  gradiometer  response  during  a clockwise 
rotation  followed  by  a counter  clockwise  rotation  about  its 
axi^.  Rotations  are  made  stepwise  in  sixteen  increments  of 
22.5°,  pausing  at  positions  numbered  1 through  16  that  cor- 
respond to  numbered  positions  in  Figure  6.  Spikes  mark  move- 
ments between  incremental  positions,  and  flats  mark  values 
recorded  while  pausing. 

Figure  6 Orientation  of  positions  1 through  16  with  respect 
to  magnetic  north.  Rotations  begin  with  the  X2  axis  of  the 
gradiometer  basis  pointing  northward  and  proceed  through  posi- 
tions 1,  2,  3,  ...,  16,  1 during  a clockwise  rotation  and 
through  positions  1,  16,  15,  ...,  2,  1 during  a counter  clock- 
wise rotation. 
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Figure  7 Signature  of  hysteresis  observed  during  rotation  at 

"5  -1 

an  imbalance  somewhat  less  than  10  m . The  trace  records 
gradiometer  response  during  a stepwise  oscillation  of  the  X2 
axis  of  the  gradiometer  basis  about  position  5 marked  in 
Figure  6,  beginning  with  a clockwise  movement  to  position  7. 
Numbered  flats  mark  values  recorded  while  dwelling  at  cor- 
responding positions.  Spikes  mark  movements  between  positions 
in  the  sense  denoted  by  vertical  arrows. 

Figure  8 Data  depicting  the  noise  spectrum  of  the  gradiometer 

operating  in  a magnetically  quiet  environment.  We  use  a 

calibration  of  4.5  mV/(pT/m)  to  obtain  spectral  density  in 
2 

units  of  (pT/m)  /Hz,  marked  on  the  left  hand  coordinate  scale, 
and  a calibration  of  12.7  V/(l)^  to  obtain  spectral  density  in 
units  of  marked  on  the  right  hand  coordinate  scale. 

Points  marked  by  a circular  dot  • represent  data  obtained  from 
digital  records.  Points  marked  by  a triangular  dot  ▼repre- 
sent data  obtained  from  a real-time  spectrum  analyzer.  Error 
bars  mark  limits  corresponding  to  two  standard  deviations. 

Figure  9 Data  depicting  the  spectrum  of  intrinsic  noise  of 
the  SQUID  sensor  compared  to  the  noise  spectrum  of  the  gradio- 
meter operating  in  a magnetically  quiet  environment.  The  left 
hand  coordinate  scale  gives  spectral  density  in  units  of  (pT/ 

m)  /Hz,  and  the  right  hand  scale,  in  units  of  <p  /Hz.  The 

-4  ° ■ 

conversion  factor  is  3.5  X 10  (f>^/  (pT/m)  . The  solid  line  de- 

lineates the  spectrum  of  inherent  gradiometer  noise  shown  in 
Figure  8,  as  represented  by  the  spectral  form  S(f)  = (2.1  X 
10"^)/f^'^(pT/m) ^/Hz,  for  f‘0.1  Hz,  and  S(f)  = 0.03(pT/m)^/ 

Hz,  for  f>0.1  Hz. 
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Figure  10  Expected  enhancement  of  instrument  performance  from 

use  of  a thin-film  dc  type  SQUID  sensor.  Data  points  depict 

the  spectrum  of  intrinsic  noise  of  a thin-film  dc  type  SQUID 

in  units  of  /Hz  marked  on  the  right  hand  coordinate  scale, 

and  the  heavy  dashed  line  marks  a corresponding  expected 

2 

spectrum  of  gradiometer  noise  in  units  of  (pT/m)  /Hz  marked 

on  the  left  hand  coordinate  scale.  The  light  dashed  curve 

delineates  the  spectrum  of  intrinsic  noise  of  a torodial 

point-contact  type  sensor  as  used  in  the  instrument,  in  units 
2 

of  $ /Hz,  and  the  solid  line  marks  the  corresponding  spectrum 

° 2 
of  inherent  gradiometer  noise  in  units  of  (pT/m)  /Hz. 

Figure  11  Magnetic  dipole  describing  the  induced  gradient 
field  at  a radial  position  r"  outside  an  iron  sphere  magnetized 
by  the  earth's  magnetic  field.  Eigenvectors  and  e^  mark 
principal  axes  in  a meridian  plane  for  several  values  of  the 
anale  x = ^ngle,  0^,  of  the  earth's 

magnetic  field  at  the  site  is  about  55°,  and  the  dip  angle 
e,,  of  the  position  vector,  ?,  is  +5°  or  so.  The  angle  a 

^ A O 

between  eigenvector  e^  and  the  position  vector  is  about  68 
for  = 0 *^d  ~ ’ 

Figure  12a  Spectrum  of  gradient  fluctuations  from  an  iron 
sphere  placed  1.8  m north  of  the  gradiometer  with  pickup 
loops  facing  nominally  north-south  so  that  4>  = -n/2.  The 
corresponding  time  series,  called  POSTA  9,  begins  4 February 
1976  a"-  16:22  hours  (PST)  and  ends  5 February  1976  at  09:40 
hours  (PST) . The  heavy  solid  line  delineates  the  spectrum  of 
inherent  instrument  noise  shown  in  Figu^'e  8,  and  light  and 
heavy  dashed  lines  mark  spectra  expected  with  and  x-’ituout 
suppression  of  noise,  respectively,  owing  to  aligning  for  a 
null  response. 
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Figure  12b  Spectrum  of  aradient  fluctuations  from  an  iron 
sphere  placed  1.8  m north  of  the  gradiometer  with  pickup  loops 
facing  nominally  east-west  so  that  (|)  = 0.  The  corresponding 
time  series,  called  POSTA  10,  begins  5 February  1976  at  15:15 
hours  (PST)  and  ends  6 February  1976  at  16:37  hours  (PST). 

The  heavy  solid  line  delineates  the  spectrum  of  inherent  instru- 
ment noise  shown  in  Figure  8,  and  light  and  heavy  dashed  lines 
mark  spectra  expected  with  and  without  suppression  of  noise, 
respectively,  owing  to  aligning  for  a null  response. 

Figure  13  Section  of  time  series  record  POSTA  10  showing  a 
class  P::  5 micropulsation  event  recorded  on  6 February  1976. 
Pulsations  begin  at  about  0400  hours  (PST)  and  continue  for 
^•5  "'0  seconds  as  a damped  sinusoidal  oscillation  having  a 
period  of  roughly  500  seconds.  (Gradient  magnitudes  marked 
on  the  left  hand  coordinate  scale  are  relative  to  an  arbi- 
trary zero.) 

Figure  14  Fiducial  basis  used  to  define  measurements  of 

gradients  in  the  vicinity  of  the  oceanographic  tower.  The 
z^  axis  points  vertically  downward;  the  axis,  northv/ard; 
and  the  Z2  axis,  eastward.  Polar  angles  (<})j^,  0^^;  '^2’  ®2'  *3' 

0^)  specifying  directions  of  the  set  of  eigenvectors  ^2' 

at  each  position  are  listed  in  Table  I.  Eigenvector  ^2 
points  nearly  northward  and  10°  to  20°  above  horizontal,  and 
eigenvector  e^  points  nearly  due  west  and  5 to  15  below 
horizontal . 

Figure  15  Location  and  orientation  in  the  vertical  plane  of 
symmetry  of  an  equivalent  dipole  for  gradients  at  each  position 
near  the  oceanographic  tower.  Heavy  dashed  lines  mark  loca- 
tions of  dipole  images.  Lengths  of  vectors  indicating  orien- 
tation of  dipoles  are  proportional  to  respective  dipole 
moments  and  decrease  with  increasing  horizontal  distance. 
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from  the  tower  centerline  in  proportion  to  with  n = 2.55. 

Coordinates  Z^)  specify  a position:  gives  its  horizontal 

distance  in  meters  due  west  of  the  centerline  of  the  tower  and 
gives  its  height  in  meters  above  the  ocean  bottom.  Positions 
are  nominally  spaced  at  2 m intervals  vertically  and  3 m inter- 
vals horizontally. 

Figure  16  Elevation  and  top  view  of  the  cantxlever  support  for 
the  gradiometer. 
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APPENDIX 


GRADIENTS  OF  MAGNETIC  FIELDS  IN  FREE  SPACE 


We  represent  gradients  of  a magnetic  field,  b(R,t), 
by  a matrix  G having  elements  g^^ (R,t) , where 

■ “^ij  x^.V(Xj-b),  (i,j  = 1,2,3)  (la) 


in  an  orthogonal  basis  | x^  | . An  element  g^^  represents 
the  gradient  in  a direction  x.  of  the  component  of  magnetic 

A 1 


field  in  a direction  x. 

1 


The  relation 


r(u,0)  = ^ 

ifi 

then,  gives  the  gradient  in  a direction  u of  the  component 
of  magnetic  field  in  a direction  v,  where  and  v^  are 
components  of  the  unit  vectors  u and  v in  the  basis | x^  | . 
We  express  I’(u,v)  as  a matrix  product  by  using  the  relation 


r(u,v)  = UGV 


(Ic) 


where  U and  V are  column  matrices  representing  unit  vectors 
u and  V,  respectively,  and  tilde  marks  a transposed  matrix. 

Equation  la  shows  that  the  sum  of  diagonal  elements  of 
a gradient  matrix  vanishes  because  the  divergence  of  a 
magnetic  field  vanishes.  Moreover,  the  vector  identity 

g.  . - g . . = (x.  X X . ) • (V  X b)  (2) 

tells  us  that  matrices  representing  gradients  of  an  irrota- 
tional  magnetic  field  are  symmetric;  namely,  G = G. 
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A magnetic  field  in  free  space  is  both  nondiv/ergent  and 
irrotational  and  so  is  the  gradient  of  a harmonic  function, 
iji;  nameiv,  b = where  = 0.  Conseauently , gradients 

of  a magnetic  field  in  free  space  are  represented  by  a 
symmetric  matrix  that  has  a vanishing  trace*,  'vhose  elements 
are  given  by 


' ax.  ax  . 

1 3 

Only  five  of  the  nine  elements  of  a matrix  representing 
gradients  of  a magnetic  field  in  free  space,  then,  are 
independent . 

An  axis  of  a perfect  gradiometer  that  is  aligned  in  a 
direction  x^  and  senses  the  component  of  magnetic  field  in 
a direction  x^  measures  the  element  g^j  of  a gradient  matrix. 

A gradiometer  comprising  two  orthogonal  axes  that  measure 
longitudinal  elements  g^^^  and  together  with  three  mutu- 

ally orthogonal  axes  that  measure  transverse  elements  9^2, 
g^2'  ^23'  determines  gradients  of  a magnetic  field 

in  free  space  in  the  orthogonal  basis  defined  by  its  axes. 

Because  of  Equation  lb,  matrix  elements  measr.red  in  the 
qradiometer  basis  in  turn  determine  gradients  in  everv 
direction  (a  of  the  component  of  maqnetic  field  in  arv  direction  v. 

A.  GRADIOMETER  ROTATIONS 

As  a gradiometer  rotates,  matrix  elements  measured  in 
the  gradiometer  basis  change  continuously.  We  describe 

gradiometer  rotations  by  a matrix  R that  reoresents  a rotation 

/ A ' • ( ^ J 

from  a basis  j ^ ^ basis  ^ 


* The  trace  of  a matrix  is  the  sum  of  its  diagonal  elements. 
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The  matrix  R is  orthoaonal * , and  its  elements  are  aiven  by 


r . . 
1] 


A /' 

X . • y . 
1 


(4) 


For  a vector  represented  by  a column  matrix  in  a 
basis  I , the  relation 

U = RU 

y X 

qivps  its  representation  U in  a basis  1 / where  R rep- 

y . i A 1 ' ‘ / A ) 

resents  rotation  from  the  basis  < x^  ^ to  the  basis  < y^^  > • 

Because  the  matrix  product  expressed  by  Eauation  Ic  is  in- 
variant, we  find  that  the  relation 


G = RG  R 

y X 

specifies  the  matrix  representing  gradients  of  a magnetic 
field  in  a basis  in  terms  of  their  representation  G^ 

in  a basis  |x^|  . Elements  of  a gradient  matrix  measured 
in  one  basis  then  determine  the  gradient  matrix  in  every 
other  basis  obtained  by  rotating  the  gradiometer. 


We  express  a matrix  representing  a rotation  as  the  pro- 
duct of  three  matrices  representing  three  successive,  inde- 
pendent rotations  through  the  Eulerian  angles  4’,  0,  and  f 
defined  in  Figure  Al  by  the  product  (Indritz,  1963) 


* The  relation  RR  = RR  = I , where  I denotes  the  unit  matrix, 
defines  an  orthogonal  matrix. 
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R((P,  G,  t) 


Z(1))  X(G)  Z(>f) 


(6a) 


where 


Z($) 


cos 

sin 


<I) 


0 


-sin  0 

cos  <1>  0 

0 1 


(6b) 


and  represents  a rotation  through  an  angle  about  a z-axis, 
and 


X(0) 


10  0 
0 cos  0 -sin  0 

0 sin  0 cos  0 


(6c) 


and  represents  a rotation  through  an  angle  0 about  an  x-axis. 
We  note  that  Z(<I>)  = Z(-'I')  and  )T(0)  = X(-0),  so  that  R(<I>,  0, 

H*)  = R(-4' , -0,  -<I>)  . 


A rotation  specified  by  Euler  angles  $ , 0»  and  4',  however, 
is  equivalent  to  rotation  through  an  angle  4)  about  some  axis 
specified  by  polar  angles  and  0^  with  respect  to  the  fi- 
ducial basis  (Goldstein,  1959).  Namely, 

Z(4>)  = P(4)q,  0q)R(4>,  0,  H')P(4)^,0j^)  , 


so  that 

R(4)o,  e^;  4.)  = P(4>q,  9^)Z(4))P(4>o,  e^)  = 0,  H-)  , (7! 

where  P(4)  , 0^)  = Z (4>q+tt/2)  X(0^)  and  represents  a rotation 
from  the  fiducial  basis  to  a basis  having  its  z-axis  aligned 
with  the  axis  of  rotation.  From  Equation  7b  we  find  that 
angles  4>  , 0^,  and  4>  are  expressed  in  terms  of  Euierian  angles 
specifying  the  same  rotation  by  the  relations 
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'1' 


(8a) 


o 


(j)  - 


2 


tan  0^ 
o 


tan  0/2 

4 + S' 


sin  I 


(8b) 


and 


cos  (f/2  = cos  0/2  cos 


<t>  + S' 


(8c) 


To  ma)ce  the  dependence  on  the  angle  of  rotation  ())  explicit, 
we  write  a rotation  matrix  as  the  sum 

R(4>o»  Oq'  " ^l^*^o'  ®o^  ^ Q((()Q,e^)sin^  (9a) 

where 

Ei«>o'  V = '■'♦o'  «o>  ' 

(9c) 


E2(4>o.  ®o^  ' 


and 


Q(*„,  9„)  = P(4'o'  “o' 


(9d) 


with 


0 

0 

0 

0 

0 

1 

0 

-1 

0 

= 

0 

0 

-1 

e„  = 

0 

0 

0 

, and  - 

1 

0 

0 

0 

1 

0 

y 

-1 

0 

0 

0 

0 

0 

The  matrix  Q corresponds  to  a sum  of  differential  rotations 
about  x,y,  and  z axes  of  the  fiducial  basis,  as  expressed  by 
the  expansion 
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i 


! 

; 

. 


and  for  an  infinitesimal  rotation  Arp, 


Rlip  , 0 ;A(}))  = I + Q(<ti_/  0^)A(j) 

00  00 

Having  presented  means  for  describing  gradiometer  rota 
tions,  vre  now  investigate  changes  in  gradient  signals  pro- 
duced by  changes  in  orientation  of  a gradiometer.  We  first 
consider  rotation  about  an  axis  of  the  gradiometer  basis 
and  then  about  an  arbitrary  axis. 


1 . Rotation  about  a gradiometer  axis 

Naming  of  axes  of  a gradiometer  basis  is  arbitrary,  and 
so  we  call  the  axis  of  rotation  the  z— axis.  Following 
rotation  of  a gradiometer  through  an  angle  <p  about  its  2- 
axis,  the  gradient  matrix  measured  in  the  gradiometer  basis 
r((|))  , is  expressed  in  terms  cf  its  initial  value  G by  the 
relation 


r(4))  = Z((b)GZ(4))  . 


Consequently,  its  elements  are  given  in  terms  of  their 

initial  values  g.  , by  the  relations 


^11 

^22 

^33 

'^12 


-g,,/2  + I <511+933/21^ 


'33' 
-933/2  - 


1/2 


cos2  (4)]^ 


^12  ^^ll’*'^33'^^^ 


1/2 


cos2  {<)), 


.33 


^12  ^^ll'*'^33/^^^ 


1/2 


sin2(4)j^  - <(i) 


<t>) 

<t>) 
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(10b) 


(11a) 

(12a) 

(12b) 

(12c) 

(12d) 


I .-IJII 


^13  = ^^13  923^ 

and 

^23=  OL  + 

where 


(12e) 


(12f) 


tan  2<{)^  = ^ *^33^^^ 


(12g) 


and 


tan  <p  2 ~ ^13'^*^23 


(12h) 


We  observe  that  the  longitudinal  gradient  Y33  is  con- 
stant, that  longitudinal  gradients  Y33  and  Y22  sinusoi- 
dally about  the  mean  value  {g^j^  + ^22^/^'  transverse 

gradients  vary  sinusoidally  about  zero.  As  a result,  we  can 
always  find  a null  position  for  a transverse  gradient  by 
rotating  a gradiometer  about  any  one  of  its  axes,  since 
specification  of  a z-axis  is  arbitrary.  We  can  not^find 
a null  position  for  longitudinal  gradients  unless  9^2^ 

2 „ r>T-  a . Because  the  sum  of 

9ii522'  923  >922933'  °"  9l3  >9ll933‘ 

diagonal  elements  vanishes,  however,  at  least  two  of  the 
conditions  are  satisfied.* 


2 . Rotation  about  an  arbitrary  axis 

During  rotation  of  a gradiometer  about  an  axis  other  than 
one  of  its  three  mutually  perpendicular  axes,  each  element 

* Nonetheless,  determination  of  a null  position  for 
tudinal  qradients  is  impractical  because  mean  values  of 
sinusoidal  oscillations  are  experimentally  indeterminant . 
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of  the  gradient  matrix  oscillates  abov t a nonzero  mean  value. 
Oscillations  are  a sum  of  sinusoids  having  arguments  c|)  and 
2c|),  where  <p  is  the  angle  of  rotation.  Amplitudes  of  sinusoi- 
dal oscillations  depend  on  initial  orientation  of  the  gradio- 
meter  and  polar  angles  and  9^  specifying  the  axis  of 
rotation  with  respect  to  its  initial  orientation. 

The  relation 

r(c()o,  0^;  (}))  = R(4)q.  4))GR((j)^,  6^;  t|):)  (13a) 

gives  the  gradient  matrix  obtained  by  rotating  a gradio- 

meter  through  an  angle  cf  about  an  axis  specified  by  polar 

angles  d)  and  0 with  respect  to  its  initial  oricn+'ation  in 
^ o o 

terms  of  the  initial  value  G.  When  the  rotation  matrix 

R((b  , 6 ; (b)  is  expressed  explicitly  in  terms  of  the  rota- 
o o 

tion  angle  by  Equation  9a,  we  find  that  9^/  <P) 

is  expressed  as  a sum  of  sinusoids  with  arguments  d)  and  2(J) 
by  the  relation 

r(d)Q,  9^;  (b)  = E^GE^  + 1/2(E2GE2  - QGQ)  + 1/2  (E2GS2  + 

QGC)cos  2(b  + 1/2(E2GQ  - QGE2)sin  2(p  + (13b) 

(Ej^GE2  + E2GE^)coS(b  + (E^GQ  - QGE^);:in(b  , 

where  matrices  ' ^2^‘^o'  ®o^'  ^^*^0'  ®o^  given 

by  Equations  9b,  9c,  and  9d. 

Because  mean  values  and  amplitudes  of  sinusoidal  com- 
ponents depend  on  polar  angles  <b^  and  0^  and  on  the  initial 
matrix  G,  we  can  not  always  find  a null  position  for  trans- 
verse gradients  by  rotating  a gradiometer  about  an  axis 
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other  than  one  of  its  three  mutually  perpendicular  axes. 
Nonetheless,  for  axes  of  rotation  lying  within  a narrow 
cone  about  a gradiometer  axis,  defined  by  «o«  ^ " 

z-axis,  we  find  from  Equation  9a  that 


‘’o' 


♦ ) S'  Z(4.) 


I + 6^  E(t^, 


6o«  1 . 


where 


e(())q,  4>)  = 


sin  +sin<|)^ 


"x 


cos  {(J)-(t)Q)  -COS({)^ 


e / 

y 


so  that 


r (4>o/  QqJ  4>)  = ^ 


r {(j))E{(j)^,  (t))-E{(})^,  <j))r{(j)) 


where  1(4))  is  given  by  Equation  11a.  We  conclude  from  Equa- 
tion 14c  that  transverse  gradients  vanish  during  rotation 
of  a gradiometer  about  an  axis  lying  within  a narrow  cone 
about  any  one  of  its  three  mutually  perpendicular  axes. 


3.  Infinitesimal  rotations 


For  an  infinitesimal  rotation  through  an  angle  A4>  about 
an  axis  specified  bv  polar  angles  4>q  and  9^,  we  find  from 
Equations  13a  and  10b  that  the  differential  of  the  gradient 
matrix  is  given  by 


1(4) 


o 


o 


A4>)  - G = A4i 


00(4)^,  9^)  - d^)G 


where  the  matrix  Q(4>q.  0^)  is  expressed  as  a sum  of  differ 
ential  rotations  about  x,y  and  z axes  of  the  fiducial  basi 
by  Equation  10a.  From  Equations  15a  and  10a,  we  then  find 
that  differentials  of  the  matrix  elements  are  given  by 


(14a) 


(14b) 


(14c) 


(15a) 
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(15b) 


Agij 

= 2(9^2  - 

^13 

*922 

' 2(923  ■ 

^12 

*■333 

= 2(9^3  - 

^23 

''♦x* 

= 9i3  A*x  ■ 9 

23 

y *'222 

93^3^)  A(j) 

&9i3 

= ■‘3i2  ^'I’x 

(911 

- 933)A*y 

+ 923A<}> 

and 

*‘’23 

= <933  ■ 922* 

■ '3i2  ‘♦y  ■ 

where 

A^x 

= A4)sin  0Q  cos 

^0 

>1 

-e- 

< 

= A(})sin  0^  sin 

^0 

and 


(15c) 

(15d) 

(15e) 

(15f) 


(15g) 


(15h) 

(15i) 


A(})^  = A4)COS  0Q  (15j) 

A.  PRINCIPAL  AXES  OF  A GRADIENT  FIELD 

Because  a gradient  matrix  is  both  real  and  synunetric, 
it  is  diagonal  in  an  orthogonal  basis  comprised  of  principal 
axes  with  diagonal  elements  or  eigenvalues  X2»  and  X^* 

Since  the  trace  of  a gradient  matrix  vanishes,  + ^2 
X^  = 0.  The  diagonal  matrix  A and  a matrix  R specifying 
rotation  from  a gradiometer  basis  to  the  basis  comprised 
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of  principal  axes  determine  the  matrix  G measured  in  a 
gradiometer  basis;  namely,  A = KGR,  so  that  G ~ RAR.  As 
a result,  two  eigenvalues  and  three  Euler  angles  specifying 
orientation  of  principal  axes  with  respect  to  a gradiometer 
basis  are  a set  of  five  independent  scalars  that  determine 
gradients  measured  in  a gradiometer  basis. 

If  orientation  of  both  the  gradiometer  basis  and  princi- 
pal axes  is  specified  with  respect  to  an  arbitrary  fiducial 

basis  by  matrices  R and  R , respectively,  then  R = R„R, 
^^g  p pg 

and  so  G = R R AR  R . On  the  one  hand,  if  principal  axes 
g p p g 

are  taken  as  the  fiducial  basis,  then  R = I and  G = R AR  , 

p g g 

and  on  the  other  hand,  if  the  gradiometer  basis  is  taken  as 
the  fiducial  basis,  then  R = I and  G = R AR  . For  a gradio- 

g p p 

meter  aligned  along  principal  axes,  R^  = R^  and  so  G = A. 
Transverse  gradients  then  vanish,  and  longitudinal  gradients 
equal  the  eigenvalues  A2»  and  A^. 

Eigenvalues  of  a gradient  matrix  are  roots  of  the  charac- 
teristic equation 

det(G  - AI)  = (A  - A^) (A  - A2) (A  - A3)  = 0 (16a) 

or 

A - ^ (g]^]^  g33  *^11^33  ^ ^12  ^13  g23^~‘^®^  (16b) 

Roots  of  the  characteristic  equation  are  ordered,  and  we 
name  them  so  that  Aj^>A2>A3»  Because  the  sum  of  eigen- 
values vanishes,  at  least  one  root  is  positive  and  one  is 
negative,  so  that  and  A3  are,  respectively,  the  greatest 
positive  root  and  smallest  negative  root  of  the  characteristic 
equation. 


116 


Direction  cosines  of  corresponding  eigenvectors  e.^, 

which  define  directions  of  principal  axes,  form  the  columns 

of  the  matrix  R specifying  rotation  from  a gradiometer  basis 

to  the  basis  comprised  of  principal  axes;  namely,  r = 

X -e  Elements  of  the  kth  column  of  R are  a solution  of 
i k ' 

the  system  of  three  homogenous  equations 


% ' q..r.,  = X,r.,  , i ~ 1/2,3  , 

]k  k ik 

j 

obtained  from  the  matrix  equation  GR  = RA.  T 
tic  equation  assures  existence  of  a solution. 


(17) 


characteris- 


Consequently , we  find  from  Equation  17  that  polar 
angles  and  6^  specifying  the  direction  of  with  respect 
to  the  gradiometer  basis  are  determined  by  the  relations 

tan0j^cos({)j^  = *^13  ^^k  ^22^ ^12^23 

^^k  ■ ^11^  ^^k  ■ ^22^"^12 


and 


tan0j^sin({)j^ 


^23^^k  ~ ^11^  ^12^13 

^^k  “ ^11^  ^^k  “ ^22^“^12 


(18b) 


We  observe  that  polar  angles  (f>j^  = <j>^+Tr  and  0j^  = iT-0p.  as 
well  as  and  3,^  are  solutions  of  Equations  18a  and  18b  for 
each  eigenvalue,  so  that  both  ejt  and  -e^  are  eigenvectors 
corresponding  to  eigenvalue  Xj^,  as  is  also  evident  from 
Equation  17.  Because  the  sense  of  an  eigenvector  is  inde- 
terminant,  we  can  combine  three  eigenvectors  in  eight  ways 
to  form  eight  orthogonal  bases  comorised  of  principal  axes. 
Four  of  the  eight  bases  are  right-handed,  and  four  are  left- 
handed.  We  exclude  left-handed  bases  and  find  that  four 
right-handed,  orthogonal  bases  comorised  of  principal  axes 
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are 


defined  by  the  four  sets  of 


/\ 

eiqenvectors  (e,  , 

. , A A A : 

and  (-e^ , e^ , -e^)  . 


To  determine  a first  set  of  eigenvectors,  we  choose 
e^^  and  e^  so  that  0^6j^^Tr/2  and  O^O^^tt/2  and  require 
right-handedness,  namely  that  62  = X e^^,  so  that 

cos  02  = sin  0^  sin  0^^  sin((|)j^  - (f)^) 


From  Equations  18a  and  18b,  we  then  find  a unique  set  of 
polar  angles  (cj)^,  0^; 


({)_,  0„;  (p-,  0,)  defining  directions 

^ -^A  "^A 

A'  ''  - ' Polar  angles 


of  the  set  of  eigenvectors  (e^^,  62,  . 

defining  directions  of  the  other  three  sets  of  eigenvectors 
are  then  directly  determined.  For  example,  pol^r  angles 

1' 


for  the  set  -^2'  given  by  the  set  (4)^,  0. 


<(,2  + TI 


IT  - 0 


2' 


4,3  + T 


7.-63). 


By  expressing  elements  r^j^  of  the  matrix  R specifying 
rotation  from  the  gradiometer  basis  to  the  basis  defined 
by  the  set  of  eigenvectors  (e^,  02,  33)  in  terms  of  Euler 
angles  <I) , 0,  and  T,  we  find  that  the  Euler  anqles  specify! 
the  rotation  are  determined  by  the  relations 


^)  = 4)^  + ti/2  , 

0 = 63  , 

cost  = sin0  j^sin  (4)^-4)3)  , 


and 


sinT  = sin02sin  (<})2-(1>2)  • 


(19) 


(20a) 

(20b) 

(20c) 

(20d) 
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Euler  anqles  specifying  rotation  from  the  gradiometer  basis 

, A A Ax 

to  bases  defined  by  the  sets  of  eigenvectors  (e^,  -e^,  ~e^)  , 
(-^1,  ^"^1'  ^2'  given  in  terms  of 

those  of  the  first  set  by  the  sets  ($  + tt,tt-0,  tt-'F), 

(<5>,  0,  TT  + H-)  , and  (!>  + ti,  ti  - 0,  2tt  - H')  , respectively. 


As  we  show  subsequently,  eigenvalues  of  a gradient  field 
at  a point  are  equal  to  eigenvalues  of  the  gradient  field 
produced  by  a magnetic  dipole  located  on  a sphere  of  unit 
radius  about  the  field  point.  Because  four  sets  of  eigen- 
vectors correspond  to  the  same  set  of  eigenvalues,  a magneti 
dipole  located  at  any  one  of  four  po-’itions  on  the  unit 
sphere  produces  the  same  magnetic  field  gradients  at  the 
field  point. 


C.  GRADIENT  FIELD  OF  A MAGNETIC  DIPOLE 

We  examine  the  gradient  field  of  a magnetic  dipole  both 
to  illustrate  means  of  describing  gradients  of  magnetic 
fields  in  free  space  and  to  provide  the  basis  for  demonstra- 
ting that  gradients  of  a magnetic  field  at  a point  in  free 
space  are  equal  to  gradients  produced  by  a magnetic  dipole 
located  on  a sphere  of  unit  radius  about  the  point. 


The  gradient  of  the  harmonic  function 


-U  iti-r'  _7 

^ = o...—  , = 4ir  X 10  H/m  , 

4tt  r 

gives,  in  MKS  units,  the  magnetic  field  produced  at  a posi- 
tion "t  by  a magnetic  dipole  of  moment  m.  The  relation 


b = V4; 


Moin 


4tt  r 


3 


2(m.r)r 


r X (m  X r) 


(21a) 


(21b) 
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then,  expresses  the  magnetic  field  b as  a sum  of  comoonents 

. A 

parallel  and  perpendicular  to  the  unit  vector  r,  where  r - 
r'/r  and  m = ni/m.  Equation  21b  tells  us  that  the  magnetic 
field  produced  at  a position  r by  a dipole  of  strength  m 
is  a sum  of  fields  produced  by  one  dipole  of  strength 
m cos  Y directed  along  r and  another  of  strength  m sin  x 

A A 

directed  perpendicular  to  r,  where  cos  x - 

To  examine  gradients  of  the  magnetic  field  of  a dipole, 
we  consider  the  gradient  field  in  three  right-handed,  ortho- 
gonal bases:  (1)  a basis  with  its  x^  axis  alignea 

along  the  dipole  axis  m,  which  we  call  the  dipole  basis; 

(2)  a basis  I y^  I with  its  aligned  along  the  poL  Ition 

vector  r,  which  we  call  the  position-vector  basis;  and 

(3)  a basis  |^^|with  its  axis  aligned  in  an  arbitrary 
direction . 


1 . Dipole  basis 


We  chooTM  a basis  |x^^with  the  x^  axis  directed  along 

the  dipole  axis  m;  the  X2  axis,  along  m X r;  and  the 
axis,  along  {m  X r)  X m.  From  Equation  3 and  the  harmonic 
potential  function_of  a dipole  expressed  by  Equation  21a, 
then,  we  find  that  the  matrix  Gx  representing  gpdients  of 
the  magnetic  field  of  a dipole  in  the  basis  | expressed 

as 


cosxd  - Ssin^x)  0 sinxd  - 5cos"x) 


Gx  = g 


cosx 


. 2 


sinxd  - 5cc3  x)  0 cosx(-2  + 5sin  x^ 


(22a) 
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(22b) 


where  cosx  = m^r  and 


We  then  find,  from  Equation  16b,  that  its  eigenvalues 
are  given  by 

(23a) 
(23b) 


^ 2 


, C 2 ,1/2 

(4  + 5cos  x)  ~ cosx 


= g cosx 


and 


(4  + 


c 2,1/2 

5cos  x) 


I 


+ cosx 


(23c) 


Figure  AF  delineates  ^3/g  functions  of 

X and  shows  that  is  positive,  negative,  and  A^ 

is  positive  for  0<x'^.  ff/2  and  negative  for  •n/2<x'^Ti. 

Moreover,  the  set  of  Euler  angles  (tt/2,  0,  tt/2),  where 
0 = X + 

tan  2a  = (2/3)tanx  (0<><tt/2)  (24a) 

specifies  a rotation  from  the  dipole  basis  to  the  basis 
comprised  of  principal  axes  defined  by  the  set  of  eigen- 
vectors e^,  e^) , as  illustrated  in  Figure  A3.  The 

angle  a between  the  eigenvector  position  vector 

r ranges  from  0 to  tt/2  as  x ranges  from  0 to  it.  The  relation 


tanB  = (l/2)tanx 


(24b) 


dstsrminss  ths  angle  B between  the  magnetic  field  vector  and 
the  position  vector,  which  ranges  from  0 to  and  as  x 
ranges  from  0 to  it  . 

The  sets  of  Euler  angles  (3tt/2,  tt  - 0,  tt/2)  , (tt/2,  0, 

3-rr/2)  , and  (3-rr/2,  tt  - 0 , 3;./2)  specify  rotations  to  bases 
defined  by  the  sets  of  eigenvectors  (e^,  ~e^) , 

63),  and  (-e^,  ^2'  ' respectively.  Figure  A4 

depicts  orientation  of  principal  axes,  marked  by  eigen- 
vectors 63^  and  ^3,  for  several  values  of  the  angle  x* 
Principal  axes  corresponding  to  eigenvalues  and  A3  lie 
in  a meridian  plane  of  the  dipole,  the  plane  defined  by  the 
dipole  axis  and  the  position  vector,  and  the  principal  axis 
corresponding  to  eigenvalue  A 2 is  normal  to  a meridian  plane 

2 . Position-vector  basis 

We  next  choose  an  orthogonal  basis  jy^j  in  which  the 
axis  is  directed  along  r and  orientation  of  the  dipole 
is  specified  by  polar  angles  and  x,  as  shown  in  Figure 
A5.  The  matrix  R(fi,  x)  defined  by 

R(fl,  X)  = + tt/2)  X(x)Z  (tt/2) 

specifies  a rotation  that  brings  the  basis  |y^|  into  coin- 
cidence with  the  basis  | x^  |« 

Consequently,  the  relatio.i 

G = R(fi,  X) 

y ^ 

determines  the  matrix  representing  gradients  of  a dipole 
field  in  the  basis  | y. | in  terms  of  their  representation 


G in  the  basis  |^|  , as  given  by  Equation  22a.  We  then 

5C  ^ 

find  that 


y 


cos  X 0 

0 cos  X 

sin  X cos  Q.  sin  x ^ 


sin  x cos 
sin  X sin 

-2  cos  X 


Expression  26b  tells  us  that  longitudinal  gradients  in 


the  basis  | | result  from  a dipole  of  strength  m cos  x 

directed  along  and  that  transverse  gradients  result  from 
dipoles  of  strength  m sin  x cos  Q and  m sin  x sin  directed 
along  the  y^^  and  ^2  axes,  respectively.  Eigenvalues  and 
corresponding  principal  axes,  of  course,  are  the  same  as 
given  before.  The  set  of  Euler  angles  (JJ  + tt/2,  -a,  ~v/2)  , 
however,  specifies  rotation  from  the  position-vector  basis 
to  the  basis  comprised  of  principal  axes  defined  by  the  set 
of  eigenvectors  62^  ^3)  • 


3 . Arbitrary  basis 

Finally,  we  choose  an  arbitrary  orthogonal  basis  |z^[ 

in  which  orientation  of  the  dipole  is  specified  by  polar 

angles  d and  9 and  direction  of  the  position  vector  is 
^ m m 

specified  by  polar  angles  and  9^.,  as  shown  in  Figure  A6 . 
The  matrix  9^)  defined  by 

0J,)  = Z(<Pr  + tt/2)X(9^) 

specifies  a rotation  that  brings  the  basis  | \ into  coin- 

cidence  with  the  basis  jy^J  defined  in  Figure  A5  provided 

cos  Q sin  x = sin  9j^sin(({)j^  - ip^.) 


(26b) 


(27a) 


(27b) 
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(27c) 


sin  n sin  x sin  “ cose^sin9j^cos  - (|)^) 


and 


cos  X = COS0  cos0_  + sin0_sin0  cos  ((J)^  - 


r m 


r m 


m ’ r 


(27d) 


The  relation 


= R{^^,  0^)R(J^,  X)Gj^R(Q,  x)R(4)j.,  0J.)  / 


(28a) 


then,  determines  the  matrix  G representing  gradients  of  a 

M I 1 

dipole  field  in  an  arbitrary  basis  z . 

(a  ’ ^ 

representation  in  the  basis  I x 
that 


in  terms  of  their 
. As  a result,  we  find 


= g cos  x ji  - 3/2  (c()^,  0^)j 

jcos  A2(4)^/  0J.)  “ sin  0 A^((j)^,  0^)j  , 


(28b) 


+g  sin  X 


where  I denotes  the  unit  matrix. 


Aj^((fj.,  0J.) 


2 2 

2cos  d>  sin  9 
^r  r 

2 

sin2(t)^  sin  0^ 

coscf)  sin20^ 
r r 


. 2 


sin2tt>  sin  0^  cosc})^sin20^ 

2 2.. 

2sin  4>j.sin  0^  sin(J)^sin20^ 


sint})  sin20 
r r 


2 

2cos  0 


r -J 


(28c) 


6^) 


and 


A3(*^,  0^) 


■sin2(t)^sin0j. 
cos2(t)^sin0^ 
[_-sinct)^cos0^ 


cos2(})^sin0^ 

sin2(})^sin0^ 

COS(}>j.COS0^ 


-sin4>^cos0^ 

CCS  6 COS0 

^r  r 


(28d) 


*~cos^(t^sin20j.  (l/2fein24>j.sin20^  coscf)^cos20^ 

(l/2)sin2(i)j.sin20^  sin^(J)^sin20^  sincf)^cos20^ 

cosd)  cos2'6  sin4>„cos20  -sin20 

^ r r r r ^ 


(28e) 
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In  an  arbitrary  basis,  then,  either  of  two  sets  o'- 
four  independent  angles,  0^.;  x)  or  (({)^,  9^;  ' 

and  the  scale  factor  g determine  gradients  of  the  magnetic 
field  of  a dipole.  Because  five  independent  matrix  elements 
determine  gradients  of  a magnetic  field  in  free  space,  gra- 
dients at  a point  in  free  space  are  eguivalent  to  the  gra- 
dient field  of  a magnetic  dipole  located  on  a sphere  of 
unit  radius  about  the  point.  By  requiring  G = for  a 
matrix  representing  gradients  of  a magnetic  field,  v/e  can 
in  principle  determine  a set  of  four  angles  ((p^,  6^;  i})^,  0^^)  , 
which  specify  location  and  orientation  of  a dipole,  and  a 
scale  factor  g by  solving  the  system  of  five  equations  obtained 
by  equating  matrix  elements.  Choosing  a unit  radius  in  the 
scale  factor  g then  determines  the  strength  of  an  equivalent 
dipole.  Nonetheless,  the  direct  approach  soon  becomes  fouled 
by  algebra,  and  we  choose  another  approach  to  obtain  equiva- 
lent dipoles  of  a gradient  field  in  free  space. 


D.  EQUIVALENT  DIPOLES  OF  A GRADIENT  FIELD 


We  seek  location,  orientation,  and  strength  of  a magnetic 
dipole  that  produces  gradients,  represented  by  a matrix  G^, 
at  a point  in  free  space  that  are  the  same  as  gradients 
represented  by  a matrix  G;  namely,  we  ask  that  the  differen- 
tial matrix  AG  vanish,  where  AG  = G - G(j.  If  elements  of  a 
matrix  vanish  in  one  basis,  they  vanish  in  all  bases,  and  so 
we  require  elements  of  the  differential  matrix  to  vanish  in  a 
basis  in  which  G is  diagonal. 

The  relations 


^d 


^d^d^d 


and 


G = RAR 


(29a) 


(29b) 
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express  matrices  and  G in  terms  of  their  respective  dia- 
gonal matrices  and  A and  matrices  and  R that  specify 
i-Qtations  from  a fiducial  basis  to  bases  comprised  of  their 
respective  principal  axes.  Rotation  to  principal  axes  of 
the  dipole  field  is  the  resultant  of  two  consecutive  rotations 
expressed  by  the  matrix  product 

^d  ^^6  ' 

where  R^  represents  a rotation  from  principal  axes  of  the 
gradient  field  represented  by  G to  principal  axes  of  the 
gradient  field  of  the  dipole.  Consequently,  we  express 

the  differential  matrix  AG  as 
AG  = R(A  - R 

Expression  29d  tells  us  that  the  differential  matrix  vanishes 
if  both  the  principal  axes  of  the  two  gradient  fields  coin- 
cide, Rg  = I,  and  their  eigenvalues  are  equal,  A^  = A . 

1 . Dipole  orientation  and  location 

To  determine  orientation  and  location  of  an  equivalent 
dipole,  then,  we  first  find  eigenvalues  ^2'  ^3^ 

gradient  field  represented  by  the  matrix  G and  a corresponding 
set  of  Euler  angles  $ , 6,  and  f that  specify  rotation  from  a 
fiducial  basis  to  principal  axes  defined  by  the  set  of  eigen- 
vector • ^2'  Namely,  we  use  Equations  16b,  18a,  20a, 

20b,  20c,  and  20d  to  find  the  set  of  five  independent  scalars 
(A  , Xj'  from  the  five  independent  elements  ^^11' 

^33'  ^12'  ^13'  ^23^  matrix  G.  We  recall  that  the  set 

of  Euler  angles  specifies  one  of  four  possible  rotations  to 
a basis  comprised  of  principal  axes  because  the  sense  of  an 


(29c) 


(29d) 
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eigenvector  is  indeterminant . 


Because  eigenvalues  of  the  matrix  G are  ordered 
A2  > and  their  sum  vanishes,  we  find  from  Equations 
23a,  23b,  and  23c  that  they  are  equal  to  eigenvalues  of 
the  gradient  field  of  a dipole  at  a position  for  which 


tan  X 


, 0 < X < TT 


(30a) 


and 


g 


4Trr^ 


(30b) 


The  set  of  five  scalars  (g,  x»  ©#  '*')  » then,  specifies  the 
gradient  field  represented  by  the  matrix  G. 

Figure  A7  shows  orientation  and  location  of  the  dipole 
that  produces  a gradient  field  whose  eigenvalues  are 
X^,  X2,  and  X^  and  whose  principal  axes  coincide  with  the 
set  of  eigenvectors  (e^^,  62#  e^)  • The  dipole  axis  and 
position  vector  r^^  lie  in  the  plane  normal  to  the  eigen- 
vector 62*  Equation  30a  determines  the  angle  x between  the 
dipole  axis  and  position  vector.-  and  Equation  30b,  magnitude 
of  the  dipole  moment.  Dashed  lines  in  Figure  A7  delineate 
orientation  and  location  of  dipoles  that  produce  gradient 
fields,  having  eigenvalues  X^^,  X2»  and  X^,  whose  principal 
axes  coincide  with  the  sets  of  eigenvectors  ( e^,  ~^2'  ”^3^ ' 
4^,  -4^),  and  (-4j^,  -e^,  e^)  as  noted.  As  is  evident, 
the  four  equivalent  dipoles  lie  in  the  plane  normal  to  the 
eigenvector  and  they  form  the  two  pairs  (m^,  r^^;  "m^^, 

-F^)  and  (m2,  f2'  with  the  first  pciir  associated 

with  the  eigenvector  62  and  the  second, with  the  eigenvector 

A 

-®2- 
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Wg  obsGrvG  from  FigurG  A7  that  colunin  matricGS 


sin  (a+x) 

M = 0 

P 

cos  (a+x) 


and 


-sina 

0 

-cosa 


( 31) 


givG  dirGCtion  cosines  of  the  dipolG  axis  and  position  vgc 
tor  m and  r",  rGspGctivGly , in  thG  basis  defined  by  the  set 
of  eigenvectors  (e^,  e^)  . Consequently, in  the  fiducial 

basis  uhe  matrices  M = RM^  and  p = Rp^  give  orientation  ana 
location  of  the  dioole  whose  principal  axes,  defined 
by  the  eigenvectors  (e^,  62,  e^) , coincide  with  those  of  the 
gradient  field  represented  by  the  matrix  G.  The  matrix 
R describes  rotation  from  the  fiducial  basis  to  principal 
axes  of  G as  specified  by  the  Euler  and  angles  <I> , T , and  0. 


In  the  fiducial  basis,  we  then  find  that  polar  angles 

4)  and  6 defining  direction  of  the  dipole  rn,  are  determined 
^m  m 

by  the  relations 


cose  = COS0  cos (a+x)  + sinT  sin0  sin(a+x) 
m 

sin0  cos((l)  -<i')  = cosH*  sin(a+x) 
m m 


(32a) 

(32b) 


and 


sin0  sin((j>  -<I>)  = sinH*  cos0  sin  (a+x)  ”sin0  cos  (a+x) 
m m 

and  that  polar  angles  and  6^  defining  direction  of  the 
position  vector  r^  are  determined  by  the  relations 

-cos0^  = COS0  cosa  + sinH*  sin0  sina 

-sin0^cos (c|>j,“'I’)  = cosH*  sina 


(32c) 


(32d) 

(32e) 
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and 


-sin9  sin(<{)  -'!>)  = sin'P  cos0  sina  - sinOcosa  . (32f) 

r r 

Equations  32a  through  32f  give  a prescription  for 

transforming  the  set  of  five  scalars  (g,  x»  to 

the  set  (g,  9^^^)  that  specifies  the  gradient 

field  represented  by  the  matrix  G in  terms  of  the  gradient 

field  of  its  equivalent  dipole  m^^  located  at  position  r^^. 

The  set  (g,  tt+cI)^,  tt-9^,  corresponds  to  the 

equivalent  dipole  -m,  located  at  position  -f, . We  obtain 

''■*■11  ■*■  .10. 
the  set  (q,  d)  , 9 , (1)  , 9 ) corresponding  to  the  equivalent 
^ r r m m 

dipole  in*2  located  at  position  by  adding  tt  to  the  angle 

4'  in  Equations  32a  through  32f.  Finally/  the  set 

(q,  TT+*  , TT-9  , iT+d)  , tt-9  ) corresponds  to  the  equivalent 
' ^ ' r r m m 

dipole  -m2  located  at  “^2* 

2 . Alignment  errors 

A dipole  placed  at  the  position  of  an  equivalent  dipole 
cancells  gradients  at  the  corresponding  field  point.  Because 
positioning  is  imperfect,  however,  errors  in  aligning  a can- 
celling dipole  leave  residual  gradients  at  the  field  point. 

We  describe  residual  gradients  resulting  f 'm  errors  in  align- 
ment in  terms  of  infinitesimal  rotation  c principal  axes 
from  true  alignment  and  infinitesimal  deviation  of  eigen- 
values from  true  values. 

The  matrix 

= I -I-  A4)Q((f)Q,  9q)  (33) 

specifies  an  infinitesimal  rotation  from  principal  axes  of 
the  gradient  field  represented  by  the  matrix  G through  an 
angle  A(t>  about  an  axis  specified  by  polar  angles  (p^  and  9^ 
with  respect  to  the  principal  axes,  where  Q(4>q»  6q) 
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expressed  by  Equatio.i  10a.  Eigenvalues  of  the  gradient  field 
of  the  cancelling  dipole  deviate  from  those  of  the  gradient 
field  represented  by  the  matrix  G by  an  amount  AA^  = 


From  Equation  29d,  we  then  find  that  the  relation 


AG  = R 


-AA^  + A4)(AQ  - QA) 


R 


(34a) 


gives  residual  gradients  resulting  from  alignment  errors. 
The  difference  in  eigenvalues  results  from  deviations  in 
the  scale  factor,  Ag,  and  in  the  angle  x as  expressed  by 
the  relation 


(34b) 


In  the  basis  comprised  of  principal  axes  of  the  matrix 
G,  then,  residual  gradients  are  given  by 


(35a) 


(35b) 


(35c) 


Ag^^2  ~ - (2X^+X2)  Acl)^ 

(35d) 

Agi3  = (X  j^-X^)  A({)2 

(35e) 

and 
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f 


(35f) 


where 


A(f)j^  = A(() 

sind  cos(t) 
o 

(35g) 

2 ~ 

sin0  sind) 
o 

(35h) 

cosG 

o 

(35i) 

and  and  are  eigenvalues  of  G. 
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TABLE  OF  ROTATIONS 
ROTATION  ABOUT  AN  x ^ AXIS 


The  matrix 


Z (4>) 


cos(})  -sine))  0 

sin4)  cost})  0 

0 0 1 


represents  a positive  rotation  through  an  angle  4>  about  the 
axis  of  a basis.  The  relation 

r = z(4>)Gz(4)) 


gives  gradients , rotated  basis  in  terms  of  gra- 


S, 

g . . , in 

the 

initial  basis,  and  so 

11 

+ g 

22^  2^^11 

g22)cos2({)  + 

gj^2sin2({) 

22 

2*^11 

+ g 

22^  “ 2^^11 

g22)cos2(})  - 

gj^2sin2(}) 

33 

^33 

12 

= g^2COs24> 

I^^ll  “ ^22 

) sin2(}) 

13 

= g^^cosc})  + 

g23sin4) 

23 

= g22COS(|)  - 
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ROTATION  ABOUT  AN  x ^ AXIS 


The  matrix 


X(0)  = 


10  0 
0 COS0  -sin0 

0 sin0  COS0 


represents  a positive  rotation  through  an  angle  0 about  the 
x^  axis  of  a basis.  The  relation 


r = X{0)GX{0) 


gives  gradients 
dients , g ■ ■ , in 


in  the  rotated  basis  in  terms  of  gra- 
the  initial  basis,  and  so 


^22 

2^'^22 

^33^  2^^22 

- g^2)cos20 

+ g22sin20 

^33 

2^^22 

^33^  " 2^'^22 

- g22)cos26 

- g22sin20 

^12 

+ 

^13 

= g^^^cosS 

^23 

= g22cos0 

2^^22  ~ ^33 

) sin20 
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ROTATION  ABOUT  AN  x ^ AXIS 


The  matrix 


Y(iP)  = Z {TT/2)X(t|;)  Z(7i/2) 


cosip  0 sirn|» 

0 10 
-sini|»  0 cosijj 


represents  a positive  rotation  through  an  angle  ijj  about  the 
A 

x^  axis  of  a basis.  The  relation 


r = Y(ijj)GY(iJj) 


gives  gradients,  ^ij ' rotated  basis  in  terms  of  gra- 


dients , 

g..,  in  the  initial  basis,  and  so 
^1 1 

^11 

° + 933) 

+ " 933)00521;  - g^^sin2ij> 

^22 

^22 

^33 

2^^11  ^33^ 

■ 933'“"2»  + 

^12 

= g^2<^osi|j  - g23sim|j 

^13 

= gj^^cos2i|;  + j 

<9ll  - 933)sin2«, 

^23 

= g23ccsif;  + 

I! 
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FIGURE  CAPTIONS 


Figure  Al  Three  success ivf  ependent  rotations  through 

the  Eulerian  angles  <t  , 0 , and  F specifying  a rotation 
from  the  basis  to  a basis  . Rotations  are  made  con- 
secutively in  the  following  order:  (1)  rotation  through  the 

angle  4>  about  the  x axis,  represented  by  the  matrix  Z{^>); 

(2)  rotation  through  the  angle  0 about  an  x^  axis  in  the  di 
rection  X y^,  represented  by  the  matrix  X(0);  and  (3)  ro- 
tation through  the  angle  f about  the  axis,  represented  by 
the  matrix  Z(T).  Positive  rotations  are  right  handed. 

Figure  A2  Eigenvalues  A^,  A^,  and  A^  of  the  gradient  field 
of  a magnetic  dipole  having  a unit  gradient  strength,  g,  shown 
as  functions  of  the  polar  angle  x between  the  axis  of  the  di- 
pole and  a position  vector  emanating  from  the  dipole.  Eigen- 
values are  indexed  so  that  Aj^>A2>A^-  The  determinant,  det 
G,  of  a matrix  representing  gradients  at  a position  in  a di- 
pole field  is  positive  for  x > '^/2  and  negative  for  x*'t/2. 

Its  trace  vanishes  and  so  + ^2  ^ ^3  ~ 

Figure  A3  Basis  comprised  of  principal  axes  defined  by  the 
set  of  eigenvectfjrs  (e^,  ^ position  r in  the  field 

of  a magnetic  dipole  m.  Eigenvectors  e^^  and  e^  c.re  in  the 
meridian  plane  normal  to  the  X2  axis,  and  eigen /ector  e2 
points  along  the  negative  $2  axis.  The  angle  a between  eigen- 
vector ^2  position  vector  T ranges  from  0 to  tt/2  as 

the  polar  angle  x ranges  from  0 to  ti  , and  the  angle  B between 
the  magnetic  field  vector  b and  the  position  vector  ranges 
from  0 to  TF . 

Figure  A4  Orientation  of  principal  axes,  marked  by  eigen- 
vectors and  ^3,  depicted  for  several  headings  in  a meri- 
dian plane  of  the  field  of  a magnetic  dipole  m. 
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Figure  A5  Definition  of  the  position-vector 
axis  points  along  the  position  vector  r,  and 
and  X specify  orientation  of  the  dipole  m in 


basis 

polar  angles 
the  basis  • 


Figure  A6  Polar  angles  0^;  4>j,/  specifying  di 

of  the  dipole  axis  and  position  vector  in  a basis 
Polar  angles  4)^^^  and  6^  give  direction  of  the  dipole  m, 
polar  angles  ^ and  0^,  give  direction  of  the  position 


rection 

and 

vector 


r . 


Figure  A7  Location  and  orientation  of  the  magnetic  dipole 
that  produces  a gradient  field  having  principal  axes  coincident 
with  the  set  of  eigenvectors  ^2'  ^3^'  dipole  axis 

and  position  vector  are  in  the  plane  normal  to  the  eigenvector 
^ . Dashed  lines  mark  location  and  orientation  of  image  dipoles 
that  produce  the  same  gradient  field.  They  form  the  two  pairs 
(m^,  r^;  -m\ , -r\)  and  (1^2.  with  the  first  pair 

associated  with  the  eigenvector  ^2  second,  with  the 

eigenvector  ”^2* 


FIGURE  A3 


METRIC  SYSTEM 


BASE  UNITS; 

Quant  itj^  Unit 


length 

metre 

mass 

kilogram 

time 

second 

eiectric  current 

ampere 

thermodynamic  temperature 

kelvin 

amount  of  substance 

mole 

luminous  intensity 

candela 

SUPPLEMENTARY  UNITS; 

plane  angle 

radian 

solid  angle 

steradian 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

activity  (of  a radioactive  source) 

disintegration  per  second 

angular  acceleration 

radian  per  second  squared 

angular  velocity 

radian  per  second 

area 

square  metre 

density 

kilogram  per  cubic  metre 

electric  capacitance 

farad 

electrical  conductance 

siemens 

electric  field  strength 

volt  per  metre 

electric  inductance 

henry 

electric  potential  difference 

volt 

electric  resistance 

ohm 

electromotive  force 

volt 

energy 

joule 

entropy 

joule  per  kelvin 

force 

newton 

frequency 

hertz 

illuminance 

lux 

luminance 

candela  per  square  metre 

luminous  flux 

lumen 

magnetic  field  strength 

ampere  per  metre 

magnetic  flux 

weber 

magnetic  flux  density 

tesla 

magnetomotive  force 

ampere 

power 

watt 

pressure 

pascal 

quantity  of  electricity 

coulomb 

quantity  uf  heat 

joule 

radiant  intensity 

watt  per  steradian 

specific  heat 

joule  per  kilugram-kelvin 

stress 

pascal 

thermal  conductivity 

watt  per  metre-kelvin 

velocity 

metre  per  second 

viscosity,  dynamic 

pascal-second 

viscosity,  kinematic 

square  metre  per  second 

voltage 

volt 

volume 

cubic  metre 

wavenumber 

reciprocal  metre 

work 

joule 

Si  PREFIXES; 

Multiplication  Factors 

1 000  000  000  000  = lO'J 
1 000  000  000  = 10* 

1 000  000=  itr 
1 000-  10’ 
100  - 10* 

10  = 10' 
0.1  = 10-' 
0.01  = 1t)-> 
0.001  - 10-’ 
0.000  001  - 10-‘ 
O.OMt)  000  001  = 10-* 
0.000  tHK)  tKKl  001  ^ 10" 
O.OOO  000  000  (Kit)  001  » 10->’ 
0.000  000  000  tKK)  tKK)  001  « Itr'* 

* To  be  avoided  where  possible. 


